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Chapter 1 


General Introduction 


1.1 Dynamics of Epidemic Diseases 


Epidemiology is a discipline which deals with the study of disease in a population. It has 
been defined as ‘the study of the distribution and determinant of health related states or 
events in specified population, and the application of this study to the control of health 
problems’ (Last 1988). 

It is well known that the spread of various infectious diseases in human populations, is 
one of the most serious problems that modern society faces today. It causes mortality of 
millions of people as well as expenditure of enormous amounts of money in health care 
and control of diseases. Thus it is essential that adequate attention must be paid to the 
study and the control of such diseases. 

It is observed that the transmission of communicable diseases mainly depends upon: (i) 
the susceptible population, (ii) the infective population, (iii) the reservoir population, 
(iv) the mode of transmission. In the following, we provide introductory information 
about transmission of infectious diseases as given by Park (1997). 

Susceptibles are those who are prone to be infected and infectives are those who are 
already infected. The source or reservoir of infection is defined as ‘the person, animal, 
arthropod plant, soil or substances (or any combination of these) in which an infectious 
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perature humidity etc, (iv) ecological factors such as vegetation, biomass density and (v) 
geographical factors. But the study of the spread of an infectious disease using math- 
ematical models by considering these factors has not been done, particularly wh(-n the 
density of the human and carrier, bacteria or vector populations are variables. By con- 
sidering the example of malaria, it may be noted that in highly endemic areas, such as 
in parts of Africa or the north-eastern part of India, persons who have been rejieatedly 
infected with parasites of malaria acquire a degree of immunity which suppresses most 
clinical symptoms. These people may carry gametocytes in their blood that infect the 
mosquitoes biting them and form a separate class of reservoir population which helps in 
spreading malaria without being affected themselves. The effect of the reservoii popula- 
tion on the spread of malaria or any other infectious disease has also not been st iulifnl 
using mathematical models. Both the effects are important and need to be modeled and 
analyzed. 


A common phenomenon is the migration of population within the same region due to eco- 
nomic, social, religious, political, environmental or other considerations. The migrating 
population carries with it all its traditional values, cultural heritage and so on, including 
diseases if any are present in the population. After movement of the population into the 
new habitat the susceptibles join the new susceptible class and the infectives join tlie new 
infective class and the usual interaction begins. The same phenomenon exists between 
two socially structured populations (rich and poor) living in two environmentally differ- 
ent regions, one is cleaner and the other is affected by various discharges caused by t he 
population in that region and a fraction of this population moves to the environmentally 
better region. Further, it may happen that the poor population living in the environ- 
mentally degraded habitat spread the disease by interacting with the rich population by 
working in their habitat as service providers (without migration). All these phenomena 
need to be modeled and analyzed. 


In the modelling of infectious diseases, as suggested above, one can use simple mass action 
incidence for interaction of susceptibles with infectives or the true mass action incidence. 
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It can also be proved that, this equilibrium point is globally asymptotically siab.e by 
using the following Liapunov function about P*{Y,N) as follows, 


V=iY-Y-Ylnl)+lh {N-Nf, 


(2.5) 


where h is to be suitably chosen. From (2.3) and (2.5), we get 

fN N\ 

V = {Y -Y)[P{N - N) - P{Y -Y) + XCm^Y yj^ 

+h{N-N)[-d{N-N)-a{Y-Y)] 

= - L + (Y - Yf - hi(N + '■ ^’)('’ - n- 

Choosing h = it is seen that 1/ is negative definite, implying glol.al stal.ility 

of this equilibrium point P*. The above theorem implies that the s}>re;ui oi a earner- 
dependent infectious disease increases as the carrier population density incieases due to 
household discharges. 


2.2.2 Case II: Q = Qo + IN 

In this case, using X -f 7 = iV, the model (2.1) is rewritten in the following form: 

7 = /3(iV-7)7 + A(iV-7)C-(i^-l-a-l-d)7, 

N = A-dN-aY, 

C = sC{l-y)-5C + siEC, (2.6) 

E = QiN)-SoE = Qo + lN-5oE, 

X(0) = Xo > 0, 7(0) = 7o > 0, N{0) = No >0, C{0) = Co > 0 and E{0) = Eo > 0. 

Now we discuss the equilibrium analysis of the above system and results of tlu^ analysis 

are summarized in the following theorem. 

% 

THEOREM 2.S There exist the following three eguilibrio^ notncly 

(*)P^ (o. i 0. ®#). 

(ii) P 2 (Y, N, 0, , which exists if p^-[u + a + d) > 0, 



Acknowledgement 


I express my sincere gratitude to Professor Peeyush Chandra and Professor Prawal 
Sinha, for their guidance and support over the years. I am extremely grateful to them 
for introducing me to epidemiological modelling. 

Special thanks are due to Professor J. B. Shukla who had always time to listen and 
to give me his valuable comments and suggestions. I thank him for his help. 

I am indebted to all the faculty members of the department of Mathematics, IIT 
Kanpur, particularly to Professor V. Raghavendra for helping me in various matters. 

I thank the staff members of the Department of Mathematics, IIT Kanpur for their 
help during this period. 

I extend my thanks to Professors at S. M. College, Bhagalpur, Professor Aruna Roy 
and Professor Mill Ganguly for encouraging and inspiring me to pursue higher study in 
Mathematics. 

I owe my thanks to all my colleagues especially Bishwaranjan, Dipak, Kailash, Manish, 
Muthu, Parasar, S. B. Rao and Sajid for their enthusiastic help throughout my stay at 
IIT Kanpur. I acknowledge Amrapali, Chaitali, Kusum, Malabika, Monideepa, Reeta 
and Shruti for helping me in various matters. My friends Sushmita and Swagata have 
been always there whenever I needed them for help and advice. I remain indebted for 
their love, encouragement and support. 

At last, but by no means the least, I wish to express my deep gratitude to my parents, 
brothers Prashant and Rajiv and sister Mili, whose love and support have been an in- 
valuable source of inspiration to me. Tlunr faith in my abilities and their encouragement 
have been instrumental throughout this period. 


Mini Ghosh 



30 


Modelling the Spread of Carrier- Dependent Infectious Diseases with Environmental and 

Demographic Effects 



Figure 2.9: Existence of equilibrium point. 

where Cm = 7(5 — ^ 4 - Si which increases as the household di.scharge raft* Qa 
increases. 

The result of an equilibrium analysis is stated in the following theorem. 


THEOREM 2.5 There exist the following two equilibria, namely (i) £’t((J,0) and 
( a) E 2 (Y, N), which exists ifv + a + d> ^ XCm- 


Proof: Existence of E^ is obvious. The existence of second equilibrium point i.s .shown 
as follows: 

Setting the right hand side of (2.11) to zero, we get 


r N 
Y=-(l--)N 


PY^ - [{P - (1 - a)^}N -{u + a + d + XCm)]Y ~ XCmN 


0 . 


( 2 . 12 ) 

(2.13) 


It may be pointed out that in the N-Y plane (2.12) gives a parabola with vertex (^: -ti^) 


aiid passing through the points (0,0) and (iT.O). while (2.13) gives a hyperbola with a 

branch m the first and fourth quadrants and passing through (0,0). The two curves will 

rntersect at a point (K,iV) provided the slope of parabola at (0,0) is more than that of 

the hyperbola branch in first quadrant at (0,0), he. slope of (2.13) at (0,0) is le.,., than 
that of slope of (2.12) at (0,0), 


i' + a + d > IaO. 


r 


(2.14) 



Synopsis 


It is well known that the spread of various infectious diseases in human populations, is 
one of the most serious problems that modern society faces today as it causes mortality 
of millions of people, and costs enormous amount of money in health care and control 
of diseases. In general the spread of an infectious disease is governed by factors such 
as : (i) the density of the human population, (ii) the density of the carrier, bacteria 
or vector population, (iii) various environmental factors such as household discharges, 
rain, temperature, humidity etc, (iv) ecological factors such as vegetation and biomass, 
(v) geographical factors. It is noted here that the study of the spread of infectious 
diseases using mathematical models by considering some of these factors has not been 
done, particularly when the density of the carrier, bacteria and/or vector populations are 
variables. For example, household discharges may be very conducive to the growth of 
the carrier population, such as flies, cockroaches as well as the bacterial population and 
the vector population (such as mosquitoes). This increase helps in enhancing the contact 
rate of infectives with susceptibles, leading to fast spread of the disease. 

Keeping the above in view, we have modeled the spread of infectious diseases by using 
simple mass action incidence. The effects of indirect contact with carriers, bacteria or 
vectors have been included in the models proposed. Two types of human demography 
have been considered, (i) constant immigration and (ii) logistic growth. The thesis con- 
sists of seven chapters and is divided into two parts. The first part deals with modelling 
infectious diseases in a single population whereas the second part deals with the spread of 
infectious diseases into two populations living in two different environmental conditions. 

The first chapter provides the introduction regarding the dynamics of infectious diseases 
including a brief literature survey relevant to modelling and analysis of infectious diseases, 
so that the work done in the thesis can be seen in its proper perspective. 

In the second chapter, we study the effects of the following factors on the spread of a 
carrier-dependent infectious disease using SIS mathematical models: 
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Figure 2.13: Variation of infective population with susceptible poiiulation. 



Figure 2.14: Variation of infective population with time for different cumulative environ- 
mental discharge rates. 
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not satisfied simultaneously, so tliis equilibrium is 


We note that both the conditions are 

unstable. - -■ . . , 

ix M 2 at equilibrium point BtiY, N, ) Rio - 

, AgN\ fly + AB A(JV - Y] 

^ -d ,0 ,, 

0 -(fB + ‘’t^) 


The variational matrix 

M2 = 


(-{0Y . y 

-a 
s\ 


V 

The characteristic polynomial is 

•0^ + + <22^ + “3 = 0, 


where' 


ai 


= /3y + — ;^ + (i + T^ + ' 


Y 


B 




B 


a, = <i 1 (31' + -^ j + T J '* '■' ) ^ ' 


B , 


Here ai > 0, as > 0 and also 0102 - 03 is 
stable if it exists. 


positive. HeiK:(‘ £’•' hnnilly asymptotically 


Nonlinear Analysis and Simulation: 

Before proceeding to computer simulation, we first analyze the nmdel (.b,l) vlun the 
disease related death rate a is zero. In this case the model (:i.:5) reduces l.o the loliowing, 


Y = {^Y + XB){N -Y)-{v + d)Y, 

N = A-dN, 

B = sb{i-j^+siY-SoB + 5B^. 


(3.6) 


Equating the right hand sides of the above equations to zero we get the following two 
equations for finding equilibrium points: 


I3Y^ -(p^--v-Xb]y- XB 


A 


sB +SiY-soB-^SB 


d 

Qa 


(3.7) 

(3.8) 



VI 


(i) variable carrier population, the density of which follows a logistic model and increases 
further due to household discharges in the habitat, 

(ii) human demography. 

The following two types of demographic factors are considered in the analysis . 

(i) a population with constant immigration and death rate, 

(ii) a population with logistic growth. 

It is shown that the endemic equilibrium is globally stable when the rate of cumulative 
environmental discharges, conducive to the growth of carrier population, is a constant. 
When the rate is a function of human population density, the endemic equilibria is 
also locally and globally stable but under certain conditions. The global stability is 
also demonstrated by computer simulation. It is shown that if the human population 
increases due to demographic changes, the spread of the infectious disease increases. It 
is further noted from the analysis that the spread of the infectious disease increases as 
the growth rate of the carrier population (caused by conducive environmental discharges 
due to human sources) increases. 

In the third chapter, we have considered a bacterial disease model, which is very similar 
to the models of Chapter 2 involving carriers. Here the only difference is that bacteria 
is also released by infectives continuously into the habitats and its growth rate increases 
not only because of household discharges but also by continuous release from infectives. 
Here-the SIS- models for bacterial infectious diseases, (like T.B., typhoid, etc.), which 
are caused by direct contact of susceptibles with infectives and also by bacteria, are 
proposed and analyzed quantitatively as well as by simulation. It is found again that the 
endemic equilibria is globally stable when the rate of cumulative environmental discharges 
conducive to the growth of bacterial population is a constant. But when the rate of 
discharge is a function of human population density, the endemic equilibrium is locally 
and globally stable only under certain conditions. It is further shown that under various 
demographic situations mentioned above, the infective population increases with the 
increase of household discharges and the spread of the disease increases and it becomes 
more endemic. 
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Figure 3.9: Variation of infective population with susceptible population. 



Figure 3.10: Variation of infective population with time for different intrinsic growth rate 
of bacteria population. 
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In Chapters 4 and 5, malaria models without and with reservoir have been proposed and 
analyzed by considering, 

(i) the effect of cumulative discharges from household sources, and (ii) the effects of de- 
mographics of the human and mosquito populations. 

The cases of constant as well as variable discharge rates are considered. The threshold pa- 
rameter for spread of the diseases is derived in each case. If this threshold is greater than 
one, the corresponding nontrivial equilibrium is always feasible and is locally asymptoti- 
cally stable to small perturbations. It is also shown by analysis and simulation that this 
equilibrium is globally stable under some conditions. The effect of demographic increase 
in the human population is found to increase the spread of malaria. It is observed that 
as the cumulative environmental discharges increase, the mosquito population increases, 
leading to fast spread of malaria. 

In the Chapter 6, a model is considered to study the effect of immigration of population 
from an environmentally degraded habitat to an environmentally cleaner habitat on the 
spread of an infectious disease. Constant as well as variable household discharge rates 
are considered. In each case the conditions for existence of endemic equilibrium points 
have been obtained. It has been shown that under certain conditions, these equilibria are 
locally stable indicating that the disease may become endemic. Further by simulation, 
it is observed that the endemic equilibrium points are in fact globally stable under local 
stability conditions. It is concluded that as the migration rate increases, the total and 
the infective population densities of the region with cleaner environmental conditions 
increase, where as the total and the infective population densities of the region of degraded 
environmental condition decrease as expected. 

Lastly in Chapter 7 a mathematical model is proposed by keeping in view that a poor 
population living in an environmentally degraded habitat spreads the disease by inter- 
acting with the rich population by working in their habitat as service providers (without 
migration) . The model is presented under similar environmental conditions as in Chapter 
6 and as in the previous chapters, here also linear and nonlinear analysis of the equilibria 
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and also, 


{s-so + 5^)[pK-{l-a)r-{u + a + d)] > XKsi- 


( 3 . 26 ) 


Clearly both conditions are not satisfied simultaneously, so this equilibrium is unstable. 

Now the variational matrix M4 at equilibrium point P4 is given by 

/-(/ 3 y+^) ^f + A.B-(l-a)f X{N-Y) ^ 

M 4 = 


—O' 


r — 


2rN 

K 


V 


0 


Si 0 

The characteristic polynomial corresponding to the above matrix is 


+ CLi'lp'^ + CL2'4^ + — O 7 


where 


ai - PY + ^ + ^i2N-K) + ^B+"-^, 


<^2 


WY + ^) 


^{2N - K) + jB + 


+ -{2N-K) 


L B 


+a\^Y + \B-{ 1 - a) 


rY' 


03 = [^Y + 


K ^ 

XBN\fs^ siF' 


xs^{N-n 


+ a\pY + XB-il-a) 


K 


-{2N - K) 

L B y 


— (2iV-p:)Asi(N - Y). 


Hence by the Routh-Hurwitz criteria, the system is locally stable if ai > 0 , 03 > 0 
and ai 02-03 > 0 . However since oi > 0 , the system is stable if 03 > 0 and ai - 03 > 0 
and it is unstable if 03 < 0 or Oi 02 ~ 03 < 0. 

Remark; It is easily seen that the local stability conditions are satisfied if iV > 
Hence for numerical purposes, the set of parameters are chosen in such a way that 
iV > f. 


Nonlinesu- Analysis and Simulation: 

We first analyze the model ( 3 . 18 ) for a = 0 . As before it can be shown that the nontrivial 
equilibrium point of the model ( 3 . 18 ) is globally stable by taking the following Liapunov 
function, 

y = i(y - Yf +pi(Ar - JV - Nln^) + \p^{B - B)\ 

iV ^ 



have been studied for various caaes. It is seen that as the interaction between rich and 
poor classes increases, the disease spreads more rapidly in the rich class as expe 

It is suggested that it is in the interest of the rich class to not only clean their own 
environment but also to clean the environment of the districts in their neighborhood, 

where poor people live. 
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+ 

<^0 

+ 

+ 


f2rN 

V K 


— r 


B . 


I 

2rN 

K 


a + A5 - (1 - a)^Y | - .SiA(iV - Y) 
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By Murata (1977), conditions for local stability of the system arc 



as ai 

1 ^2 


as ai 0 


as > 0, 

> 0, 

1 as ao 

0 as ai 

> 0, 


03 «i 0 0 

1 «■> «() 0 

0 0.3 Oi 0 

0 1 02 Oo 


> 0 . 


First inequality is obviously true, so if other inequalities are satisfied, then this equilib- 
rium is locally asymptotically stable. 

Remark: We note that the second condition is satisfied for .V > So for simulation 
we choose set of parameters such that N > ^. 


Nonlinear Analysis and Simulation: 

As before, we speculate that the system (3.27) may be globally stable for A > f . To 
show this, the system (3.27) is integrated by the fourth order Runge-Kutta method using 
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Modelling the Spread 

The result ot an equilibrium analysis is stated in the following theorem. The proof is 
obvious. 

THEOREM 4.1 There exist following two equilibria, namely 
(i) Eo ( 0 , 0 ) and (ii) Ei{Yi, Ni, Y 2 ), where 

y . iV. = d^. 

and , -a , , 

. 0^024;N2 - 1^2 + + (1 - (^1 + Ofi + di) 

“ A'[Ai + {«2 + <i2 + (i-o')ft:^e} (1 + S )1 

El exists if 

A A d jvj 

A— - = Ro (say) > 1. (4.4) 

| o ;2 + ^2 + (1 - (^^1 + 0^1 + <^ 1,1 

It is easy to see that Rq is the threshold parameter of the system (4-3 )■ 

Remark: We note here that % increases with N 2 - Hence Yi increases as Qa or K 2 
increases, i.e. the infective human population increases as the cuinulative rate of envi- 
ronmental discharges increases. 

4.2. 1.1 Stability Analysis 

In the following, we study the linear analysis of these equilibria and nonlinear analysis of 
the nontrivial equilibrium Ei. We state the local stability of these two equilibria in the 
following theorem. 


THEOREM 4.2 (i) The equilibrium point Eq is locally asymptotically stable if Rq < I, 
otherwise if Rq > 1 it is unstable and then the second equilibrium Ei exists. 

(ii) The second equilibrium Ei, if it exists, is locally asymptotically stable. 
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Figure 4.4; Varisttion of infective population Yi with time foi diffeieiit giowt.h rate coef- 
ficients of the mosquito population. 



Figure 4.5: Variation of infective population Yj with time for different growth rate coef- 
ficients of the mosquito population due to environmental discharges. 




CONTENTS 


XI 


5 Modelling the Spread of Malaria with Human Reservoir: Environ- men- 
tal and Demographic Effects 120 

5.1 Introduction 120 

5.2 Malaria Model with Constant Immigration 121 

5.2.1 Case I : Q = Qa, a Constant 123 

5.2.2 Case 11: Q = Qo + lNi 128 

5.3 Malaria Model with Logistic Population Growth 138 

5.3.1 Case I; Q = Qa, a Constant 139 

5.3.2 Case 11: Q = Qo + lNi 146 

5.4 Conclusions 157 

6 Modelling the Spread of a Carrier- Dependent Infectious Disease in 

Two Neighboring Habitats with Migration in Between 158 

6.1 Introduction 158 

6.2 Mathematical Model 159 

6.3 Case I: Q is a Constant Qa 161 

6.3.1 Stability Analysis 163 

6.4 Case II: Q is a Variable 169 

6.4.1 Stability Analysis 172 

6.5 Conclusions 179 

7 Modelling the Spread of Bacterial Disease in a Population: Effect of 

Service Providers from an Environmentally Degraded Region 181 

7.1 Introduction 181 

7.2 SIS Model with Immigration 182 

7.2.1 Case I: Q is a Constant Qa 184 

7.2.2 Case II: Q is a Function of iV 2 193 



110 


Modelling the Spread of Malaria: Environmental and Demographic Effects 


For uniqueness of this equilibrium we need an additional sufficient condition > o, 

for all Ni > 0 , where Yi is given by 


PII 32 NI N 2 — 

k + d2 + (1 - a')^ ^^ 2 ] {ni + oq + di + (1 a)^} 

y.- j 

with W 2 = ^ |r 2 — a 2 + ( 

^2 [I'l + CKi + di + (1 - a) 
12^}. 

iiAi ' 
Ki J 




Proof: The existence of Ei, E 2 , E 3 or E 4 is obvious. We prove here the existence of E^. 
Setting right hand side of (4.18) to zero, we get the following equations, when Ni ^ 0, 
Ni ^ Ki, N2 ^ 0 , iVs # N2, N2 # N^, 


E = 






5o r2 

+ oi + di + (1 — a)^Ni 


Y, 


Ol Ai 


(Al) 

(A2) 

(4.19) 


_ P1P2N1 N2 - [02 + ^2 + (1 — iV2]{ni + Oi + di + (1 - a)-);^} ^ 

1 - [/32{i/i + + di + (1 - a)^} + ^ ^ 

where N 2 is given in equation (Al). 

Now we use the isocline method for finding the equilibrium point Er, from (4.19) and 
(4.20). Clearly (4.19) is a parabola in the Ni - Yi plane, passing through (0,0) and 
{Ki, 0) with vertex (^, ^). 

(i) From (4.20), we get when Ni — 0, 

^ ^ (ni + oi + di) 




02 + d 2 + (1 - a') |r 2 - 02 + 


P2 (r'l +0:1 + di) + ffiP2 1^2 — Q!2 + 

, , ^2 1 do J 

(ii) Also when Ni = Ki, equation (4.20) gives 


< 0. 


provided 


- |a2 + ^2 + (1 - |i/i +o;^ + dj + (1 _ a)^Ki^ 

^2 [^iAf2 + {r'l + Oil + di + (1 — a)ri}] 


> 0 , 


ffil 32 KiK 2 f , ,Qq + IK, 

ir2-a2 + d2- 


r2 


^0 


|a2 + d2 + (1 ® + oil + di + (1 — a)ri} 


^ 1. 


(4.21) 
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Chapter 5 


Modelling the Spread of Malaria 
with Human Reservoir: Environ- 
mental and Demographic Effects 


5.1 Introduction 


The parasites of malaria are normally transmitted from one person to another via female 
mosquitoes of the Anopheles species. The parasite which infects mosquitoes, one stage of 
which is called gametocyte, develops in capillaries of the inner organs of infected persons 
after the invasion of the blood by merozoites. Mature gametocytes, which are infective to 
mosquitoes, appear in the peripheral blood some three days later in the case of P. vivax 
and after about ten days later in the case of P. falciparum and P. rnalariae. The female 
Anopheles mosquito ingests malaria gametocytes when it takes a blood meal from an 
infected person. Once a parasite enters into a mosquito, it needs a period of development 
before it can affect another person again. The length of this period (the sporogonic cycle) 
depends on the Plasmodium species and the ambient temperature, which is about two 
weeks. 

In highly endemic areas, such as in parts of Africa, persons who have been repeatedly 
infected by malaria acquire a degree of immunity which suppresses most clinical symp- 
toms. These people carry gametocytes in their blood that infect the mosquitoes biting 
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From ( 5 . 8 ), we note the following, 

(i) f ( 0 ) = ^i(/?2 -I- A2^)^A(2* - {0:2 4 - c?2 + (1 - o!)^N2}{yi +(y.i + di + ( 5 i), 

(iii) F(oo) > 0. 

So if F(0) > 0, then there exists one root say Yi between 0 and and another 

root between and 00. Corresponding to this Fi, the values of iVi, F21 N2 and E 

are determined using (5.5) and (5.6) and we get a nontrivial equilibrium point Ez under 
following condition 


{vi -|- a;i -h di -h 5i){a2 -I- 0^2 + (1 “ 


= F;(say) 


> 1 . 


As in previous cases here also i?' is the threshold of the system ( 5 . 4 ). Using ( 5 . 5 ) and 
(5.6), values of Zi, Ni, N2, E can be derived corresponding to Fi. 


5.2.2. 1 Stability Analysis 


We state the stability result of these equilibria Pi, P2 and P3 in the following theorem. 


theorem 5.4 The equilibrium Pi is unstable, the equilibrium P2 is stable if P' < 1, 
otherwise if P' > 1 it is unstable and the equilibrium Pz exists, which is stable under 
the following conditions: 
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‘OJhere 04, 03, 02, Oi and oq are given explicitly in the proof of the theorem. 


Proof: The variational matrix at (Fi, 2 'i, JV'i,F2, A2,P) is given by 
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Proof: The existence of the first two equilibria is obvious. We prove the existence of the 
third nontrivial equilibrium by the isocline method as follows : 

Setting the right hand sides of (5.10) to zero, we get 
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From (5.13) and (5.14), we get (for Yi ^ 0) the following: 
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We note from (5,15) that for a given Ni, we have only one value of Ti- 
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Proof: The variational matrix M at the equilibrium point (Fi, Zi,Ni, Y 2 , N 2 , E) is given 

by 

/ mil ^^13 ^14 0 ^ \ 
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E 2 (0, 0, 0, 0, ^{ra - 02 + <^2^ }, and £;3(0, 0, 0, 0, have one positive char- 

acteristic root, thus implying instability of these equilibria. 


The variational matrix M 4 corresponding to the equilibrium point E 4 is given by 
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Clearly three roots of the above matrix 
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the other roots are given by following cubic equation 
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which is also affected by various household discharges. In the degraded region the effect 
of the carrier population in the model is considered along the same lines as described in 
Chapter 2. Thus a mathematical model is proposed as follows (Fig. 6.1), 

= Ai-diXi-^iX{{i + viYx+dX2, 

Yi = ^xX\Yi — {yi -h Q!i + di)Yi -f 6Y2i 
Ni = Ai — diNi — aiYi + 9N2, 

X2 = A2 — 1^2X2 — {P2Y2 + X 2 C)X 2 — 9X2 + ^ 2 Y 2 i 

Y 2 = (/ 52 I 2 -f A2C')X2 - ( 1^2 + 0^2 + ~ (6-1) 

JV2 = A2 — ^2X2 ~ 0:212 ~ 9N2, 

C = sC — SflC "t" SiEC, 

E = Q{N2)-5oE, 
where s > Sq, ^ > 0, 

Xi( 0 )=Xio > 0 , Yi( 0 ) = rxo> 0 , X2(0) = X20 > 0 , F2(0) = r20>0, 

(7(0) = Co > 0, E(0) =Eo > 0. 

Here, C is the density of the carrier population and E is the cumulative concentration of 
household discharges which are conducive to the growth of the carrier population; Ai and 
A 2 are the constant rates of recruitment/immigration into the respective populations of 
susceptible class; di and ^2 are the natural death rates of the richer and the poorer classes 
respectively; o;i and a 2 are the disease related death rates of the richer and the poorer 
classes respectively; Pi and V 2 are the recovery rates of richer and poor populations; 9 
is the constant rate of migration of people from the environmentally degraded region to 
the cleaner region of the habitat; and P 2 are the rates of contacts of infectives with 
susceptibles; A2 is the contact rate of poor class susceptible with the carriers. As before, 
the cumulative rate of environmental discharges Q{N 2 ) is such that 5^ > 0. We take 
Q{N 2 ) = Qo + IN 2 , where I is a constant. 

In writing the model (6.1), we have assumed that the susceptibles and infectives from 
the degraded environment join the susceptibles and infectives of environmentally better 
habitat and the rate of immigration is proportional to the respective density of the 
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The result of equilibrium analysis is stated in the following theorem. 


THEOREM 6.3 There exist the following four equilibria, 
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disease in both the populations. The later result is also confirmed by computer simu- 
lation. In case II, the condition for existence of nontrivial equilibrium point has been 
obtained. It has been shown that this equilibrium point is locally asymptotically stable 
under certain conditions. By simulation it is shown that the nontrivial equilibrium is also 
globally asymptotically stable under local stability conditions for the set of parameters 
considered. In both the above cases, it is seen that as the migration rate 6 increases, 
the total population and the infective class of the region of better environmental condi- 
tions increase, whereas the total population and the infective population of the region 
of degraded environmental condition decreeise as expected. The analysis suggests that it 
is in the interest of rich population to clean the habitat of the poor population to avoid 
immigration on environmental consideration. 
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Figure 7.2: Variation of Ni with Yi. 


bounded by the solution of (7.3) with af = 0 , i — 1, 2 . So we speculate that the equi- 
librium Ez of (7.3) may be globally stable, provided we start away from other equilibria. 
To illustrate this and to see the effects of various parameters on Yi, the system (7.3) is 
integrated using the fourth order Runge-Kutta Method and using the following parame- 
ter values in the simulation (Greenhalgh 1990, 1992), 

A = 0.00000051, Ai = 0.000000011, = 0.012, di = 0.0004, 

ai = 0.0005, = 10, 132 = 0.00000051, A 2 = 0.0000000002, 

a 2 = 0.00052, V 2 = 0.011, ^2 = 0.0004, A 2 = 10, s = 1 , S 20 = 0.65, 

S 2 = 10, S 3 = 0.000002, L = 5000000, = 20, 60 = 0.001. 

The equilibrium values of Yi, Ni, Y 2 , N 2 and B 2 are found as 
Yi = 741.710, Ni = 24072.893, Y 2 = 2893.195, N 2 = 21238.835, B 2 = 2021558.647. 
Here all parameters are in per day except L and I which are constants. Simulation is 
performed for different initial positions 1 , 2 , 3 and 4 as shown in Figs. 7.2 and 7 . 3 . In 
these figures, we have plotted the infective population against the total population of 
the respective classes. From the solution curves, we conclude that it is plausible that 
for the set of parameters considered and provided that we do not start at Ei or E 2 the 
system tends to the endemic equilibrium (Tj, Ni, iVj, . 82 ). The effects of various 
parameters on Yi and Y 2 are shown in Figs. 7.4-7.7. The following is concluded from 
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General Introduction 


agent lives and multiplies, on which it depends primarily for survival, and where it 
reproduces itself in such a manner that it can be transmitted to a susceptible host 
(Benenson 1981, Park 1997). 

The reservoir may be of three types : 

1. Human reservoir, 2. Animal reservoir, and 3. Reservoir in non-living things. 

All communicable diseases are transmitted from the reservoir or source of infection to the 
susceptible host. They may be transmitted in many different ways, depending upon the 
infectious agent, portal of entry, environmental and ecological conditions. The mode of 
transmission of infectious disease may be direct or indirect as described below in Table 1 
and Table 2 (Last 1988, Benenson 1981). 


Table 1. Direct Transmission 


Type 

Pathway for transmission 

Diseases 

Direct contact 

Direct contact from skin to 
skin, mucosa to mucosa 
or mucosa to skin of 
the same or another person, 
sexual contacts 

STD and AIDS, Lepro.sy, 

Skin and Eye Infection 

Droplet 

infection 

1 

Direct projection of a spray 
of droplets of saliva and 
naso-pharyngeal secretions 
during coughing, sneezing, 
spitting or talking 
into the surrounding 
atmosphere 

Respiratory Infection 
Eruptive Fevers, Common Cold, 
Many Infections of the Nervous 
System, Diphtheria, 
Whooping Cough, Tuberculosis, 
Meningococcal, 
Meningitis, etc. 

Contact with 
soil 

Direct exposure of susceptible 
tissue to the. disease 
agent in soil, compost or 
decaying vegetable matter 

Hookworm Larvae, 
Tetanus, Mycosis, etc. 

Inoculation into 
skin or mucosa 

“T* 1 i 

Disease agent may be 
inoculated directly into 
the skin or mucosa 

Rabies, Hepatitis B. 

Transplacental 
or vertical 
transmission 

Disease agents can be 
transmitted transplacentally 
from mother to child 

AIDS, Syphilis, Rubella, 
Measles, etc. 
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Figure 7.11: Variation of Vi and Yi with time for different rates of release of bacteria 
from the infective population. 



Figure 7.12. Variation of Y\ and Y 2 with time for different rates of growth of bacteria 
population due to the environmental discharges. 
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Table 2. Indirect Transmission 


Type 

Pathway for transmission 

Diseases 

Vehicle-borne 

Transmission of the infectious 
agent through the agency of 
water, food, ice, blood, serum, 
plasma or other biological 
products such as 
tissues and organs 

Acute Diarrheas, 
Typhoid Fever, Cholera, 
Polio, Hepatitis A, 

Food Poisoning 
Intestinal Parasites, etc. 

Vector- borne 

Vectors transmit infection 
by inoculation into the 
skin or mucosa by biting 
or by deposit of infective 
material on the skin or on 
the food or other objects 

Diarrhea, Dysentery, 
Typhoid Fever by House 
Fly, Malaria 

Air-borne transmission 
(a) Droplet nuclei 

Small droplets evaporate 
rapidly leaving behind a 
minute residue or nucleus 
which may be a virus particle 
or bacteria and being very 
light they may remain air- 
borne for considerable periods 

Influenza, Chicken-pox, 
Measles, Tuberculosis 
Typhoid Fever, etc. 

(b) Infected Dust 

When a person coughs or 
sneezes, larger drops of 
moisture settle on the floor 
and become part of the dust, 
so when the dust is inhaled, 
the person acquires infection 

Tuberculosis, 
Pneumonia, etc. 

Fomite-bome 

Soiled clothes, towels, door 
handles, chains, syringes, 
instruments and surgical 
dressings 

Diphtheria, Typhoid Fever, 
Bacillary Dysentery, 
Hepatitis A, Eye 
and Skin Infection 

Unclean hands and 
fingers 

Hands are the most 
common medium by which 
pathogenic agents are 
transferred to food from the 
skin, nose, towels, etc. as 
well as from other foods 

Streptococcal Infections, . 
Typhoid Fever, Dysentery, 
Hepatitis 
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General Introduction 


1.1.1 Stages for Parasitism in a Susceptible Host 

There are four stages in successful parasitism (Last 1988) . 

(a) The infectious agent must find a portals of entry by which it may enter the host. 
There are many portals of entry, e.g. the respiratory tract, the alimentary tract, tlie 
genito-urinary tract, skin etc. Some organisms may have more than one portal of entry, 
e.g. hepatitis B, Q. fever. 

(b) On gaining entry into the host, the organism must reach the appropriate tissue in 
the body of the host where it may find optimum conditions for its multiplication arifl 
survival. 

(c) The disease agent must find a way out of the body (portal of exit) so that it may 
reach a new host and propagate its species. If there is no portal of exit, the infection 
becomes a dead end infection as in rabies, bubonic plague, tetanus, etc. 

(d) After leaving the human body, the organism must survive in the external environment 
for a sufficient period until a new host is found. In addition, a successful disease ag(mt 
should not cause the death of the host but produce only a low-grade immunity so that the 
host is vulnerable again and again to the same infection. The best example is coinnion 
cold virus. 

1.1.2 Effects of Demographic, Environmental and Ecological 
Changes on the Spread of Diseases 


One of the most important factors responsible for the spread of diseases is the demo- 
graphic changes that have occurred in our society on a global basis. Population growth 
in the past 50 years, since the end of World War II, has been phenomenal, especially 
in tropical developing countries. Cities have just expanded in an unplanned fashion, 
causing a deterioration in the housing, water, sewage and waste management systems. 
This creates ideal conditions for the spread of various infectious diseases (e.g. water 
borne, mosquito-borne, rodent-borne diseases) in these urban areas. Further household 
discharges dumped in open space on the road, particularly in areas of the city where poor 
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people live also contribute to the growth of carriers, bacteria and vectors, causing the 
spread of infectious diseases. These factors have contributed to the resurgence of infec- 
tious diseases in general, but vector-borne diseases in particular. Changes in agricultural 
practices have also influenced the resurgence or the emergence of vector-borne diseases. 
Building new dams which flood areas gives more area for mosquitoes to breed. Some 
irrigation schemes increase mosquito breeding and their population. Clearing forests and 
moving into previously unoccupied areas put people in greater contact with potential 
carriers, bacteria and vectors for infectious diseases. 

1.1.3 Effect of Societal Changes on the Spread of Diseases 

There are lots of societal changes as well as modern living conditions that cause the 
spread of infectious diseases, e.g. AIDS, dengue fever, malaria and so on. 

(i) Plastics: Many of our consumer goods are now packaged in non-biodegradable plastics 
and cellophane. Because of mismanagement, several of these get into the environment, 
collect water, and make ideal larval habitats for all kinds of mosquitoes. 

(ii) Automobiles: There has been an explosion in the automobile industries in the last 
several decades. The tires of these automobiles are discarded after use. They are non- 
biodegradable, and are not only good breeding grounds for mosquitoes, but also they are 
rat harborages. 

(iii) Commerce: Shipping containers provide a vehicle for moving mosquito vectors and 
vermin around the world. 

(iv) Air Travel: It provides the ideal mechanism to constantly move pathogens between 
population centers, especially urban diseases such as dengue, from one part of the coun- 
try or the World to the other. 

(v) Migration: Habitat to habitat migration, due to social, economic, ecological, envi- 
ronmental or political causes, brings infectives in contact with susceptibles in the new 
habitat and thus spreading the disease. 

It is therefore essential to model and analyze the effects of various factors on the spread 
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of infectious diseases and predict the consequences. 

In the following we give a brief literature survey of mathematical modelling of the spread 
of infectious diseases so that research work on modelling and analysis of the problems 
presented in the thesis can be seen in its proper perspective. 


1.2 Brief Literature Survey on Modelling of Infec- 
tious Diseases 


The mathematical theory of infectious diseases, pioneered by Ross, Kennack and M<.K- 
endric, has been an important applied tool, especially for the establishment ot vaccination 
strategies. There have been many notable contributions involving modelling of specific 
infectious diseases such as influenza (Liu and Levin 1989, Castillo-Chavez et al. 1988, 
1989), rubella (Hethcote 1989), Japanese Encephalitis (Tapaswi et al. 1995) and AIDS 
(Anderson and May 1987, Castillo-Chavez et al. 1990, Castillo-Chavez 1989), (set; also 
Hethcote 1976, Hethcote and Yorke 1984, Bailey 1957, 1975, 1979, 1980, 1982, Anderson 
and May, 1979 and May and Anderson 1979). 

In general the spread of such diseases in human populations depends upon various factors 
such as the numbers of infectives, susceptibles, modes of transmission (carriers, v(K;tors 
etc.), social, cultural and economic factors, environmental, ecological and geographical 
conditions (Fuzzi et al. 1996, Dufer 1982). An account of the modelling and study of 
epidemic diseases can be found in lecture notes by Waltman (1974), Hethcote (1974) and 
in the monographs by Bailey (1975, 1982) and the book by Levin et al. (1989). 

As pointed out earlier many infectious diseases are spread by direct contact lx* tween 
susceptibles and infectives, while others get transmitted indirectly, for examphi tuber- 
culosis, typhoid, cholera, malaria, etc. Some diseases spread in the environment and 
are transmitted to the human population by carriers, insects or vectors, which grow in 
the environment due to various household discharges (Cooke 1979, Marcati and Pozio 
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1980, Cairncross and Feachem 1983). The diseases such as tuberculosis, typhoid etc. are 
transmitted both ways i.e. directly as well as indirectly. In direct transmission suscepti- 
bles get infected through meeting or mixing with infectives, while in the case of indirect 
transmission bacteria enter the environment and then contaminate food or water which 
may be later consumed by susceptibles (Hethcote 1976, Gonzalez-Guzman 1989, Shukla 
1986). 

The asymptotic behavior and stability of non-linear epidemic models have been investi- 
gated by many researchers Bailey (1979, 1980, 1982), Cooke (1979), Dietz (1979, 1982), 
Hethcote (1973, 1974, 1976, 1981), Hethcote et al (1973), Hamper (1978) and Wich- 
mann (1979). In particular, Hethcote (1976) presented a qualitative analysis of several 
non-linear models with vital dynamics and with carriers of constant density. Greenhalgh 
(1990) considered SIS models with density dependent death rate and Greenhalgh and 
Das (1992) extended this model for a variable contact rate. In most models quoted 
above the latent period has been assumed negligible but there are diseases in which the 
latent period is important. This aspect can be considered by including an exposed class 
in the usual epidemic models (Hethcote and Driessche 1991, Hethcote 1994). Trawis and 
Lenhart (1987) and Sexena et al. (1993) have analyzed an SIR epidemic model with a 
heterogenous population with the assumption that individuals of different communities 
have different rates of contact. Some models related to nonlinear contact rates (dependent 
upon susceptibles, infectives or population size) have also been proposed and analyzed 
(Capasso and Serio 1978, Hethcote 1994, Liu et al. 1986, 1987). In most of the models of 
infectious diseases it is assumed that the total population size remains constant as birth 
and death rates are equal. But this does not happen if the death rate caused by disease 
is significantly large as in the case of cholera. A similar situation also arises if there is no 
balance between incoming and outgoing population from the region under consideration 
(Hethcote 1994). In view of this, models with demographic structure have been proposed 
and analyzed by considering variation of total population size, which involves birth rate, 
death rate, immigration etc. (Brauer 1995, Hethcote 1994, Mena-Lorca and Hethcote 
1992, Zhou and Hethcote 1994). In particular Hethcote and Van den Driessche (1995) 
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analyzed an SIS model with variable population size and a delay. Gao and Het hcote 
(1992) and Gao et al. (1995) discussed some epidemic models with density d<‘pendent 
birth and death rates and also models with periodicity. Zhou and Hethcote (1904) an- 
alyzed epidemic models with population size dependent incidence for diseases without 
immunity. 

In the papers mentioned above, indirect transmission of diseases is not considtu c’d. Ihu.h- 
cote (1976) analyzed an SIR model with constant density of carriers. Gonzalez-Guznian 
(1989) discussed an SIS model for typhoid fever by considering the effect of flow of bac- 
teria from infected population in the environment through sewage which then affect s tire 
susceptible population by contaminating drinking water. 

Modelling of the transmission of malaria started in the early part of the ninet<ienth 
century (Ross 1911, 1929, Lotka 1923). Ross (1911, 1921) constructed an epidenniological 
model to show that if the mosquito population density is reduced below a thre.sliold level, 
the rate of getting new infections would fall below the rate at which infected per.son.s 
recover leading to elimination of malaria. Since then, various investigators have .studi<‘d 
the spread of malaria using both deterministic and stochastic models (Bailey 1979, 1982, 
Dietz et al. 1974, Macdonald 1953, 1957, Manoharan et al. 1996, Molineaux et al. 
1978 and Radcliffe 1973, 1974). In particular Radcliffe (1973, 1974) has studied the 
utility of catalytic models in the estimation of incidence and prevalence of malaria in 
a hyper-endemic situation. It has been suggested that the models for malaria involve 
two populations of host (human) and vector (mosquito), and are similar in matlumuiticai 
structure to gonorrhea epidemic models (Bailey 1979, 1982, Hethcote and Yorke 
Lajmanovich and Yorke 1976, Nallaswamy and Shukla 1982). 

It may also be noted here that spatial migration of population also plays a very important 
role in the spread of infectious diseases. In countries like India having large population 
density and population distributed over a large number of villages with small inter-village 
travelling distances, deterministic models for infectious diseases involving a spatial aspect 
may be quite suitable. Various attempts have been made to study geographical spread 
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of infectious diseases by considering dispersive migration of susceptibles and infectives 
(Capasso 1978, Capasso and Maddalena 1981, Marcati and Pozio 1980, Webb 1981). 
In particular, Marcati and Pozio (1980) have investigated the global behaviour of a 
vector disease model by considering spatial spread and hereditary effects which may be 
applicable to growth and spread of malaria. A theoretical expression giving the velocity 
of propagation for geographical spread of host-vector and carrier borne epidemics has 
been developed by Radcliffe (1973). The effects of cross dispersal, which arise due to 
spatial influence of one species on another species, have also been considered in the study 
of infectious diseases (Bailey 1980 and Radcliffe 1973). More details regarding spatial 
spread of epidemic diseases, both deterministic and stochastic, can be found in the review 
papers by Bailey (1980) and Mollison (1977, 1995). Effects of dispersal on the linear and 
non-linear stability of the epidemic equilibrium state of the system governing the spread 
of gonorrhea have been investigated by Cooke and Yorke (1973) and Nallaswamy and 
Shukla (1982). 

1.2.1 Environmental and Demographic Effects on Modelling of 
Spread of Infectious Diseases 

Prom the above survey it can be noted that little effort has been made to model the 
effect of environmental changes on the spread of infectious diseases transmitted by car- 
riers, bacteria or vectors, whose population may be dependent on such changes. For 
example, household discharges may be very conducive to the growth of carrier popula- 
tions, such as flies, cockroaches as well as bacterial populations and vector populations 
such as mosquitoes. This increase helps in enhancing the contact rate of infections with 
susceptibles leading to fast spread of the disease (Shukla 1986, Shukla et al. 1987, Misra 
1987). 

It is noted here that the spread of an infectious disease is governed by factors such 
as: (i) the density of the human population, (ii) the density of the carrier, bacteria or 
vector population, (iii) environmental factors such as household discharges, rain, tern- 



12 


General Introduction 


(V) Modelling the Spread of a Carrier-Dependent Infectious Disease in Two Neighboring 

Habitats with Migration in Between, 

(vi) Modelling the Spread of Bacterial Disease in a Population: Effect of Service Providers 
from a Environmentally Degraded Region. 

The density of the carrier, bacteria or vector population has been assumed to follow a 
logistic growth model, the growth rate of which further enhances as the cumulative rate 
of household discharges increases. These problems have been modelled and analyzerl hy 
using the variational matrix method, Liapunov’s second method and computer simulat.ion 
(La Salle and Lefschetz 1961). 


1.3 Summary of the thesis 

The thesis consists of seven chapters and is divided into two parts. The first part deals 
with modelling infectious diseases in a single population whereas the second part deals 
with the spread of infectious diseases into two populations living in two different envi- 
ronmental conditions. 

The first chapter provides the introduction regarding the dynamics of infectious diseases 
including a brief literature survey relevant to modelling and analysis of infectious diseases, 
so that the work done in the thesis can be seen in its proper perspective. 

In the second chapter, we study the effects of the following factors on the spread of a 
carrier-dependent infectious disease using SIS mathematical models: 

(i) A variable carrier population, the density of which follows logistic model and increiises 
further due to household discharges in the habitat, 

(ii) human demography. 

The following two types of demographic factors are considered in the analysis: 

(i) a population with constant immigration and death rate, 

(ii) a population with logistic growth. 

It is shown that the endemic equilibrium is globally stable when the rate of cumulative 
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environmental discharges, conducive to the growth of carrier population, is a constant. 
When the rate is a function of human population density, the endemic equilibrium is 
also locally and globally stable but under certain conditions. The global stability is also 
illustrated by computer simulation. It is shown that if the human population increases 
due to demographic changes, then the spread of the infectious disease increases. It is 
further noted from the analysis that the spread of the infectious disease increases as the 
growth rate of carrier population (caused by conducive environmental discharges due to 
human sources) increases. 

In the third chapter, we have considered a bacterial disease model, which is very much 
similar to the models of Chapter 2 involving carriers. Here the only difference is that 
bacteria are also released by infectives into the habitats continuously and their growth 
rate increases not only because of household discharges but also by continuous release 
from infectives. Here the SIS models for bacterial infectious diseases, (like tuberculosis, 
typhoid, etc.), which are caused by direct contacts of susceptibles with infectives and 
also by bacteria, are proposed and analyzed quantitatively as well as by simulation. It is 
found again that the endemic equilibria are globally stable when the rate of cumulative 
environmental discharges is constant. But when the rate of discharge is a function of 
human population density, the endemic equilibrium is locally and globally stable under 
certain conditions only. It is further shown that under various demographic situations 
mentioned above, the infective population increases with the increase of household dis- 
charges conducive to the growth of bacteria population and the spread of the disease 
increases and so it becomes more endemic. 

In Chapters 4 and 5, malaria models without and with reservoir have been proposed and 
analyzed by considering, 

(i) the effect of cumulative discharges from household sources and (ii) the effects of 
demographics of human and mosquito populations. 

The case of constant as well as variable discharge rates are considered. The threshold 
parameter for spread of the diseases is derived in each case. If this threshold is greater 
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than one, the nontrivial equilibrium in the corresponding case is always feasible and is 
locally asymptotically stable to small perturbations. It is also shown by analysis and 
simulation that this equilibrium is globally stable under certain conditions. I he effect of 
demographic increases in human population is found to increase the spread of malaria. 
It is observed that as the cumulative environmental discharges increase, the mosquito 
population increases, leading to fast spread of Malaria. 

In Chapter 6, a model is considered to study the effect of immigration of population 
from an environmentally degraded habitat to an environmentally cleanei luihit .ii on t he 
spread of an infectious disease. The constant as well as the variable rates of household 
discharges are considered. In each case, the conditions for existence of endemic ec|uilil)- 
rium points have been obtained. It has been shown that under certain conditions, t,h<‘S(^ 
equilibria are locally stable showing that the disease is endemic. Further by simulation, 
it is observed that the endemic equilibrium points are in fact globally stalile iimler local 
stability conditions. It is concluded that as the migration rate increases, ih(> totfii and 
the infective population densities of the region with cleaner environmental condit.ion.s in- 
crease, whereas the total and the infective population densities of the region of degraded 
environmental conditions decrease. 

Lastly in Chapter 7, a mathematical model is proposed by keeping in view I, hat. a poor 
population living in an environmentally degraded habitat spreads the dis<uis(! by inter- 
acting with the rich population by working in their habitat as service providers (without 
migration). The model is presented under similar environmental conditions as in pn!vious 
chapters, here also linear and nonlinear analysis of the equilibria have becm studi(>d for 
various cases. It is noticed that as the interaction between rich and poor cla.sK increases, 
the disease spreads more rapidly in the rich class as expected. It is sugg(!st<id that it is 
in the interest of the rich class to not only clean their own environment but also to chum 
the environment of the habitat, in their neighborhood, where poor people live. 
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Modelling the Spread of Carrier- 
Dependent Infectious Diseases with 


Environmental and Demographic 
Effects 


2.1 Introduction 


Many infectious diseases spread by carriers such as flies, ticks, mites and snails, which 
are present in the environment (Harold 1960, Harry and John 1962, Harry and Kent 
1961). For example, air-borne carriers or bacteria spread diseases such as tuberculo- 
sis and measles; while water-borne carriers or bacteria are responsible for the spread 
of dysentery, gastroenteritis, diarrhea, etc. ( Cairncross and Feachem 1983, Taylor and 
Knowelden, 1964). Various kinds of household and other wastes, discharged into the envi- 
ronment in residential areas of population, provide a very conducive environment for the 
population growth of some of these carriers (Ludwig 1975, Purdpm 1980) . , This enhances 
the chance of carrying more bacteria from infectives to the susceptibles in the population 
leading to fast spread of carrier dependent infectious diseases. Thus unhygienic environ- 
mental conditions in the habitat caused by humans population become responsible for 
the fast spread of an infectious disease. 

In recent decades, there have been several investigations of infectious diseases using deter- 
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ministic mathematical models with or without demographic change (Bailey 1979, 1980, 
Gao and Hethcote 1992, Gonzalez-Guzman 1989, Greenhalgh 1990, 1992, Greenhalgh 
and Das 1995, Hethcote 1976, 1994, Mena-Lorca and Hethcote 1992, Zhou and Het- 
hcote 1994). In particular Greenhalgh (1992) has studied an infectious disease model 
with population dependent death rate using computer simulation. Gao and Hethcote 
(1992) analyzed an infectious disease model with logistic population growth. Zhou and 
Hethcote (1994) have studied a few models for infectious diseases using variou-s kinds of 
demographics. Hethcote (1976) has discussed an epidemic model in which the (.ariier 
population is assumed to be constant. But in general the size of the carrier poiiulation 
varies and depends on the natural conditions of the environment as well as on various 


discharges into it by the human population. 

Thus in this chapter, the effects of the following factors on the spread of an infectious 
disease governed by SIS mathematical models are studied: 

(i) a variable carrier population caused by environmental discharges in the habitat, 

(ii) human demography. 

The following two types of demographic factors are considered in the analysis: 

(i) a population with constant immigration, and 

(ii) population with logistic growth. 


2.2 SIS Model with Constant Immigration 

In this chapter an SIS model with immigration is considered, where the population density 
divided into two classes: susceptibles X(t) and infectives Y(t). It is assumed 
that all susceptibles living in the habitat are affected by a carrier population of (hmsity 
^ C(t), which grows logistically with given intrinsic growth rate and carrying capacity. 
The growth rate of its density is further assumed to increase with the increase in the 
cumulative density of discharges by the human population into the environment. Keeping 
the above in mind and by considering simple^mass action interaction, a mathematical 



2.2 SIS Model with Constant Immigration 


17 


model is proposed as follows, 

X = A-dX - j3 XY -XX C + vY, 

Y = ^ XY + \X C -{vYa + d)Y, 

N = A-dN-aY, 

C = sC(l-j'^-6C + siEC, ( 2 . 1 ) 

E = Q{N)-5oE, 

X + Y = N, 

X(0) = Xo > 0, y(0) = Fo >0, A^(0) = ATo > 0, C(0) = Co > 0, and E{0) =Eo>0. 

Here E{t) is the cumulative density of environmental discharges conducive to the growth 
of carrier population; A is the constant immigration rate of the human population; d 
is the natural death rate constant; (3 and A are the transmission coefficients due to the 
infectives and the carrier population respectively; a is the disease -related death rate 
constant and u is the recovery rate constant i.e the rate at which individual recovers 
and moves to the susceptible class again from the infective class. The constant L is the 
carrying capacity of the carrier population in the natural environment; s is its intrinsic 
growth rate; 5 is the death rate of carriers due to control measures, where s > (5; Si is 
the per capita growth rate coefficient of the carrier population due to the cumulative 
environmental discharges rate Q(N) which is human population density dependent (an 
increasing function of N) and (5o is the depletion rate coefficient of the environmental 
discharges. In writing the model (2.1), we use the term transmission coefficient in the 
sense as used by Anderson and May 1983 (see also de Jong et al., 1995), which means that 
new cases of disease occur at the rates PXY and XXC due to interaction of susceptibles 
with infectives and carriers respectively. 

It can be seen that the region of attraction 

T = |(y; iv, c, E) : 0 < y < iv < :^, o<c<j(^s-5 + si , o<e< , 

is positively invariant and all solutions starting in T stay in T. The continuity of the 
right hand sides of (2.1) and their derivatives imply that a unique solution exists (Hale 
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1969). 

The model (2.1) is analyzed in the following two cases: 

(i) the rate of cumulative environmental discharges Q is a constant, 

and 

(ii) the rate of cumulative environmental discharges Q is a function of the population 
density. The function Q(N) is such that it satisfies following conditions: 


Q(0) = > 0, <3'(iV) > 0, 

i.e Q is an increasing function of N. We consider the form of Q{N) as Q{N) ^ Q„ -I IN, 
where Z > 0 is a constant. 


2.2.1 Case 1: Q = Qa, a Constant 

To analyze the model (2.1) in this case, we consider the following equivalent .system, since 
X + Y = N. 

Y = p{N-Y)Y + X{N-Y)C-{v + a + d)Y, 

N. = A- dN -aY, 

C = sc(l-j'^-5C + siEC, (2.2) 

E = Qa — ^0 E/. 

From the last two equations of (2.2), it is seen that 

lim sup E{t) = ^ and lim sup C{t) = — (s - 5 + Si > 0. 

Thus to see the global behavior of the system it is reasonable to refornmlate syst.em (2.2) 
as follows: 

Y = P{N-Y)Y + X{N-Y)Cm-{i^ + a + d)Y, 

N = A-dN-aY. (2.3) 

Remark: It may be noted that Cm increases as Si, Qa and L incrruise or as S and Sq 
decrease. 

The result of equilibrium analysis is stated in the following theorem. 
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THEOREM 2.1 There exists a nontrivial equilibrium P*{Y, N) regulating the population 
density, where 

_ if + + 4/3 (1 + a) A A-aV 

2^(l + f) ’ <i 

and H =\^-\(l + ^)C„-(v + a + d)^. 

Proof: Setting the right hand side of (2.3) to zero, it is noted that Y is the positive root 
of the following quadratic: 

F{Y) = + y2-|/3^-AC^(|l + ^)-(i^ + a + d)} y-AC^^ = 0 

(2.4) 

and N = As F{0) = - < 0 and F{-^) > 0, Y lies between 0 and 

which implies N is positive. 

Remark: From (2.4), it can be seen that > 0. This implies Y increases as Cm 
increases or as any of Qa, L or Si increases. It may be pointed out that increase in Qa 
gives a more conducive environment to carrier population growth and so spread of the 
epidemic increases. 

2.2. 1.1 Stability Analysis 


The stability result of the equilibrium P* is stated in the following theorem. 


THEOREM 2.2 The equilibrium P* is locally as well as globally asymptotically stable. 


Proof: Let M be the variational matrix corresponding to the system (2.3) at the equi- 
librium point (Y,N), then 


M 


- /Sy-H 


XNCr, 


Y 


—a 


y -f A C, 

-d 


\ 


J 


From the above variational matrix, it is noted that equilibrium point P*{Y, N) is locally 
asymptotically stable. 
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Figure 2 . 1 : Existence of equilibrium point P3 and effect of Qo 011 it. 


2. 2. 2.1 Stability Analysis 


Now we discuss the linear stability of equilibria Pi, P2 and P3 and nonliiKuir stabilit y of 
the only non-trivial equilibrium P3. 

The local stability results of these equilibria are stated in the following tlniorem. 


THEOREM 2.4 The equilibria Pi and P2 are locally unstable and the equilibrium P3 is 
locally asymptotically stable provided 


03 ai 0 
1 02 Oq 

0 03 oi 


> 0 , 


where 03, 02, oi and oq are given in the proof of the theorem. 



2.2 SIS Model with Constant Immigration 


23 


Proof; To prove the theorem, let Mi be the variational matrix corresponding to equilib- 
rium points Pj, for i = 1,2,3, then 


Ml = 


{1/ + a + d) 0 ^ 0 

-a -d 0 0 

0 0 s-^ + U(Co + if) 0 

0 1 0 — < 5^0 


Mo = 


P N - 20 Y -u + a + d p? X{N-Y) 

—a —d 0 

0 0 s — 5 -1- (Qo T IN) 

0 I 0 


/-{PY + ^) pY + XC X(N-Y) 0 \ 

M3= 0. . 

^ 0 0 -j;C Si C 

Vo I 0 -60 J 

From Ml and M2, it is clear that the equilibrium points Pi and P2 are locally unstable. 
To study stability of P3, the characteristic polynomial corresponding to M3 is obtained 


•0^ -f 03 ■0^ 4- 02 -f- aitp -h oo = 0, 


where 


^ » XNC s - . ^ 

O 3 = P Y -\ :: \- — C + 5o + d, 

Y -b 

02 = (^pY+^^^'^(^^C + So + d^ + ^C 5 o + d(^^C + 5 o'^+a{pY + XC), 

Oi = Y H ^ J C Sq + d {— C + 5 o d-{- o:{pY -b A C) C 4- 1 

XNC\ s ^ . , . . . 


ao = ^C^od + oAl siCCiV-y). 

Using the Routh-Hurwitz criteria, the conditions for local stability of this equilibrium 
are as follows: 


03 > 0, 


O3 Oi 

1 02 


03 Oi 0 

>0, 1 02 do 

0 03 oi 


> 0 , 


03 fli 0 0 

1 O2 0(3 0 

0 03 oi 0 

0 1 02 Oo 
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It is noted that first two inequalities are obvious and if the third ineciuaiity is satisfied 
so is the fourth one. Thus the equilibrium P 3 is locally asymptotically sialdo if only the 
third inequality is satisfied. Hence the theorem follows. 

Nonlinear Analysis and Simulation 

We first note that the system (2.6) is bounded by its corresponding system with n = 0 
and Q{N) replaced by Qij). Further, using comparison theorems (LakshmikaiiHiani and 
Leela 1969), we observe that the non-trivial solution of the system (2.fi) is Ijounded by 
its corresponding system with a = 0 and Q{N) replaced by <f£i( ;|)! wliieh can l>e shown 
to be globally stable following similar analysis as in Case I. Hence w(' si)et'ulati- that 
the system (2.6) to be globally stable under the local stability condition in the int <>iior 
of the region of attraction. To support this conjecture and to see the (‘fleet s of vai ious 
environmental and other parameters on the growth of the epidemic (i.e. oti infective 
density) we integrate the system (2.6) by the fourth order Rimge-Kutt.a method using 
the following set of parameters, which satisfy the local stalfility condit ion stated in the 
previous theorem. 

^ = 0.00000051, A = 0.000000021, 1 / = 0.012, a = 0.0005, h 0 , 0 , 4 O.OOl, 

A = 10, (i = 0.0004, s = 0.9, Qo = 20, Si= 0.000002, i =. 0.00005, 0 lOiUiOO. 

We may note here that all the parameters are in per day except tlu! carrying, capacity L 
which is just a number. The equilibrium values of Y, N, C and E have b(‘en found as 
Y = 3407.490, N = 20740.596, C = 38008.224 and E = 21037.027. 

The computer simulation is performed for different initial positions in the following, four 
cases, 

1. r(0) = 100, iV(0) = 12000 , C(0) = 5003 and .S(O) = 300, 

2. y(0) = 1800, iV(0) = 24000, C(0) = 11000 and .B(O) = 9000, 

3. 7(0) = 8330, W(0) = 23000, 67(0) = 4000 and £?(0) = 100, 

4. 7(0) = 4000, Ar(0) = 15500, C(0) = 9500 and S(0) = 455. 

In Fig. 2.2, the infective population is plotted against the susceptible population. 
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Figure 2.2: Variation of infective population with susceptible population. 



Figure 2.3: Variation of infective population with time for different cumulative environ- 
mental discharge rates. 
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Figure 2.4: Variation of infective population with time for different carrying caparituis 
of carrier population. 



Figure 2.5: Variation of infective population density with time for different intrinsic 
growth rates of carrier population. 
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Figure 2.6: Variation of infective population with time for different growth rate coeffi- 
cients of carrier population due to the cumulative environmental discharges. 



Figure 2.7: Variation of infective population with time for dijfferent immigration rates of 
human population. 
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Figure 2.8: Variation of infective population with time for (lifh'renl. 1. 

We see from this figure that for any initial starting condition, the solution curves tend 
to the equilibrium P3. Hence we infer that the system (2.6) is global!} Hiaiih? alHiut 
this endemic equilibrium point Pz{Y, N, C, E) under local staliility condit ions for 
the set of parameters considered, provided that we start away from otiier eijuilibria. In 
Figs. 2. 3-2. 8, the effects of various parameters, i.e. Qo, L, s, sj, and / on the infective 
population have been shown. It is noted from these figures that as thes(‘ param«‘ters 
increase, the population of infectives increases. These results inqily that as t he carrier 
population increases due to household discharges, the spread of the <iis(*asc incn'asi's and 
it becomes more endemic. Also an increase in population density due to immigration 
enhances further the spread of infectious diseases. 


2.3 SIS Model with Logistic Population Growth 

In this section, an SIS model with logistic growth of human populaticm is considered so 
that both the birth as well as the death rates are density dependent in such a manner 
that the birth rate decreases and death rate increases as the population density increases 
towards its carrying capacity (Gao and Hethcote 1992). Keeping in view the consideration 
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of the previous sections, a mathematical model is proposed as follows: 
r riVl r rNl 

X = b-a— N- d+{l-a)~ X - p XY - X XC + uY, 

L A J L A J 

Y = /JXy + AXC- [z/ + a + d+(l-a) ^1 F, 

N = r 1-^ N-aY, (2.10) 

C = sC{1-y)-^C + siEC, 

Jj 

E = Q{N)-doE = Qo + lN-5oE, 

X + Y = N, s>5, 0<a<l, 

A:(0) = Xo > 0, F(0) = Fo > 0, iV(0) = A^o > 0, C(0) = Co > 0 and E{0) = Eq > 0. 

Here b and d are the natural birth and death rates; r = b — d>0 is the growth rate 
constant; K is the carrying capacity of the human population density in the natural 
environment. All other parameters are as defined in the previous section. 

For 0 < a < 1, the birth rate decreases and the death rate increases as N increases to 
its carrying capacity K. When a = 1, the model could be called simply a logistic birth 
model as all of the restricted growth is due to a decreasing birth rate and the death rate 
is constant. Similarly, when a = 0, it could be called a logistic death model as all of the 
restricted growth is due to an increasing death rate and the birth rate is constant. It is 
easy to note that the above model is well-posed in the region of attraction Ti given by, 

Ti = i^{Y,N,C,E):0<Y <N <K, 0<C <^(^s-5 + si. , 0 < E < ■ 


2.3.1 Case I; Q is a Constant Qa 


Similarly to Section 2.2 here also it is sufiicient to consider the following subsystem of 
system (2.10). 



N 


r 


1 


( 2 . 11 ) 
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(2.16) 


+a{pY + XCm-0-- o)-^} = 0 - 

In order that the above quadratic has roots which have negative real parts, it is necessary 


that 


ISY + ^^^ + ^(2N-K)>0, 


(2.17) 


and ^i 2 N-X) + a{fi- ( 1 -a) ■^}Y + 01 \C,„> 0, (2.18) 

The first inequality is obviously true in view of the equilibrium conditions i.e. (2.12) 
and (2.13). Thus only the second inequality gives the condition for linear stability of 
E 2 {Y,N). Hence the theorem. 

Remark: It may be noted that the inequality (2.18) is satisfied for iV > f and 
^ The condition iV > y is also compatible with the condition of ^ > 0. 

Hence from now onward we assume the above mentioned condition. 


Nonlinear Analysis and Simulation 

Using the same Liapunov function as in (2.5) and the system (2.11), we get 

\ ]\T^ r 

V = -(^ + - y? - + N-K){N- Nf 

YY 

\n r 

-{*ia -)}(y -Y)(N- N). 

After choosing ki = + p — {1 — a)^}, we note that V is negative definite in the 

region ^ < N < K such that y < iV < A, where it is assumed that P > (1 — a)y. 
Hence £' 2 ( 1 ^, N) is globally asymptotically stable in a subregion of Tj. 

Also we note that the model (2.11) is globally stable when the disease related death 
rate a is zero and the system (2.11) is bounded by its corresponding system (2.11) with 
a = 0. So using comparison theorems (Lakshmikantham and Leela 1969), it is concluded 
that solution of the system (2.11) is bounded by the solution of the system (2.11) with 
Q! = 0. Therefore, it is speculated that the nontrivial equilibrium point E 2 of model 
(2.11) may be globally stable for iV > y. To illustrate this and to see the effects of 
various parameters on the spread of the disease, the system (2.11) is integrated by the 



Infective population ^ Infective population 


2.3 SIS Model with Logistic Population Growth 


33 



Figure 2.10: Variation of infective population with susceptible population. 



Figure 2.11: Variation of infective population with time for different Cm- 
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fourth order Runge-Kutta method using following set of parameters m the simulation, 
which satisfies the local stability condition (2.17) of equilibrium £’ 2 - 


^ = 0.00000031, A = 0.000000021 ,v = 0.012, a - 0.0005, 

a = 0.3, d = 0.0004, r = 0.0003, K = 50000, Cm - 100000. 

The equilibrium values of Y and N are obtained as T = 7299.645 and N 29087.309. 
Simulation is performed for different initial positions 1, 2, 3, 4 as shown in fig. 2.10. In 
this figure, the infected population is plotted against the susceptible population, from the 
solution curves, we observe that the system is globally stable for this set of parameters, 
provided that we start away from other equilibria. In Fig. 2.11, the infective population 
is plotted against time for different Cm and from this we observe that the infective 
population increases as Cm increases. 


2.3.2 Case II: Q is a Variable 


In this case we consider the following equivalent system of system (2.10) 
(using X + Y = iV), 


Y = ld{N -Y)Y + X{N -Y)C - 
N = N-aY, 

C = sC (1 - - 6 C + Si EC, 

E = Q{N)-6oE = Qo + lN-5oE. 


V + a -f- d -f (1 — a)r 


N 


Y, 


(2.19) 


The result of equilibrium analysis is stated in the following theorem. 


THEOREM 7 There exist the followinQ five eQuilibria, netnely 

(i) Hi(0, 0, 0, ), (a) E2{0, K, 0, ^), (Hi) Ez{Y*, N', 0, E*), which exists if 

pK > {l-a)r + 1 / + a + d, 
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Figure 2.12; Existence of equilibrium point. 


where N* 


-{^(1 - S) - (1 - «)i} + \/Wi-;)-(i-“)tP + 4iF(*' + o + <i) 


oh- 

'^aK 


y* = - ( 1 

a 


N* 

Ik 


N* > 0, 


(iv) EiiO, 0, C, E), where C = j{5 — 5 + 5i^}, 

(v) E^CV, N, C, E), which exists if{i' + a + d) 


S ^ 


E = ^ and 


Proof: The existence of the first four equilibria is obvious. The existence of the fifth 
equilibrium Es is shown as follows. Setting the right hand side of system (2.19) to zero 
and simplifying we get 


r N 

^ ^ (1 - J7)N, 

a K 


- [{p - (1 - a) —}N -{u + a + d+ AC)]!^ - XCN = 0, 


( 2 . 20 ) 

( 2 . 21 ) 


where C = -^{s — d + 

.A.s before, we see that in N-Y plane (2.20) is a parabola and (2.21) is a hyperbola unless 

3 = (1 — a);^ when it is a parabola passing through origin and a branch in the first 

quadrant for A > 0. 

From (2.21), 

'dY' 


the slope 


dN, 


= Y 


kiY + AC + - Y) 


ipY^ + XCN) 


> 0, for y > 0, AT > 0, 


where ki = P — {1 — a)j^, assumed positive. 
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Using these aspects and plotting (2.20) and (2.21) in the first quadrant (Fig. 2.12), we 
see that for the existence of nontrivial Y and N, the slope of (2.20) at (0, 0) must be 

greater than the slope of (2.21) at (0,0), i.e. 


a ^ p-Pa-bd-b^(s-<5-(- f^Qo) 


or (u-bo-Kf) > 

which is same as (2.15) for Qo = Qa- Thus, after knowing Y and N, corresponding values 
of C and E can be calculated as follows; C = j{s — S + SiE} and E Here 

the inequality (2.22) is the ufficient ondition for existence of the fifth equilibrium point 
Es{Y,N,C,E). 

Remark: In both the cases when carrier population is absent, for disease to grow we 
must have a threshold condition as ^ H which is same as mentioned in Gao 

and Hethcote (1992). 


2.3.2. 1 Stability Analysis 

Now we discuss the linear stability of these equilibria and nonlinear stability only of the 
nontrivial equilibrium £ 5 . 

The local stability results of all equilibria are stated in the following thcorein. 


THEOREM 2 . 8 The equilibria Ei, £2 and E^ are unstable. The fourth equilibria E^ is 
stable if 


— f III, y' I 

XC + 1 / + a + d > r and ^XQ y (iz ^ 


otherwise if 


{a ~ r) 


XC < (i/ “h Of + d), 


AC + 1 / + a -h d < r 


or 


r 
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it is unstable and the fifth equilibrium E5 exists. The fifth equilibrium is locally asymptot- 

(23 cti 0 

> 0 and I 1 02 ao 

0 03 oi 


ically stable provided 


03 Oi 

1 02 


are given in the proof of the theorem. 


> 0, where Oq, Oi, 02, and Oz 


Proof: The variational matrices Mi, M 2 , M 3 , M 4 and Ms corresponding to system 
(2.19) at equilibrium points .Ei, £' 2 , £3, £4 and £5 respectively are given by, 


Ml 


(—{y-\-a + d) 0 0 

—a r 0 


0 0 s- <5 + Si 2a 


V 


0 


I 


So 


0 \ 
0 
0 


Mo = 


( fiK — {u + a + d+{l — o)r} 0 
—a r 

0 0 

V 0 I 


XK 

s-s's,^ 
0 


0 \ 
0 
0 

- 5 oJ 


Mz = 


/Ki 

-a r — 

0 

V 0 

N' 


m'n AiV*{l - ^(1 - 0 ^ 

2 rNl 0 o' 

Q(N*) 


K 


0 

I 


s — (5 + Si" 
0 


<5o 


0 

-So/ 


where m'n = -/ 3 ^N *{1 - and m'12 = {/3 - (1 - - ^)N*, 


M4 = 


f — (AC* + + Q; + d) AC 0 0 \ 

—a r 0 _ 0 

0 0 SiC 


V 


0 


I 


0 -SoJ 


and 


Mo = 


/-(/ 3 y + 

—a 

0 

V 0 


XNC 

Y 


) {^-(l-a)j,)Y + \C A(JV-y) 


r — 


K 

2rN 

K 


0 

I 


0 

-§G 

L 

0 


0 \ 
0 

siC 

-So/ 


respectively. Since the matrices Mi, M2 and M3 have positive eigen values so the equi- 
librium points corresponding to these matrices are unstable. 

The characteristic polynomial corresponding to matrix M4 is 
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TT • tVifi "Rniith-Hurwitz criteria, this equilibrium point 
Clearly two roots are negative. Using the Routn tiurwirz 

is locally asymptotically stable if the following conditions are satisfied 


r 


XC + . + a + d > r and + a + (2.23) 


otherwise if 


XC + u + a + d < r or < (p + <t + d). 


it is unstable and the fifth equilibrium exists as mentioned earlier in condition (2.22). 
It is clear from above that the second inequality in (2.23) may not be satisfied even if 
the first is satisfied, but when the first inequality is not satisfied, the second will not be 
satisfied. So violation of any one of the above inequalities gives existence ol the fifth 

equilibrium. 

The characteristic polynomial corresponding to matrix M5 is 

+ 031/1® + a2'4>^ + aiV’ + oo = 0, 


where 


X(N-Y)C . , . jx Y 

a, = ^ + + ^ + 


aU® 4.fl a) 

{N-Y)'N{N-Y) ^ ‘'h<{N-Y) 


+ '-^ + ^0 + 50 > 0, 

K L 


0^2 


4.4) 


4-— C^o + ^ 

Ju 


^{L(^2N - K)} + + S, 

^-(l-a)^|v' + Aa 


+ 


K 


(2N ~ /v) I -“6’ f Aoj 


Ol 


'dY + L(2N -K){\c + S,) + ^(2iV - 


L 


+ {dY + 


Y JK 

+ al{0 - (1 - a)^}Y +\C]Cf + «„), 


a„ = (iiV-t-^)^(2W-A:)^+a^d„l{l3-(l-n)^}f-bAC’! 
+aX{N - y)siCL 

By the Routh-Hurwitz criteria, conditions for local stability of the system are 

az Ol 0 0 


03 > 0, 


1 CI2 


> 0, 


03 Ol 0 
1 O2 Oq 

0 03 Ol 


> 0 and 


1 0,2 oo 0 

0 03 Ol 0 

0 1 02 Oq 


> 0. 
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Clearly the first inequality is obvious. If the second and the third inequalities are satisfied, 
so is the fourth one. Hence the equilibrium point E5 is locally asymptotically stable if 
the second and the third inequalities are satisfied. 

Remark; It is seen that second is satisfied for ^ > Y and ^ So in this case 

only third inequality is the condition for local stability of E5. 

Nonlinear Analysis and Simulation 

As before we speculate that the nontrivial equilibrium point E5 of the model (2.19) may 
be globally stable under the local stability conditions. To illustrate this and to see the 
effects of various parameters on the spread of the disease, the system (2.19) is integrated 
using the fourth order Runge-Kutta Method by taking Q{N) = Qo + IN and using the 
following set of parameters in the simulation, which satisfies the local stability condition 
mentioned above. 

^ = 0.00000031, A = 0.000000021, u = 0.012, a = 0.0005, 

5 = 0.6, 5o = 0.001, r = 0.0003, d = 0.0004, a = 0.3, K = 50000, 

s = 0.9, Qo = 20, Si = 0.000002, I = 0.00005 and L = 100000. 

All the parameters are in units of per day except the carrying capacity L, which has the 
same dimension as C. 

The equilibrium values of Y, N, C and E have been found as 

F = 6121.797, A = 35716.819, (7 = 38174.626, .^ = 21785.839. 

In Fig. 2.13, the infected population is plotted against the susceptible population and 
from the solution curves, it is concluded that the system appears to be globally stable 
for this set of parameters. In Figs. 2.14-2.19 the effects of various parameters, i.e Qo, L, 
s, Si,r and I on the infective population have been shown. It is noted from these figures 
that as these parameter values increase, the infective population increases and we have 
similar conclusions regarding the spread of the infectious disease as discussed earlier. 
Remsurk; When AT — > 00 , the model (2.10) coincides with the model with exponential 
growth of the population, though here it is not considered explicitly. 
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Figure 2.17: Variation of infective population with time for different growth late coeffi- 
cients of carrier population due to the cumulative environmental dischargcis. 



Figure 2.18: Variation of infective population with time for different intrinsic growth 
rates of human population. 
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Figure 2.19: Variation of infective population with time for different 1 . 

2.4 Conclusions 

In this chapter SIS models for carrier dependent infectious diseases, like cholera, diar- 
rhea, etc. caused by direct contact of susceptibles with infectives as well as by carriers are 
proposed and analyzed. The following two types of demographics are considered, (i) con- 
stant immigration, and (ii) logistic population growth. For both the models, equilibrium 
analysis is presented. It is shown that in the case of constant immigration the endemic 
equilibria is globally stable when the rate of cumulative density of household discharges 
is a constant. When it is a function of total human population, the endemic equilibria is 
again shown to be locally and globally stable under certain conditions. The later result 
is shown by computer simulation. But in the case of logistic population growth, the local 
stability of the nontrivial equilibria in both the cases is guaranteed only under certain 
conditions. By computer simulation it is shown that under the local stability conditions, 
the nontrivial equilibrium appears to be globally stable in both the cases. It is concluded 
from the analysis that if the growth of carrier population caused by conducive household 
discharges increases, the spread of the infectious disease increases. Also when the human 
population increases due to demographic changes, the infectious disease spreads even 
further and becomes more endemic. 




Chapter 3 


Modelling Bacterial Disease with 
Environmental and Demographic 
Effects 


3.1 Introduction 


The widespread occurrence and increasing incidence of various infc'ctious diseases such 
as gonorrhea, tuberculosis, hepatitis, measles, influenza etc. are bec^omiiig major public 
health problems in most tropical countries and elsewhere (Bailey 1957, 1979). In general, 
spread of infectious diseases in human population depends upon various fact ors such as 
the numbers of infectives and susceptibles; modes of transmission (carritu's, vectors, etc.); 
socio-economic factors; environmental, ecological and geographical and similar factors 
(Fuzzi et al. 1997, Dufour 1982). A detailed account of modelling and study of epidemic 
diseases can be found in literature in the form of lecture notes, monograplis and similar 
survey. (Waltman 1974, Bailey 1975, 1982, Hethcote 1976, Hethcote et al. 1981). The 
population biology of infectious diseases has also been presented by Anderson and May 
(1979). While several infectious diseases are spread by direct contact between susceptibles 
and infectives, there are some diseases (such as tuberculosis and typhoid) which are also 
transmitted to the human population indirectly by the flow of bacteria from infectives 
into the environment, for example by contaminated water used by susceptibles. Gonzalez- 
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Guzman (1989) analyzed an SIS model for the spread of typhoid by including the effect 
of indirect transmission with the flow of bacteria from infective in the environment. It is 
well known that bacteria grow and survive in large numbers in almost every environment 
on earth including polluted water in ponds, lakes, rivers, in acidic and alkaline waters, 
in the presence of toxic substances, and also in air and soil (Dufour 1982). Also when 
various kinds of household wastes are discharged into the environment, some disease 
causing bacteria find these discharges very conducive to their population growth. For 
example, cholera bacteria can survive by sheltering beneath the mucus outer coat of 
various algae and zooplankton. In the case of typhoid, the typhie bacteria multiply in milk 
products and other food wastes. Warmer water may increase algae blooms, helping vibrio 
cholera to multiply and perhaps even promote the emergence of new genetic strains. Thus 
increase in population of bacteria in the environment enhances the spread of infectious 
disease in human habitat. 

Further the growth of human population in a habitats plays an important role in the 
spread of infectious diseases (Hethcote and Zhao, 1994, May and Anderson 1979). In 
general human populations vary and population change takes place in a hal)itat due to 
immigration, growth rate, etc. 

In this chapter, therefore SIS models for the spread of infectious diseases are proposed 
and analyzed by considering environmental and demographic factors. 


3.2 SIS Model with Immigration 

We consider here, an SIS model with immigration. The disease is assumed to be sjjread 
directly by infectives as well as by the flow of bacteria into the environment from in- 
fectives. The total population density N(t) is divided into a susceptible class X(t) and 
an infective class Y(t). It is assumed that all susceptibles living in the habitat are af- 
fected by the bacteria population, whose density grows logistically with a given intrinsic 
growth rate and carrying capacity. The growth of the bacteria density is further assumed 
to increase as the cumulative density of discharges into the environment by the human 
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population increases. Keeping this in mind, a mathematical model is proposed as follows: 


X = A-dX -{PY + XB)X + uY, 
Y = {^Y + XB)X -{u + a + d)Y, 
N = A-dN-aY, 

B = sb(i- 
E = Q{N)-5oE, 

X + Y = N, 


B 

I 


+ S1Y-S0B + 6 BE, 


X(0) = Xo > 0, y(0) = Fo > 0 , x(0) = Xo > 0, B(0) = Bo > 0 and E{0) = > 0. 


Here E{t) is the density of cumulative environmental discharges conducive to the growth 
of bacteria population; A is the immigration rate of human populataon assumed to be a 
constant; d is the natural death rate constant; /3 and A are the transmission coefficients 
due to the infectives and bacteria population respectively; a is the dis(!ase rcdated death 
rate constant; v is the recovery rate constant; s is the intrinsic growth ra.te of t.he bacteria 
population; L is the carrying capacity of the bacteria population; Sq is the <leath rate of 
bacteria due to control meaisures; Si is the rate of release of bacteria from t.he infective 
population in the environment; 5 is the rate of growth of the bacteria population due 
to environmental discharges; Q{N) (assumed to be increasing in N) is the cumulative 
rate of environmental discharges which may be population density de[>endent and d'o 
is its depletion rate coefficient. We note that s > Sq. In writing the model (3.1), we 
assume that new cases of disease occur at rates PXY and XX B due to the interaction of 
susceptibles with infectives and bacteria respectively (Anderson and May 1<183). 

The region of attraction corresponding to (3.1) is given by, 


T = UY,N,B,E):0<Y<N<^,0<B<B, 


max 5 


0<E< 


Qij) 


where 


2s 


s — So + 



+ 


\ 


5 — 5o + (J' 


9Al 

s. 


2 


d.S'.S'i A 


1 


is positively invariant and all solutions starting in this region T stay in it. The continuity 
of the right hand sides of equations 3.1 and their derivatives implies that a unique solution 
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exists (Hale 1969). We analyze the model (3.1) in the following two cases: 

(i) the cumulative rate of environmental discharges Q is a constant, and 

(ii) the cumulative rate of environmental discharges Q is a function of total population 
density, which we consider as Q = Qo + IN, where Qq and I are constants. 

3.2.1 Case I: Q is a Constant Qa 

Since X-\-Y = N, the system (3.1) can be studied by the following autonomous subsystem 

Y = + \B){N -Y)-{v + a + d)Y, 

N = A-dN-aY, (3.2) 

B = sb(i-j^+ siY - sqB + 6BE, ' 

E = Qa-SoE. 

It may be noted that ^ as t — >■ oo. Therefore it is reasonable to reformulate the 
above system by using the asymptotic value of E as follows: 

Y = {pY + \B){N -Y)-{u + a + d)Y, 

N = A-dN-aY, (3.3) 

B = sB(i-j'^+siY-soB + 5B^. 

The result of an equilibrium analysis is stated in the following theorem. 

THEOREM 3.1 For the system (3.3), there exist two equilibria, namely (i) Ex ^0, 
and (ii) E 2 {Y,N,B). 

Proof: Existence of the first equilibrium is obvious. Existence of the second is shown eis 
follows. Setting the right hand side of (3.3) to zero, we get the following two equations, 

andy = — + (3.5) 

Si I L I OqJJ 



48 


Modelling Bacterial Diseage 


with. Environinental and Demographic Effects 



From (3.4), we have 
(i) when 5 = 0, 


y = 0 or y = 


iA-jir + a + d) ^ j 

/ 5{1 + |) 


We see that n <0. when ^ < {u + <, + d)^ndY, > Owl.™ > (.'t»+4 

A 

(ii) B oo when Y 


(iii) the slope at (0, 0) 



-AA 

d 


iA^[u + a + d) 


< 0, if ^ + a i d) 

a 

> 0, if ~ < 4 1 d), 

a 


also the slope at (0, y) is 


/ dy \ _ A(i/ + a 4- d) 
[dBj “ ^{^-{u + a + d)} 


< 0, if ^ 1 n i d) 

d 

> 0 ,if ^ < (i/ i fv 4 d). 
d 


From (3.5), we have 

(i) when y = 0, 5 = 0 or B = ^(s - 5o 4- ^^) = Bi (say). 
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(<*5) ( 0 , 0 ) 

(PY 2s 

(iii) — — p > 0, i.e slope increases with the increase of B. 

Now plotting (3.4) and (3.5) in Fig. 3.1, we get an intersecting point {Y,B), where 
Y < irrespective of ^ greater or less than { 1 / + a + d). Then N can be derived 
from N = , which is positive as F < 

3. 2. 1.1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results are 
stated in the following theorem. 

THEOREM 3.2 The equilibrium Ei is unstable and the equilibrium E 2 is locally asymp- 
totically stable. 


-So + 5% 


So 


Si 


< 0 and 


'dY' 

MB 


(Bi,0) 


So + 

— ^ > 0 , 

Si 


Proof; The variational matrix M corresponding to system of equation (3.3) at (y, N, B) 
is given by 

-Y)-{PY + XB + u + a + d) /3Y + XB X(N-Y) \ 

M = I —a —d 

Si 0 


X{N - Y) 
0 


s - So + - 2f B/ 


Now, the variational matrix Mi at equilibrium point Ei ^0, ^,0^ is given by 


Ml = 


/ J3j — ( 1 / + a + d) 0 


Ad 


\ 


V 


—a —d 0 

Si 0 s — So + j 

We see that one eigenvalue corresponding to the above matrix is — d and others are given 


by the following quadratic, 

^s — So + ^ d^ Ip + ^s — So + — t' + Q! + d^ 

Clearly the equilibrium £^i is locally asymptotically stable if 
XQa . r . cQaA 


AAsi 


= 0 . 


s-so 


+ d^ - (i^ + o: + d) < 0, ^s - So + - 1 / + a + d| - — ^ > 0. 
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TaUng the maximnm of the right hand side of the above inequality, we get 

-j- X]j 

u ^ LA . 

^ 2d 

Thus for above choice of h and ^ 2 , is negative definite in the interior of the region of 
attraction. Hence the equilibrium El of the system (3.6) is globally stable provided we 
start away from the disease free equlibrium. 

Hence the model (3.3) is bounded by the system (3.6), which is globally stable. So using 
comparison theorems (Lakshmikantham and Leela 1969), the solution of (3.3) is bounded 
by the solution of (3.6). Hence we speculate that the nontrivial equilibrium point of (3.3) 
may be globally stable. To illustrate this, the system (3.3) is integiatcd using the. fouith 
order Runge-Kutta Method and using the following set of parameters in the simulation 
(Greenhalgh 1990, 1992). 

/3 = 0.00000031, A = 0.00000000021, v = 0.012, a = 0.0005, so = 0.G5, 8 = 0.000002, 

do = 0.001, A = 10, d = 0.0004, s = 1, Qo = 20, Si = 10, L = 5000000. 

The equilibrium values of T, N and B are found as, 


Y = 1466.742, N = 23166.57 and B = 1986910.242. 


Simulation is performed for different initial starting points. 

In 1, Yo = 3610, No = 5092, Bq = 2180 . 

In 2, To = 210, No = 24692, Bo = 5180 . 

In 3, Yo = 4222, No = 24500, Bq = 1000. 

In 4, To = 5000, No = 15000, Bo = 90. 

In Fig. 3.2, we have plotted the infective population against tlu* susceptible population. 
Prom the solution curves, we conclude that the system is globally stable about the en- 
demic equilibrium point (7, N, B) provided we start away from di.sease free equilibrium. 
Also in Figs. 3. 3-3.8, the variation of the infective population is shown for different values 
of s, Si, 5, L and A. It is concluded that with the increase of these parameters, the in- 
fective population increases. This implies that with increase in the cumulative discharges 
into the environment, the spread of bacterial dependent infectious diseases increase. 
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As before, plotting (3.7) and (3.8), we get two intersecting points (0,0) and (Y,B), cor- 
respondingly we get two equilibrium points (o, o) and E^iY, N, B). 

By using the variational matrix method, it is easy to show that the equilibrium point E^ 
is unstable and is locally asymptotically stable. 

To study the global behavior of the system about the nontrivial equilibrium point 
we consider the following Liapunov function: 


y = yf + - Nf + - B)^ 

which gives 

V = (Y-Y){l3Y{N-N) + j3N{Y -Y)-I3{Y + Y){Y-Y) + XB{N-N) 

+ XN{B -B)- XB{Y -Y)- XY{B - 5 )} - (u + d){Y - Yf - kid{N - N)^ 

-h {B - Bf + kMY -Y)(B- B) 


= -\^Y + XB + 


XB{N -Y) 
Y 


(F - Yf - krd[N - Nf + (/3y + XB){Y -Y){N - N) 

sx^ 


+ {X{N - Y) + k2S^}{Y -Y){B-B)- k^ {B - B)\ 

Using Sylvester’s criteria, V will be negative definite provided the following inequalities 
are satisfied: 


(i) 


Ako 


pY + XB XB{N-Y) 
2 F 


s r, SiF 

yB + -B- 
B B . 


> 


{X{N-Y) + k2Siy 


or 


2{PY+XB)k2 > {xiN - F) - A:2Si}' (3.9) 


Now choosing ^2 = above inequality is satisfied provided F > 0, (jB = 0 or 

B^O). 



which implies 


y/3Y + XB) 
2 


kid > {/3Y + XBf 


2kid > {0Y + XB)foTY>0kB> 0. (3.10) 

SFo A .1.3.9-6.9.0.— 
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Figure 3.4: Variation of infective population with time for different growth rate of bacteria 
due to infective human population. 



Figure 3.5. Variation of infective population with time for different growth rate of bacteria 
corresponding to environmental discharges. 




3.2 SIS Model with Immigration 


55 


L = 6000000 





I 1000 - / 

^ ; / 

- 800 if 

600 f 

400 ‘ ' ' ' 1 1 1 1 1 1 

0 5 10 15 20 25 30 35 40 45 50 

Time in years ^ 

Figure 3.6: Variation of infective population with time for different carrying capacity of 
bacteria population. 
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Figure 3.7: Variation of infective population with time for different rate of immigration 
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Figure 3.8: Variation of infective population with time for different rate of cumulative 
environmental discharges 

3.2.2 Case II: Q is a Variable 


In this case, let us consider the original model (3.1) with Q{N) = Qo + IN. 

Y = {I3Y + XB){N ~Y)-{u + a + d)Y, 

N = A- dN-aY, (3.11) 

B = s(^B- 
B = Qo + IN — SqB. 

The result of an equilibrium analysis is stated in the following theorem. 


B^ 


+ siV — SqB + 6 BE, 


THEOREM 3.3 There exist two equilibria, namely (i)E\ fo, 4,0, ^ and 

(ii)E 2 {Y, N, B, E). The second equilibrium exists if 


Lf rQo . SIA., siSod 

-{s-So + 5— + tv) < ^7^- 

-S do Sod Slot 


Proof: We prove the existence of the second equilibrium E 2 as follows. Setting the right 
hand side of the above system to zero, we get 


P_ Qo + IN 
d So 


d 


(3.12) 
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d 


{iy + a + d)-X(l + 3)sl y - = 0, 

a J a 




Here (3.13) is the same as (3.4). Prom (3.14), we get, 

(i) when y = 0, S = 0 or 5 = - js - So + <5% + f4l = > 0 (say), 

■S ( oq ood J 


/'dy\ _ _ , 


[dBj^0,0) 

i 


< 0 


and 


) {Bifi) 


+ + f^ ' 


> 0 if Si < 


Si^od 

5l(x ’ 


(3.13) 


(3.14) 


(3.15) 


(iii) (S) ^ B = 

As before plotting (3.14) and (3.15), we will get a positive intersecting point (S, y) under 
condition (3.15). Then E and N can be calculated from (3.12) as Y < 


3.2.2. 1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results of 
these equilibria are stated in the following theorem. 

THEOREM 3-4 The first equilibrium Ei is always unstable, and the second equilibrium 
E 2 is stable provided 03(0501 — OoOs) — a\ > 0, where Oq, Oi, 02 and 03 are given 

explicitly in the proof of the theorem. 

Proof: To study the stability analysis let us consider the variational matrix M corre- 
sponding to the system (3.11), 


/fiN- 

- 2fiY -XB-iu + a + d) 

fiY + XB 

x{N - y) 

0 \ 


—a 

-d 

0 

0 


Si 

0 i 

i - So + 5E-^ 

6B 


0 

1 

0 

-5oJ 
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The variational matrix Mi corresponding to the system of equation (3.11) at equilibrium 
point El (^0, f , 0, is given by 


Ml = 


V 


{y + a-\-d) 

0 

\A 

d 

0 \ 

— OL 

-d 

0 

0 

Si 

0 

j;-Q0 + X 

0 

0 

1 

0 

-So) 


Here two characteristic roots are —d and —So and the other roots are given by the 
following quadratic equation, 

'll? - |^-(i/ + a + d) + s-so + 5^^-^-^|v' 

A . I r -Qn + ^ 1 . A 


+ 


{^--(z/ + Q; + d)}|s-So + (^-^^ I Asi^ -0. 

By the Routh-Hurwitz criteria this equilibrium is unstable because in the above quadratic, 
the coefficient of tp and the constant term are not -positive simultaneously. 

The variational matrix M 2 corresponding to equilibrium point E 2 {Y,N,B,E) is given 

by 

+ pY + XB X{N-Y) 

-d 0 . 


M2 = 


—a 




Si 

0 


I 


0 


-SoJ 


The characteristic polynomial is given by 

+ 03 Ip^ + CL2lp^ + Cll'0 + Go — 0; 


where 

O.Z = 


, XNB\ ^ s - sif , 

i PY H — 7;— j + d + — jB ^ — — + So, 


02 = [PY + 


Y 

XNB' 

Y 


L B 

^ Sa siY 
d+-B+ — 
^ L B 




+ <^0 I + 


+a(^y + AB) - A(.V - 


Oi 


BY + 


XNB\ 


- 


d(iB-^~ + So] + (pY + ^ 4 ^ + d 


B 


Y 


yB + 

,L B 


75 + ^) 
.B B y 


SiT' 


^0 


+a{PY + XB) fp + fi + ^ 0 ^ _ X{N - Y)si{d + So), 


ao = {PY + 


XNB\ 


Y 




-B + 


fil' 

B , 


So + a{pY + XB)So(^B + ^ 
\L B 


-X{N - Y){dsiSQ - alSB) > 0. 
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By Murata (1977), conditions for local stability of the system are 


as > 0, 


as ai 
1 Cl2 


> 0 , 


as ai 0 
1 0,2 
0 as ai 


> 0 , 


5m 

are 



03 

ai 

0 

0 

1 

0^2 

Oq 

0 

0 

<^3 

Oi 

0 

0 

1 

O2 

Oq 


> 0. 


The first two inequalities are obvious. If the third is satisfied then so is the fourth one 
as ao > 0. Thus the equilibrium B 2 is locally asymptotically stable under the condition 
mentioned in the theorem. 


Nonlinear Analysis and Simulation: 

The global stability of B 2 can be speculated as before. To show this, the system (3.11) 
is integrated using the fourth order Runge-Kutta Method and using same parameters 
as mentioned in the previous subsection, with Qo = Qa and an additional parameter 
I = 0.000005, which satisfy the local stability conditions. The equilibrium values of 
Y, N, B and E have been found as 

T = 1467.388, iV = 23165.766, 5 = 1988063.370, 5 = 20115.828. 

In this case also simulation is performed for different initial positions 1, 2, 3 and 4, shown 
in Fig. 3.9, where the infective population is plotted against the susceptible population. 
From the solution curves, it appears plausible that system is globally stable about the 
endemic equilibrium point (y, iV, 5, E) under the local stability condition provided 
we start away from the other equilibrium point. Also in Figs. 3.10-3.15 the variation 
in the infective population is shown for different values of s, Si, 5, L, I and A. It is 
concluded that with the increase of these parameters , the infective population increases 
showing that the spread of bacterial infectious disease increases with bad environmental 
conditions as well as with immigration. 
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Figure 3.13: Variation of infective population with time for different carrying capacity of 
bacteria population. 



Figure 3.14: Variation of infective population with time for different 1 . 
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Figure 3.15: Variation of infective population with time for different rate of immigration 
of human population. 

3.3 SIS Model with Logistic Growth 

In this case, we consider an SIS model in which logistic growth of human population 
is assumed. Here both the birth and death rates are density dependent. The birth 
rate decreases and the death rate increases as the population size increases towards its 
carrying capacity (Gao and Hethcote 1992). The mathematical model is given by the 
following set of equations, 

V = j^b-ar^'^N-!^d+{l-a)'^^X-{pY + XB)X + iyY, 0 < o < 1, 

Y = (0Y + XB)X -l^u + a + d+{l-a)'^^Y, 

N = r(l-^^N-aY, (3.16) 

B = sB (l-j^+ siY - SoB + 5BE, 

E = Q{N)-5qE = Qc, + IN -5^E. 

Here b and d are the natural birth and death rates; r = 6 — d>0is the growth 
rate constant; K is the carrying capacity of the environment corresponding to the human 
population. All other parameters are as defined in the previous subsection. For 0 < a < 1, 
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the birth rate decreases and the death rate increases as N increases to its carrying 
capacity K. When a = 1, the model could be called simply a logistic birth model as all 
of the restricted growth is due to a decreasing birth rate and the death rate is constant. 
Similarly, when a = 0, it could be called a logistic death model as all of the restricted 
growth is due to an increasing death rate and the birth rate is constant. We see that the 
region of attraction T' corresponding to (3.16) is 


T' = J (y, W, E) : 0 < r < ^, 0<N <K, 0<B< 0<E< 


Q{K) 


and the model is well-posed in the region T', where B^ax is given by 


Bmzx — 




issiK 

L 


3.3.1 Case I: Q is a Constant Qc 


Since X + Y = A^, it is sufficient to consider the following equivalent system of (3.16), 

y = (/5y-hAB)(iV-y)-|z/ + a; + d+(l-a)^|y, 

/ N\ 

N = ^ kJ ^ ~ 

B = sB (l-j^+siY-soB + 5BE, 

E = Qa-6oE. 

Putting the asymptotic value of E i.e. Era ~ ^ iii the above system of equations we get 
the following subsystem; 

y = iPY + XB){N-Y)-{u + a + d+(l-a)~}Y, 

N = r{l-~)N-aY, ( 3 . 18 ) 

B = sB (1 — — ) + s^y — sqB 5B^^. 

^ On 


The result of equilibrium analysis is stated in the following theorem. 
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Figure 3.16: Existence of equilibrium point when a > r. 
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THEOREM 3.5 There exist the following four equilibria, namely 
(i)PiiOM, (n)P2{0,0,B*), whereB* = ~{s-so + 5^}>0, (Hi) Ps{0,K,0) 
and (iv) PifY, N, B), for a> r, this exists if 

u + a + d > 


Proof: The existence of the first three equilibria is obvious. We will show the existence 
of the fourth equilibrium point P4 by the isocline method. Setting the right hand sides 
of (3.18) to zero, we get the following for N 0 and N ^ K. 



T 

Y = -N 
a 



(3.19) 


jB" - (s - 'Jo + 5^) 

a 

11 

(3.20) 

htf- 

K 

[^-N{l-^) + XB + v + Oi + d\ -(1 
V a K J Of 

-^) + AB=0, 

(3.21) 


where ki = jdK — (1 — a)r. 
Also (3.21) can be written as 


B = 




h- 


^r K'' 


N 


a 


aK 

I3r 


{u + a + d) 


Xr 




or 


/?r {N - K){N - Ni){N - N 2 ) 
XaK {N - N*) 


(3.22) 


where Ni is the negative and N 2 is the positive root corresponding to the quadratic in 


the numerator of the above equation and N* = {1 — f)K. 


Also when N = 0, B = -■^-777—^^ = Bi (say). 

A(l-p) 


Clearly in the N-B plane, (3.20) is an ellipse passing through (0, 0) and (K, 0) with major 

/K L ( rQaW 

and minor axes parallel to the coordinate axes and origin at ^ 

We consider the following two cases: 

Case 1: a > r, i.e. N* is negative 
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In this case, (a) when ki>u + cx + d>0^0 < N 2 < K and 
(b) when kx is negative oi^<kx<i' + (x-\-d=^K < N 2 < 00 . 

The graphs of (3.20) and (3.22) are shown in Fig. 3.16. We note that for the existence of 
positive intersecting point with N less than carrying capacity K, we must have Bi > B 2 , 
where B 2 is given by (3.20) when W = 0. We see that B 2 = B*. Thus in this case, the 

condition for existence of (N, B) is 

, , j ^ \ R* 

i/ + a + d > A±> . 

r 

Case II: 0 < O' < r, i.e. N* is positive. 

In this case, (a) when fci>z/ + a; + d>0=4>0 < N 2 < K and 
(b) when kx is negative ovO<kx<i' + cx + d=^K < N 2 < 00 . 

The graphs of (3.20) and (3.22) are shown in Fig. 3.17. It is noted that the equilibrium 
value Ni > Y (^) > f and (ii) N 2 > j, which give rise to the following 

sufficient conditions on the parameters, 

0<«<J, + a + (3.23) 

3. 3. 1.1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results of 
these equilibria are stated in the following theorem. 

I 

THEOREM 3 . 6 The equilibrium Px is unstable, the equilibrium P 2 is locally asymptoti- 
cally stable provided 

— — —\B* > v + a + d, 
r 

otherwise if 

(cx, — t] 

-XB* < u-i-a + d 

r 

it is unstable. The equilihviuTn is unstable and the equilibrium point P 4 . is locally 
asymptotically stable if the condition aa > 0 and 0 x 02 -^ a^ > 0 are satisfied, where ax, a 2 
and as are given explicitly in the proof of the theorem. 
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Proof: The variational matrix M at (Y, N, B) corresponding to the system (3.18), is 
given by 

'p{N-2Y)-XB-{u + a + d+{l-a)^} pY + \B X{N-Y) 

-a (r - 2rf ) 0 

^1 Os — So + 

The variational matrix Mi at equilibrium point Pi(0, 0, 0) is 


M = 


Ml 


/'—(ly + a + d) 0 
— CK r 

Si 


0 

0 


0 s-So + Sf 


Clearly one eigenvalue of the above matrix is positive, this equilibrium is unstable. 
The variational matrix M2 at equilibrium point P2 is 


M2 = 


/ — {AP* + 1/ + Q; + d} XB 


V 


— O' 
Si 


0 

r 0 

0 -(s-so + d^), 


Clearly one eigenvalue is negative and the other eigenvalues are given by the following 
quadratic 

+ {XB* + i/ + Q; + d-r}^ + aXB* - r{XB* + iy + a + d)=0. 

{oi — r) 


So this equilibrium is stable provided 


-XB* > 1/ + o + d. Clearly for a < r, this 


equilibrium is unstable. 

The variational matrix M3 at the equilibrium point P3 is 
/ /8K — {1/ -h O' -h d + (1 — a)r} 0 


M, 


V 


—a 

Si 


XK \ 
—r 0 

0 s - So + ) 


The characteristic polynomial corresponding to the above matrix is 

(r -j-'ip) 1/)^ — |s - So + + PK -(t' + Q; + d+(l- a)r)| ip 

4" ^s — So + d — ( 1 / + q; + d + (1 — 0')^) — XKsx = 0 

So one root is negative and the other roots have negative real part if 

ySJT — (1 — o)r — (i^ + a + d) + 5 — So + d^ < 0, 

^0 


(3.24) 


(3.25) 





Figure 3.19: Variation of infective population with time for different intrinsic growth rate 
of bacteria population. 



Figure 3.20; Variation of infective population with time for different growth rate of 
bacteria population due to infective human population. 
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Figure 3.21: Variation of infective population with time for different growth rate of 
bacteria corresponding to environmental discharges. 



Figure 3.22: Variation of infective population with time for different carrying capacity of 
bacteria population. 
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Figure 3.23: Variation of infective population with time for different rate of immigration 
of human population. 

3.3.2 Case II: Q is a Variable 

In this case let us consider following system of equations 

rN 

Y = {pY + \B){N-Y)-{u + a + d+{l-a)—}Y, 

N = r(l-^'^N-aY, (3.27) 

B = sb(i-j'^-\- siV - SqB + SBE, 

E — Qq + IN — 5qE. 

The result of an equilibrium analysis is stated in the following theorem. 

THEOREM 3. 7 There exist the following four equilibria^ namely 
(i) Ei(0, 0, 0, f ), (ii) E2i0, 0, B\ ^), where B* = ^ {s - So + 6^} , 

(in) £'3(0, K, 0, and (iv) E^(Y,N,B,E), this exists if 

Assir ^ 

^ w ■ 

For a > r, we should have an additional condition as i/ + a + d > 





a > r. 
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Proof; The existence of each of the first three equilibria is obvious. We present in the 
following, the proof of the existence of the fourth equilibrium £' 4 . Setting the right hand 
sides of (3.27) to zero, we get the following for AT 7^ 0 and N 


T N 


s 

L 




B = 


{N K) J 

q;2£:2^ M 

f . r2 CX.K i 
1 ' 

i h ySr'l 

\N 


Ar j 
OiK 1 

[n 

- (1 - ^)k] 

\ 


(3.28) 

(3.29) 

(3.30) 


where ki = — (1 — a)r. 

Clearly (3.29) is an ellipse if 


4ssir 

aKL 


> 


'5V 


(3.31) 


which passes through (0,0), (K,0), (0,J5*) and (^K, j{s - sq + £(Qo + IK)}^ Also 
(3.30) is the same as (3.22). Hence as before, plotting (3.29) and (3.30) in the N-B plane, 
we will get a positive intersecting point (N, B) (Figs. 3.24 and 3.25) and corresponding 
to it we get Y and thus we have the fourth equilibrium point £ 4 . 


3.3.2. 1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results of 
these equilibria are stated in the following theorem. 

THEOREM 3.8 The equilibrium E\ is unstable, the equilibrium £2 is stable if 

\B’ 4-1/4- or + d > r and — —XB* > u + a + d, 

r 

otherwise if 

A£* + 1 / a + d < r or — —XB* < u + a + d, 

r 

it is unstable. The equilibrium £3 is unstable and £4 is locally asymptotically stable 
provided Oq > 0, aza^ — Oi > 0 and ai{aza 2 — oi) — Uq 03 > 0, where oq, Oi, 02, and az 
are given explicitly in the proof of this theorem. 



78 


BsLCtGrisJ Dis69iS6 with Environniciitsl 3*n(l D6nio§r3»phic Effects 


Proof: The variational matrix M at (Y, N, B, E) is 


mil 

mi2 

X{N - Y) 

0 

\ 

—a. 

r-^^N 

0 

0 


Si 

0 

s-so+5E-^B 

5B 


0 

1 

0 

-So 

y 


where mu = — 2^Y — XB — + a + d + {1 — a)j^N'^, mn — /3F + A5 — (1 — o,)^Y. 

The variational matrix Mi at equilibrium point Ei is given by 

f-{u + a-\-d) 0 0 0 \ 

r 0 0 

0 s — So + 5^ 0 

1 0 ° -6oJ 


Ml = 


—a 

Si 

0 


Here two characteristic roots of the above matrix are positive implying that Ei is unstable. 
The variational matrix M2 at equilibrium point E2 is given by 

/-{XB* + u + a + d) XB* 0 0 \ 

M = ^ ^ ^ 

^ Si 0 -(s-So + (Jfii) 6B' ■ 

Vo I 0 “ -(io / 

Here two characteristic roots are -60 and -(s - Si +5^), and the other roots are given 
by the following quadratic 


+ {XB* -\-u + a + d-r}ip- {XB* + i^ + a + d)r + XB*a = 0. 
Hence using the Routh-Hurwitz criteria, the equilibrium E2 is stable if 


XB*+u + a + d>r and XB* ^- > (u + a + d) 

r 

and unstable if 


XB* -\-v + Oi + d<r or XB* 


(q! — r) 


r 


< (iz + ft + rf). 


The variational matrix M3 at equilibrium point E3 is given by 


M, 


! PK — {zz + o + d + (1 — a)r} 0 XK 

-a -r 0 

5l 0 (s- So + (5^2^ 

V 0 i 0 ^ 


0 \ 

0 

0 

“5o/ 



3.3 SIS Model with Logistic Growth 


79 


Two characteristic roots of the above matrix are —Sq and — r, the other roots are given 
by the following quadratic equation, 


■ 0 ^ — — (zz + a + d+ l — ar) + s — so + 5( 


Qo + IK 
6o 


) ?0 


+ S S — So “b 


'Qo + IK^ 


-^^{l3K-{u + a + d+l-ar)}-XKsi = 0. 
By the Routh-Hurwitz criteria the equilibrium Es is locally stable if 

,3ii:-(l-a)r-(i. + c( + ci) + s-s„ + 5||^i^j <0, 
and {s ~ So + >5 \}{I3K — + a + d+1 — a r)} > XKsi, 


are satisfied and it is unstable if 


pK - {1 - a)r -{u + a + d)+s-so + 6 > 0, 


or {s - So + ^ ^ - {u + o: + d + 1 — a r)} < XKsi, 

and £'4 exists. Clearly both conditions are not satisfied simultaneously, hence £3 is 
unstable. 

The variational matrix M4 at equilibrium point £4 is given by 

f-ipY+^) pY + XB - {1 - a)j^Y X{N-Y) 


Mi 


V 


—a 

Si 

0 


r-|iV 

0 

I 


0 


0 \ 
0 


(f + ^) 

0 -5oJ 


The characteristic polynomial in this case is given by 

Ip * + 03 0^ + CI2 0^ "b <^10 + Oq — 0) 

where 

r - XR/Vl r Or- .-i f s - 

+ 5o > 0, 


o-z 


0-2 



+ 

2r - 
-r + —N 

+ 


y 

K . 


B \ 


— <^0 I 


,,, XBN 2rN sB s,r1 XBN\ 

PY + -r+-- + -^}+{pY+ — ^ 

Y K L B } \ Y J 


2rN sB siF 
- r + — + -b- 
K Lb 
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gure 3.27: Variation of infective population with time for different intrinsic growth rate 
bacteria population. 



0 5 10 15 20 25 30 35 40 45 5 ( 

Time in years ^ 


Figure 3.28: Variation of infective population with time for different growth rate of 
bacteria population due to infective human population. 
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Figure 3.29: Variation of infective population with time for different growth rate of 
bacteria population corresponding to environmental discharges. 



Figure 3.30: Variation of infective population with time for different carrying capacity of 
bacteria population. 
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Figure 3.31: Variation of infective population with time for different 1 . 



Figure 3.32: Variation of infective population with time for different growth rate of 
human population. 
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3.4 Conclusions 

In this chapter, the SIS model for bacterial infectious diseases, like tuberculosis, typhoid, 
etc., caused by direct contact of susceptibles with infectives as well as by bacteria is pro- 
posed and analyzed. When demography for human population is constant immigration, 
it is shown that the endemic equilibria is globally stable when the cumulative rate of 
environmental discharges conducive to the growth of bacterial population is a constant. 
When the cumulative rate of the environmental discharges is a function of total human 
population, the endemic equilibria is also locally and globally stable under same condi- 
tions. The later result is shown by computer simulation. For logistic growth of human 
population, in each case the nontrivial equilibrium is locally asymptotically stable under 
certain conditions. By simulation it is shown that nontrivial equilibrium is globally sta- 
ble under local stability conditions for that set of parameters. It is concluded from the 
analysis that the growth of bacteria caused by conducive environmental discharge due 
to human sources increases the spread of the infectious disease. It is noted that if the 
population increases either by immigration or logistic growth, the spread of the bacterial 
disease further increases and it becomes more endemic. 


Chapter 4 


Modelling the Spread of Malaria: 
Environmental and Demographic 
Effects 


4.1 Introduction 


It is well known that the spread of malaria is governed by the following factors: (i) the 
density of the human population and its rate of growth, (ii) the density of the mosquito 
population and its rate of growth, (iii) various environmental factors sucli as rain, tem- 
perature, humidity and so on, (iv) ecological factors such as vegetation, biomass, cattle 
population, etc. and (v) geographical factors. In most of the tropical countries including 
India, the emergence of malaria has taken place and it has become endemic in the North- 
Eastern part of India, where this disease is spread by a lethal parasite called Plasmodium 
Falciparum. Although there have been several experimental studies related to surveys of 
malaria in different regions (Sharma 1991, 1998, Das 1991), the study of the spread of 
disease using mathematical models by considering the factors mentioned above has not 
been conducted, particularly when the densities of the human and mosquito populations 
are variable. However, in the case of these populations being constants, Bailey (1979) has 
given a simple mathematical model by considering criss-cross interaction between female 
mosquitoes and the human population. Some other studies for the spread of malaria have 
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also been conducted by proposing both stochastic and deterministic models (Radcliffe 
1973, 1974, Dietz et al. 1974, Molineaux et al. 1978, Bailey 1979, 1982). It has been 
suggested that the malaria model and the Gonorrhea model are very similar in mathe- 
matical structure (Nallaswamy and Shukla 1982). In those models the environmental and 
ecological factors have not been taken into account, though these factors play significant 
roles in the spread of malaria. In this chapter, therefore, an SIS non-linear mathematical 
model is proposed and analyzed to study the effects of household and other environmen- 
tal discharges on the ground, such as waste water, food stuff and so on, which are caused 
by the human population. It is assumed that the density of the mosquito population 
follows a generalized logistic model such that its growth rate decreases but its death 
rate increases as population density increases towards its carrying capacity with respect 
to the environment. It is further assumed that the growth rate of the mosquito pop- 
ulation increases as the density of cumulative environmental discharges increases. The 
human population density is divided into susceptible and infective classes. The mosquito 
population density is also divided into susceptible mosquitoes and infective mosquitoes. 
Our main focus is to investigate the effect of cumulative household discharges conducive 
to the growth of the mosquito population, on the spread of malaria by considering the 
following two types of demographics for the human population; 

(i) a population with constant immigration, 

(ii) a population with a logistic growth rate. 


4.2 Malaria Model with Immigration 


We consider here an SIS model, where the human population density Ni (t) is divided into 
two classes namely, the susceptible class Xi{t) and the infective class Ti(t). The mosquito 
population density N 2 {t) is divided into the susceptible class X 2 {t) and the infective class 
l 2 (t)- Keeping in view the above and by considering the criss-cross interaction of the 
mosquito population with the human population, a model can be written as follows: 
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Xi = A-diXi- Pi Xi ¥2 + 1^1 Yi, 

= Pi Xi ¥2 — {I'l + Q!i + rfi) ¥1, 

= A-diNi-aiYi, 

Z2 = {b2-a'~N2)N2-{d2 + {l-a')^^2}X2-p2X2Yi-a2X2 + S2N2E, 

K2 -^2 

¥2 = p2X2Yi-{a2 + d2 + {l-a')Y^N2}¥2, 

N2 = ^2 N2 — 0:2 iV2 + <52 X2 E, 

E = Q{Ni)- 5 oE, 

Ni = Xi + Fi, N2 = X2 + Y2, 

Xi{ 0 ) = Xio > 0, yi(0) = Yio > 0, X2(0) = X20 > 0, ¥20 > 0, £;(0) = /?o > 0. 

In model (4.1), A is the constant immigration rate of the human population; di is the 
natural death rate constant; Pi is the interaction coefficient of the susceptible human 
with the infective mosquito population; ui is the recovery rate coefficient of the human 
population; ai is the disease related death rate constant; 62 and ria are the birth and 
the death rate constants corresponding to the mosquito population; r2 = 62 — <^2 is the 
growth rate coefficient of the mosquito population; K2 is the carrying capacity of the 
mosquito population in the natural environment; 0:2 is the death rate of mosquitoes due 
to control measures (r2 > 0:2); y32 is the interaction coefficient of susceptible mosquitoes 
with the infective human class; 62 is the growth rate coefficient of the mosquito population 
due to the environmental discharges of cumulative concentration E; Q is the cumulative 
rate of environmental discharges which is human population density dependent; Sq is its 
cumulative depletion rate and 0 < a' < 1 is a constant (Gao et al. 1992), which governs 
the logistic birth and logistic death of the mosquito population. 

We analyze the model (4.1) for the following two cases ; 

(i) the rate of cumulative environmental discharges Q is a constant, and 

(ii) the rate of cumulative environmental discharges Q is a function of human population 
density. We consider the form of Q(7Vi) as Q{Ni) = Qq-{- IXi, where / is a constant. 
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4.2.1 Case I: Q = Qa 3l Constant 


Since Xi+Yi — Nx and ^2 + 12 = in this case it is sufficient to consider the following 
subsystem of the model (4.1), 

ii = MNx-Yx)Y^-{ux^ax+dx)Yx, 

Nx = A-dxNx-axYx, 

Y2 = ft (iVj - n) 11 - {a, + * + (1 - a')^N2j Y^, 

N2 = ^2 ^2 (1 ~ ^^) ~ (^-2) 

E = Qa — do E. 

We see that the region of attraction of the system (4.2) 

T = ((n.JV„y2,JV2,£:):0<yi<JV. < 4 , 0<l;<iV2<:^fr2-a2 + ft%'), 

t ax r2 \ do / 

0< E<^\, 

is positively invariant and all solutions starting in this region T stay in T. The continuity 
of right sides of (4.2) and their derivatives imply that a unique solution exists (Hale 1969). 
Since the system (4.2) is autonomous, the effects of N2 and E on the spread of malaria 
can be qualitatively studied by taking their asymptotic values as t — > 00 in the last two 
equations of (4.2). Thus we have 

lim sup E{t) = ^ = E (say) 

t-AOD Oq 

Ko - - 

and if A^2(0) > 0 lim supiV2 = — (’’a — 0:2 + 62E) = N2 (say). 

t^oo T2 

Here it is noted that N2 increases as Qa and K2 increase or as a2 decreases showing the 
effect of the growth of mosquito population. Now it suffices to study the global behavior 
of the system (4.2) by the following system of equations: 

= Pi{Nx-Yx)Y2-{vi+ax^dx)Yx, 

Nx = A-dxNx-axYx, 

— h {X2 — Y2) Yx — {0:2 + ^2 + (1 — l 2 - 


(4.3) 
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We note here that h>0 and also 61&2 - &3 > 0- Hence by the Routh-Hurwitz criteria 
the equilibrium Ei{Yi, Ni, Y 2 ), if it exists, is locally asymptotically stable. 

Nonlinear Analysis and Simulation: Before proceeding to simulation, we first analyze the 
model (4.3) with the disease related death rate ai = 0. In this case the model (4.3) 
reduces to 


Yi = ^i{Ni-Yi)Y2-{vi-¥di)Y^, 

Ni = A-diNu 

Y 2 = h{N2 - Y 2 ) n - \a2 + da + (1 - ^^ 2 - (4.5) 


The equilibria of the system (4.5) are El ^0, o) and E^iYi, iVj, Va), where 

Y - A NxY 2 at _ a y - A AM N 2 - (r/i + di){a 2 + A + (1 - qQ^A} 
l3iY2 + ui+di^ ^ di A[A M + {tt 2 + A + (1 - a')|^iV 2 }] 

(4.6) 

The second equilibrium exists if 


AM -^2 

(i/i + A) {oi2 + A + (1 - aO^A^a} 


= i?oo (say) > 1 . 


(4.7) 


We note from (4.6) that Y 2 increases as iVa, i.e. as Qa increases. This leads to an increase 
in Yi and for large Fa, Y N\. It is easy to see that the equilibrium £'o(0, ^,0) is 
unstable under (4.7) and is stable otherwise. Here Roo is the threshold for the system of 
equations (4.5). 

We can also prove that the equilibrium E^{Yi, Ni, Fa), if it exists, is globally stable, by 
using the following Liapunov function and choosing suitable values of ki and A, 


= 5 (U - Y,f + (N, - (K, - 

Since our system (4.3) is bounded by the system (4.5) in the region T, using a comparison 
theorem (Lakshmikantham and Leela 1969), we conclude that the solution of system (4.3) 
is bounded by the solution of (4.3) with oi = 0 . Thus we conjecture that the equilibrium 
El may be globally stable. This result is illustrated by integrating system ( 4 . 3 ) using 
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Figure 4.1: Variation of infective population Yi with susceptible population Xi. 

the fourth order Runge-Kutta method in the region T. The following data set is used 
(Greenhalgh 1992), 

Pi = 0.00000025 = Pi, = 0.012, ai = 0.0005, A = 10, di = 0.0004, aj = 0.045, r2 = 1 

di = 0.02, a' = 0.999, Ki = 1000000, Qa = 20, Jo = 0.001, Si = 0.0000002. 

Under the above parameter values, the equilibrium point Ei is found as: 

Yi = 4690.018, Ni = 19137.477, Yi = 16750.691. 

The simulation is performed for different initial positions 1, 2, 3 and 4 as given below, 

1.. .yi(0) = 5610, Vi(0) = 24092, ^ 2 ( 0 ) = 2180. 

2 . . . Fi(0) = 4610, Vi(0) = 16092, ^ 2 ( 0 ) = 1180. ; 

3 . . . Fi (0) = 3222, Ni (0) = 12000, ^ 2 ( 0 ) = 1000. 

4.. .yi(0) = 5300, iVi(O) = 15000, 3 ^ 2 ( 0 ) = 90. 

In Fig. 4.1, the infective population is plotted against the susceptible population which 
IS consistent with the solution tending to Ei if the starting point is not Eq. Thus the 
equilibrium point Ei appears to be globally stable if we start away from Eq. Also in Fig. 
4.2, the infective population is plotted with time for different A and Qa and we note that 
Yi increases as A or Qo increases showing that the spread of malaria increases' and it 
becomes more endemic due to increased immigration and household discharges. 





94 


Modelling the Spread of Malaria: Environmental and Demographic Effpr-to 



Figure 4.2: Variation of infective population Yi with time for different immigration rates 
of human population and the rate of cumulative environmental discharges. 

4.2.2 Case II: Q = Qo + INi 


In this case using Xj + Fi = Ni, X 2 + ¥2 = N 2 and Q = Qo + IN], the model (4.1) can 
be written as follows, 


Yi = 

1 

1 

g 

t— 1 


+ ai +di) Yi, 

Ni = 

A — diNi — Oil Yi, 



F 2 = 

i 

1 

{an 

+ ^2 + (1 — F 2 

N 2 = 


an 

N 2 + 62 N 2 E, 

E = 

Qq + 1 Ni — 5 q E. 



We see 

that in this case the region of attraction is 


r = {iY^,NuY2,N2,E):0<Yr<Ni<^, 

^ di 


So j ’ 

and the model is well posed. The result of equilibrium analysis is stated in the following 
theorem. 


0<Y2<N2< 


r2 


r 2 




,0< E< 
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THEOREM 4-3 There exist the following three equilibria, namely 


Qq 


(,)E, 0, 5-, 0,0,— 


I A 
di 


(ii) E2 0, 0, iV2, 

“1 


Qo + 


di 


5o 


, where N 2 


r2 


T’2 ~ CK2 "1“ ^2 


'Qo + ij 


and (Hi) Ez ( Fi, Ni, F, N 2 , 


Qq^INx 


A^2^iV2 


[0:2 + ^2 + (1 - a')^N 2 ]{i'x + 0:1 + di) 


Ez exists if 

= R'q (say) > 1. 


Proof; The existence of either of the first two equilibria is obvious. The existence of the 
third equilibrium Ez is shown by the isocline method. Setting the right hand sides of the 
system of equations (4.8) to zero, we get the following set of equations, 


ho ^2 I ' (Jo 


(4.9) 




P 2 N 2 YX 


P 2 Y 1 + a2 + ^2 + (1 “ oI)^N2 


„ _A-a,Y, „ A(JVi-yi)U 

7 iVi = , Yi = 


di 


yi + ai + di 


(4.10) 


Using these equations, we have the following quadratic in Yi, 


F{Yx)=PoY^+PxYx+P2 = ^, 


(4.11) 


where 


PO = ftft(l + ^ 


q ; i \ 521(^1 K 2 


Pi 

Po 


(Jodi r2 ’ 

+ ft'®2(l + fWi + (>^1 + ai + di) (& - (1 - 

di 0(fd\r2 di [ Oorfi J 

— { 1^1 + 0 : 1 + di) |q:2 + ^2 + (1 "• ■ 


From (4.11), we have the following, 

(i) F(0) = Pi/32-^^2 — {<^2 + (^2 + (1 — Q0^-^2}(i^i + ai + di), 

(ii) F (^) < 0, and 

(iii) F(oo) > 0. 
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So if F(0) > 0, then there exists one root say between 0 and and another root 
between and oo. Corresponding to this Yi, the values of Nu N 2 and E are 

ai+di 

determined using ( 4 . 9 ) and ( 4 . 10 ) and we get a nontrivial equilibrium point under the 
following condition 



{a2 + d2 + {l-a')j^N2}iJ^i + o:i + di) 

Here Rq is the threshold for the system (4.8). 

4.2. 2.1 Stability Analysis 

Here, we present the stability analysis of the three equilibria and tlic local stability results 
of these equilibria are given in the following theorem. 


THEOREM 44 The equilibrium Ei is unstable, the equilibrium E 2 is stable if R'q < 1 , 
otherwise the equilibrium E 3 exists and is locally asymptotically stable provided 


C 4 

1 

0 


02 Co 
C3 Cl 


> 0 , 


C4 C2 


C4 C2 Co 0 

1 C 3 C 2 0 

0 C4 C2 Co 

0 1 C3 Cl 


where Cq, Ci, C 2 , C 3 , and C 4 are given explicitly in the proof of the theorem. 


Proof: Let Mj be the variational matrices corresponding to equilil)riurn points Ei for 
i = l,2,3. 


Ml = 


/-(j/i +01 +di) 
-Oi 
0 

0 

0 


0 

-di 

0 

0 

I 


01 A 

di 

0 


—{ 0 C 2 + ^ 2 ) 

0 

0 


0 

0 

0 

r 2 - Q 2 + 62 
0 




0 \ 
0 
0 

0 


-SoJ 


M 2 


/ ~(i^i +01 + di) 0 


-oj -di 

ft N 2 0 

0 0 

0 I 


01 A 

di 

0 


-{a2 + ft + (l -a')^iV2} 


0 

0 


0 

0 

0 

K 2 

0 


0 \ 
0 
0 

0 

-SoJ 
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M.= 


[-{Pi Y2 + i'i+ai + di) Pi Y 2 Pi{Ni — Fi) 0 0 \ 

70^1 -di 0 0 0 

P2{N2-Y2) 0 mss ;52Fi-(l-a');^F2 0 

0 0 0 So N‘ 


V 


0 I 0 0 -So J 


where ttiss — ~[P2 Fi + 0:2 + <^2 + (1 ~ ci^)^iV2]. 

From Ml it is clear that Ei is unstable. From M2 it is clear that E2 is also unstable 
under following condition, 

/5l/52^iV2 


{cx2 + d2 + {I — a')j^N2}{ui + ai + di) 

The characteristic polynomial corresponding to the matrix Ms is given by 


+ 041 /''* + Cs + C 2 + Cl ^ + Co = 0, 


where 


^2 .V 


C4 — ' 7 p ~^2 + Pi Y2 + 5 q + ui + ai + 2 di + P2 Yi + 0:2 + d 2 + {1 — a')-^N 2 i 

^2 K 2 


C 3 


/?1 O'! + N2 

rV2 


+^0 


|(5o + Pi Y 2 + i2i + cxi + 2di + P 2 Y 1 + a2 + d2 + {1 — a')-^N2'^ 


ui + PiY2 + ai+ 2di.+ P2 Yi + a2 + d2 + {l- a')^N2 

K 2 


+(di + P 2 Yi){Pi Y 2 + vi Q-i -\- di) 


+ 


q; 2 + ds + (1 ~ cOt^-^2 

K 2 


PlY 2 +di 


r2 


P 2 Fi + q; 2 + ds + (1 ~ o!)-—N2 

■tio 


C2 


Pi Y2ai 
r2 N 2 

Ko 


do + P2Y1 + CX .2 + d2 — a!) T7"-^2 + 

rt2 Jti2 

■^do 1^1 + ai + 2 di + Fi + 0^2 T ds + (1 — of) 


To '• 


+(di + P2 Yi){Pi Y2 vi + Oil -{■ di) 

+ |q; 2 + ds + (1 — P 1 Y 2 di P 2 Y 1 + q; 2 + d2 + (1 — a!) 

f „ . No - 

+do I (di + P2 Yi)(Pi Y2 + 1^1 + c^i + di) + r 2 —PiY 2 


+ |a2 + ds + (1 — a') ^-^2 \ Pi Y2 + di 


^N2 

K 2 \ 


r2 


P 2 + 0:2 + d2 + (1 — a')—N2 

A 2 


+di^2 Yi{Pi I 2 + 1^1 + ai + di) + di |q; 2 + d2 + (1 — o!)-^N-;^ Pi Y 2 , 
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Cl = ^iV2{5o(c^i+^2 5^i)(A^2 + i^i+Q=i+‘^i) 

+i„ |q 2 + d2 + {l- A Y 2 + So di [a y. + Q2 + a + (1 - 

+<ii A y(A y+i'i+ai +‘^i)+‘*' {<12+ A + ( 1 A a} 


+^o di p 2 Yi {Pi Y 2 + iyi + ai + di) + Sodi | o;2 + 4 + (1 «') ^ 2 1 Pi Y 2 , 


+ 011 P 1 Y 2 

K 2 


+ <^0 {^P 2 Yi + Q!2 + <^2 + (1 a) 


(^P 2 Yi + 0;2 + 0^2 + (1 — Y) 


,^^2^2 


Ko 


Co = 


II 

iiir2 


iVa^o rfi p2 YliPi Y2 + vi^ai + di) + Y ^2<^o c?i {a2 + ^2 


+(1 — a')^iV2| /3i y2 + ^ -^2 <^0 


P2 Yi + Q!2 + (^2 + (1 ” (l')-j^N2 


Pi Y2ai 


+52 N 2 <1 P 2 Y 1 - (1 - a')^^2 [ Pi{Ni - Yi)aiL 


Then for the roots to have all real parts negative, it is necessary and sufficient that the 
following inequalities are satisfied, 

1 C 4 C 2 Co 0 0 


C4 > 0 , 


C4 C2 
1 C3 


>0, 


C4 C2 Co 

1 C3 Cl 

0 C4 C2 


> 0 , 


C 4 C 2 Co 0 

1 C 3 C 2 0 

0 C 4 C 2 Co 

0 1 C3 Cl 


> 0 , 


1 C;{ C2 0 0 

0 C 4 C 2 C{) 0 

0 1 C3 Cl 0 

0 0 C.1 0 Co 


> 0 . 


The first two conditions are always true. If the next two inequaliti(!s arc satisfied then 
so is the fifth as Cq > 0. Hence the equilibrium point 2^3 is locally asymptotically stable 
under the conditions mentioned in the theorem. 


Nonlinear Analysis and Simulation: As in Case I, it is noted that the solution of system 
(4.8) is bounded by the solution of its corresponding system with 0!i = 0, which is globally 
stable. Hence, as before, we conjecture that the equilibrium of (4.8) may be globally 
stable in the interior of the region of attraction. To illustrate this global behaviour of 
nontrivial equilibrium point E 3 and to show the effects of various parameters on the 
spread of malaria, we use simulation. 

The system (4.8) is integrated using the fourth order Runge-Kutta method by using the 
same set of parameters as in Case I with Qo = Qa and an additional parameter value 
I = 0.0005, which satisfy the local stability conditions. 
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Susceptible popualtion 


Figure 4.3: Variation of infective population Yi with susceptible population Xi. 


The equilibrium E 3 is calculated in this case as, 

Yl = 4701.9, Ni = 19122.2, ¥2 = 16825.0, N 2 = 960912.2, E = 29561‘.0. 

In this case also simulation is performed for different initial positions as shown in Fig. 4.3 
implying the global stability of E 3 . Also in Figs. 4.4-4. 7, we have shown the variation of 
Yi with time, for different r 2 , <52, Qq and I respectively. We note that with the increase of 
any of these parameters, the infective population increases showing the effect of household 
discharges on the spread of malaria. The effect of immigration on the spread of malaria 
is the same as in Case I. 




102 


Modelling Snr...d of Malaria: Environmental and Demographic Effects 


4.3 SIS Model for Malaria with Logistic Growth of 
Human Population 

Now let us consider an SIS model where the human population growth is logistic so that 
both the birth and death rates are density dependent in such a manner that the birth 
rate decreases and the death rate increases as the population density increases towards 
its carrying capacity (Gao and Hethcote 1992). As before the human population density 
is divided into susceptibles and infectives. Also the mosquito population is divided into 
susceptible mosquitoes and infective mosquitoes. It is assumed that mosquito population 
density grows logistically in the environment and its growth rate increases due to cumu- 
lative concentration of various discharges from household sources into the environment. 
Keeping these in mind, a mathematical model is proposed as follows: 

= [h - ari ^]Ni - P1X1Y2 - [di + (1 - a)ri ^]Xi + Yy , 

Ki -^^1 

Ni 

Yi = ^iX\Y2 - [t^i + Oil + di -1- (1 — a)ri— jyi, 

iVi = ri[l - - aiKi, (4-13) 

Ki 

X2 = (62 - a'~N 2 )N 2 - {d 2 + (1 - a')^N 2 }X 2 - ^2X2^ ~ 0:2X2 + 52X2^, 

K2 ^2 

Y 2 = ^2X2 yi-{a 2 + d2 + (l-a')^iV2}r2, 

N 2 — ^2X2(1 — -^) — 0 L 2 N 2 + 62 N 2 E, 

K2 

E = Q{N)- 6 oE, 0 <a < 1, 0 < a' < 1, 

Xi{0) = AAo > 0, Fi(0) = Kio > 0, X2(0) = X20 > 0, 72(0) = 72o > 0, E{0) = Eo>0. 

Here bi and di are natural birth and death rates; rj = 61 — di is the growth rate constant 
and Ki is the carrying capacity of the human population in the natural environment. 
All other parameters have already been defined in Section 4.2. For 0 < a < 1, the birth 
rate decreases and the death rate increases as Ni increases to its carrying capacity Ky. 
When a = 1, the model could be called simply logistic birth model as all of the restricted 
growth is due to a decreasing birth rate and the death rate is constant. Similarly, when 
a = 0, it could be called a logistic death model as all of the restricted growth is due to 
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an increasing death rate and the birth rate is constant. 

In this section also, we analyze the following two cases: 

(i) the rate of cumulative environmental discharges Q is constant, 

(ii) Q is human population density dependent as described in Section 4.2. 


4.3.1 Case 1: Q = Qa 3i Constant 


In this case the model (4.13) reduces to the following form, 

ri = A(iVi-n)l"2-h+Q;i + di + (l-a)^]yi, 

Ni = ri (l - iVi - aiFi, 

y = ^2(^2 — 12)11 ~ |q ;2 + (^2 + (1 — 

N 2 = r2(l-^^N2-a2N2 + 62 N 2 E, 

E = Qa-5oE. 


Using the asymptotic values of E and N 2 , E = ^ = E (say) and 

N 2 = ^{^2 — 0:2 + <52^} = N 2 (say) in the above set of equations, the behaviour of 

(4.14) can be given by the following subsystem. 


l\ — Pi{Ni — Yi)Y2— + ai + 0?! + (1 — a) 


riiVi 


K, J 




Ni = n (1 - iVi - aiFi, 

Y 2 ~ ^ 2(^2 — Y2)Yi — |o;2 + 0^2 + (1 — 

We note here that the region of attraction of the system (4.15) is given by, 


(4.15) 


T = {(yi,iVi,y2) :0 <Yi < Ni < Ki,0 <Y2 < N 2 } . 
The results of the equilibrium analysis is stated in the following theorem. 


THEOREM 4-5 There exist the following three equilibria, namely 

(i) £^i (0,0,0) , (ii) £' 2 ( 0 , iiri,0) and (Hi) Ez{Yi,Ni,Y 2 ). The third equilibrium Ez exists 



Figure 4.8: Existence of equilibrium point. 


Plp2KiN2 


a, + d2 + {l- + (1 - 


R' > 1 . 


Proof: The existence of Ei or E 2 is trivial. The existence of tlie nontrivial equilibrium 
Ez can be proved as follows by using the isocline method. From (4.15) wc have, 

(4.16) 


Oi V Ki, 


and 


1^1 = 


P 1 P 2 N 1 N 2 - \ a 2 + d 2 + {l- a') 1 1 1^1 + Oi + di + (1 “ n) 


riN\ 


P 1 P 2 N 2 + h 1^1 + Oi + + (1 — a) 


K, 


(4.17) 


Clearly (4.16) is a parabola passing through (0, 0) and (FTbO) with vertex at (^, ^)- 
From equation (4.17), we have the following: 


for Ni = 0, Fi = 


— {q! 2 + (^2 + (1 ■“ + Oi + di) 

ih ^ 

P 1 P 2 N 2 + ^ 2(^1 + Oi + di) 


< 0 , 


(z^i + oi + di){a2 + c!2 + (1 — ci!)-^N2} 

for Yi = 0, iVi = — ^ — rr:-_ ^ ^ — ^TT-^ it^ > F 


/3i/32-^2 ~ {<^2 + <^2 + (1 — a')^iV2}(l — a)- 

K 2 Ki 
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for Ni = Ki, Yi 


M 2 K 1 N 2 - {a 2 + ^2 + (1 - + ai + di) + (1 - a)n} 

Jti?. 

P 1 P 2 N 2 + /32{^i + Q:i + + (1 — a)?"!} 


which is positive provided 


0 iP2 K^N2 

{q;2 + ^2 + (1 — a')-^N2}{i'i + Q!i + di 4- (1 — a)ri} 

Also the slope of (4.17) is given by 


^y- P1I32N2 - {0:2 + d 2 + {l- a')~N 2 }{l - a)-^ - /? 2(1 - o)-p-Fi 

Pil32N2 + l32{pi + ai+di + {l-a)^Ni) 

Ki 

which is positive using (4.17) and i^o > 1 for Ni > 0. Thus, plotting (4.16) and (4.17) in 
Fig. 4.8, it is clear that there exists a unique nontrivial equilibrium point £'3(^1, Ni, Y2). 
Here Rq is the threshold for the system (4.15), which we will also see from the instability 
condition of the disease free equilibrium point. 


4.3.1. 1 Stability Analysis 

Now we present the stability analysis of these equilibria. We state the local stability 
results of these three equilibria in the following theorem. 


THEOREM 4-6 The equilibrium Ei is unstable, the equilibrium E2 is stable if Rq < 1 
and unstable if Rq > 1 in which case the equilibrium E3 exists. The equilibrium E3, if it 
exists, is locally asymptotically stable when as > 0 and 0102 — 03 > 0, where oi, 02 and 03 
are given in the proof of the theorem. 


Proof: The variational matrix at {Yx,Ni,Y2) corresponding to the system (4.15) is 


(-{^xY2 + ux^ax+di^{l-ay-^) P 1 Y 2 l3i{Nx-Yi)\ 


■M 


-ai 

UN2 - Y2) 


Ki 


0 

^33 


V 
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where 77133 = —[^2^1 + 0!2 + ^2 + (1 " 

At the equilibrium point ( 0 , 0 , 0 ), the variational matrix Mi is given by 


/ + Q!x + 0 


Ml = 


V 


-ai 

^2^2 


0 \ 

f'x _ 

0 — {q!2 + (^2 + (1 ~ / 


Prom iW^i, it is clear that it has one positive characteristic root 7i implying instability of 
this equilibrium point. 

At equilibrium point p2(0, Ai, 0 ) the variational matrix M2 is given by 


M2 = 


( — {z^i + Q!i + cfi + (1 — o)ri} 0 

-a_i -n 

^2^2 0 


PxKr 
0 

>2+d2 + (l-0/)f;iV2}, 


Clearly one characteristic root of the above matrix is -ri and other two roots are 


■0^ + {ui + Q:i + di + (1 — a)ri + 0:2 + (^2 + (1 ~ u/) Y7-iV2}V'’ 

J-^2 

+{q ;2 + (^2 + (1 ■“ + O'! + di + (1 — 0,)ri} — /5i/?2 A 1 A 2 = 0. 

A2 

Using the Routh-Hurwitz criteria, we get the condition for stability of the equilibrium 
E2 as Rq < 1, otherwise this equilibrium becomes unstable and the equilibrium point Ez 
exists. 

At equilibrium point Ez the variational matrix Mz is given by 


Mz 


^ ~{Pi^2 + t'l + Q!i + di + 1 — a ^{¥2 

-^1 ri - l^iVi 

I32{N2-Y2) 0 


Pi{Ni-Y:)\ 

0 

mzz J 


where 77133 = -[P2Y1 + {0:2 + ^2 + (1 - a')j^N2}]. 

The characteristic polynomial corresponding to the above matrix is given by 


+ OiV>^ + a2ip + 03 = 0, 
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a-z 


IL 

Kx 


{2N\ Ki) i + Oil + 4- (1 — a) ■ ^ - + ^2^1 + Q!2 + c?: 

V Ki 


+ cKi + di + (1 — a) 

+ axPxY2 + Pxh{Nx - Yi)(N2 - Y 2 ), 


X |/32^1 + 0:2 + + (1 — o !)-^ N 2 

K2 


as = 


H 

K 


{2Nx - Kx) p2Yx + ai + rfi + (1 - a)^| 

+;^ (2iVi - Kx) \a 2 + ^2 + (1 - a0;^^2} 

^zYx + |o;2 + £^2 + (1 — • 


It may be noted here that ai > 0, hence from the Routh-Hurwitz criteria, Ez is locally 
asymptotically stable if as > 0 and 0102 — 03 > 0. 

Remark: It is noted that for iVi > ^, 03 > 0 and 0102 — 03 > 0. Hence Ez is locally 
asymptotically stable under the assumption Nx> 


Nonlinear Analysis and Simulation: 

We first note that the system (4.15) is globally stable when the disease related death rate 
Q:i is zero. Since system (4.15) is bounded by its corresponding system (4.15) with ai = 0, 
using a comparison theorem (Lakshmikantham and Leela 1969), it is concluded that 
solution of (4.15) is bounded by the solution of (4.15) with o;i = 0. Thus, we conjecture 
that the system (4.15) is globally stable in the interior of the region of attraction. To 
support this result, system (4.15) is integrated by the fourth order Runge-Kutta method 
using the following set of parameters in the simulation, which satisfy the local stability 
condition. 


/3i = 0.00000022 = /?2, yx = 0.012, o = 0.3, ocx = 0.0005, dy = 0.0004, 

rx = 0.0003, Kx = 50000, a2 = 0.045 , r2 = 1, o' = 0.999, 

d2 = 0 . 02 , S 2 = 0 . 0000002 , Qa = 20 , 5o = 0 . 001 , K 2 = 1000000 . 

The equilibrium values for this set of parameters are given by 


Yx = 7475.6392, Nx = 26424.8209, Y = 23331.6445. 
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Figure 4.9: Variation of infective population with susceptible population. 



0 10 20 30 40 50 60 70 80 90 


Time in years ^ 

Figure 4.10: Variation of infective population with time for different growth rates of 
human population. 





4.3 SIS Model for Malaria with Logistic Growth of Human Population 


109 


In Fig. 4.9, we have plotted the infective population against the susceptible population 
for different initial positions 1, 2, 3 and 4 and from the solution curves, we infer that 
system is globally stable about the endemic equilibrium point (l^i, Ni, i^,), provided 
that we start away from Ei and E^. From Fig. 4.10, we note that when the growth rate of 
human population r i increases, the endemic infective human population increases. Also 
when Qa increases, the endemic infective population increases. 

4.3.2 Case II: Q - Qo + INi 


In this case, since Xj + Fi = Ni and X 2 + Y 2 = N 2 , the system (4.13) reduces to the 
following form, 


Fi = 

Ni = 

Y 2 = 

N 2 = 

E = 


MNx - Y,)Y2 - 

Ni 


1^1 + Q!l + dx + (1 ~ o) 


riNi 


Ki 


Yu 


niVx 1 - 


Ki 


— OiYi, 


I32{N2 - Y2)Y, - 1^2 + d2 + (1 - ^2- 

7-2X2 (1 - - 0 ! 2 X 2 + S2N2E, 

Qo + ZXx - SoE. 


(4.18) 


In this case the region of attraction is given by, 

Ti = l{YuNi,Y2,N2,E) : 0 < ri < Xi < K^O < Y 2 < N 2 < ^2,0 < E < 


Qo -F IKx 


The results of an equilibrium analysis is stated in the following theorem. 


THEOREM 4-7 There exist the following five equilibria, namely 

(i)B,(0,O,O,O,^), (ii) Bi ( 0 , Ki, 0 . 0 , 22^) , 

(iii) E, (o.O.O.iV,, where N 2 = ■§ {r, - a, + , 

(iv) Et ( 0 , J!-,, 0, JVX, 22^), where iVJ = a {r, - o, + 5,^^} 
and (v) Ez(Xx,Ni,Y 2 ,N 2 ,E). E 5 exists provided 

^ i . = (say) > 1. 

|q; 2 + £^2 + (1 — + ai + di + (1 — o)ri} 
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4. 3. 2.1 Stability Analysis 

We now present the stability analysis of these equilibria. The local stability results of 
these equilibria are stated in the following theorem. 


THEOREM 4.8 The three equilibria Ei, E2 and Ez are unstable. The equilibrium £'4 
is stable if R'qq < ^- V Ko < H is unstable and the equilibrium Ez exists. The fifth 
equilibrium £5, if it exists, is locally asymptotically stable provided 


ao > 0 , 


<Z 4 

1 


az 


> 0 , 


(X4 

1 

0 


02 

az 

04 


00 

01 

02 


> 0 , 


O4 02 Oo 0 

1 03 oi 0 

0 04 02 Oo 

0 1 03 O) 


where Oq, oi, 02, 03 and 04 ore given explicitly in the proof of the. theorem. 


Proof; The variational matrix M corresponding to the system ( 4 . 18 ) is given by 

/ 

n - ^ 0 


M = 


mil 

-Oil 

fi 2 {N 2 -Y 2 ) 

0 

0 


0 

0 

I 


A(iVi-yi) 

0 

mzz 
0 . 
0 


0 

0 


{l-a')f-Y 2 

mu 

0 


0 

0 

0 

S2N2 

-SzJ 


where mu = — |/3il2 + i^i + ai + di + (1 - 

mzz = - [P2Y1 + 0:2 + ^2 + (1 - a')jf^N2^ and m44 = r2 - 0-2 - 

The variational matrix Mi at equilibrium point Ei is given by 


/ -{ui + ai + di) 

-a 


Ml = 


V 


0 

0 

0 


0 

n 

0 

0 

1 


0 

0 

— { 0 l 2 + (^2) 
0 
0 


0 

0 

0 

r2-Q:2 + ^2^ 

0 


° \ 

0 

0 

0 

-So) 


Here two characteristic roots of Mi are positive. Hence this equilibrium point is unstable. 
The variational matrix M2 at equilibrium point £2 Is given by 


( — (pi + tti + di + (1 — o)ri) 

0 

PiKi 

0 

0 \ 

—a 

-ri 

0 

0 

0 

0 

0 

-( 0:2 + da) 

0 

0 

0 

0 

0 

ra - Q!2 + S 2 

0 

0 

1 

0 

0 

-Sol 
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Here one characteristic root, {r2 — a2 + is positive implying that this equilib- 

rium point is unstable. 


The variational matrix M3 at 

equilibrium point E3 

is given by 




/ 

-{ui -1-ai -l-di) 

0 

0 


0 

0 \ 



-Q_l 

n 

0 


0 

0 

M3 = 


P2N2 

0 

— ^0:2 -l- (^2 4- (1 — 

■o') ^-^2] 

0 

0 



0 

0 

0 


-itA 

S2N2 


V 

0 

i 

0 


0 

—^0 / 


Here also one characteristic root ri is positive implying instability of this equilibrium 
point. 

The variational matrix M4 at E4 is given by 


Mi 


/ — {2^1 4- Oil -h -t- (1 — a)ri} 0 0 iKi 

-ai -ri 0 

0 - [^2 + d2 + (1 - 

0 0 o' 

VO ; 0 


0 

0 

0 


Here the characteristic polynomial is given by 

(Qo + IK,) 


0 \ 
0 
0 


S2N^ 

0 -<io y 


(^1 + '0) ^^2 — 0:2 + S2- ^ ^ + 0^ {Sq + tp) {[ui -f Oil -f rfi H- (1 — a) Ti + Ip) 


X 1^2 4- d2 + (1 - ^');^A'2 + V'} - U -a 2 +S 2 ^°^ 


(^0 


= 0 . 


Clearly three roots are negative and other roots are given by the following quadratic 


.02 ^ ^ ^ ^ ^ + (^2 + (1 ““ j ^ 

+{2^i+Q!i4-di4'(l —u) Ti} |q;2 4- ^2 4- (1 — cl)~^N2'^—PiP2 — ~ — ^^2 ~ 0:2 4- S2 ^ ^ 

Hence this equilibrium point will be stable if 

(i^i-|-Q:i-[-ci!i4-(l— a) Ti) |a2 4- 0^2 + (1 ■" ^ A^2 ~ ^^2 — q;2 4- 62 ^ ^ 

i.e. J?()o < 1, otherwise if > 1 it is unstable and the fifth equilibrium point will exist. 


The variational matrix M5 at equilibrium point Es is given by 


/ mn 


M,= 


-ai 

ft(iV2 - 1^2) 


0 

V 0 


I3iY 2 — I — O’ j^y\ 0 

ri - ^Ni 0 0 

0 ^33 ^ 2^1 — (1 — 

0 0 -^iV2 

1 0 0 


0 \ 
0 
0 


62 N2 

-So J 
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m„ = -{ftn + 1^1 + ai + di + (1 - a) nf} and ™33 = + c<2 + d, + (1 - a’)^JVj] 

In this case, the characteristic polynomial is given by 


where 

a 4 = 

as : 


^iY 2 + 1^1 + ai + di + (1 - + ^2^1 


2 rx 


Ki 

r% ^-v . ’'2 ^-v 


Kx 


+a 2 + ^2 + (1 - + ^^2 + ( 5 o > 0 , 

{^xF 2 + .X + C.X + dx + (1 - a)^iVx} {|i(iVx - f ) 4- ^2^ 

+Q;2 + 4 + (1 - o!)^N 2 + -^^2 + ( 5 o| 


K2 


Ko 


+^{Ni - {P 2 Y 1 + a2 + d2 + {l- + ~N2 + <5o 


Ki 


+ 


r 2 


P2Y1 + q ;2 + (^2 + (1 — a ) N2 

J ^2 


r 2 


K2 

7’2 


^^W 2 +d„j + ^W 2 ,S„ 


+a. [m - (1 - ~ 


a 2 = 


|/5x^ + i^x + ttx + dx + (1 — a)— iVx| ( Ni 

+{l-a')^N 2 + ^N 2 + So 

i \2 J \2 


2 


'hYi + fva + (k 




J02^i + 0!2 + d2 + (1 — a!)-rr^2 

•h .2 


+ d„ 
A2 


+ 


P2Y1 + Q ;2 + d 2 + (1 — ( 2 ')-^N 2 

j ^2 


x^-^2 |do + ^x^ + t'x + Q;x + dx + (1 — 

+^■^2^0 |a^2 + i^x 4- ax 4- dx + (1 — a)-^Ni + ^iVi — ^ j | 


^1 .-V 


4-do l^xd^ 4- i^x 4 - ax -h di 4- (1 — cl)-^Ni 


P 2 Y 1 4- a2 4" d2 4- ( 1 ~ (i')~rr^^ 

l \2 


Ox 


4-ax |Ai2 — (1 — |^2^i 4- 02 4- d2 4- (1 — 0'')-^N2 4- j^A '2 4 60 

-PMNi - rx)(iV 2 - ^2) {||- {Ni - 

1^1 ^2 4- i^x 4- 0!i 4- dx + (1 - a)^iVx| ^ ^iVi - 

P 2 Y 1 -f 02 4- d2 4- (1 — a!)-~N 2 ( t^-^ 2 4- do 
>- A 2 J \A2 


+ ;^W 2 + io}, 


X 


-4 


/52Fx 4- 02 4- d2 -I- (1 — a')~^N 2 

xVo 
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x^iV2<5, 


'o [PiY 2 


+ I/I + Oil + di + (1 - a)^Ni + } 


+Y^N2So [PiY2 + + (1 _ a)^iVx} |l (iVx - ^ 

+ai[piY 2 -{l-a)^Yi'^ 


X I -j P 2 Y 1 + 0:2 + (i2 + (1 — a')-^N2 

J ^2 


— /^2 + ^o') + 

■^2 / .^2 


-A^ 2 (^i - yx)(iV 2 - F 2 ) (1^ (iVx - (;giV2 +5o) + ^iV25o} 

Qo = + Oil +cfi + (1 - a)^iVi| ^ ^iVi - /52Fi^iV2(5o 

+Q;i — (1 — P 2 Y 1 + q;2 + £^2 + (1 ~ (i')-^N2 -^N25o 

^ Ai j L il 2 J iV 2 

+Q:xA(iVi - yi)(52 I {^2^1 - (1 - «');^^2} N 2 


(^. - f 


{ 


0^2 + <^2 + (1 ~~ a^)-~N 2 
^2 


} 


^iV25o. 

A2 


By using the Routh-Hurwitz criteria, will be locally asymptotically stable if the 
following conditions are satisfied. 


04 > 0, 


OLA OL 2 

1 03 


> 0 , 


OLA 0 L 2 0-0 

1 03 Oi 

0 04 02 


> 0 , 


<24 02 flO 0 

1 03 Ox 0 

0 04 02 Oq 

0 1 03 Ox 


> 0 , 


QjA 

0>2 

ao 

0 

0 

1 


ai 

0 

0 

0 

cla 

02 

Qq 

0 

0 

1 

03 

ax 

0 

0 

0 

cla 

02 

ao 


> 0 . 


Remark: Here 04 > 0. Also it is noted that if iVx > the second inequality is satisfied 
and oo > 0. 


Nonlinear Analysis by Simulation: 

As before, we conjecture that the equilibrium E 5 is globally stable in the interior of the 
region of attraction. To illustrate this global stability behaviour of E 5 and to see the 
effect of various parameters on the spread of the disease, the system (4.18) is integrated 
using the fourth order Runge-Kutta method by taking the same set of parameters as in 
case I with Qq = Qa and an additional parameter / = 0.0005, which satisfy the local 
stability conditions mentioned above. 

The equilibrium value for this set of parameters is 

Yi = 1408.005, Ni = 47531.467, Y 2 = 1817.685, N 2 = 385501.258, E = 43765.735. 
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Figure 4.12: Variation of infective j>oi>n!ati(»n wuJi iM,jiiiiaHon. 


In this case also simulation is perforine<i fur flilihiuiif mifial pu atiun . from the so- 
lution curves (Fig. 4.12), we conclude that thi.s ctjuiUiu uam ituni! r. i'l"l>;il!y stable in 
the region of attraction Ti, provided that we .start awav liusii uf le-i funi equilibria. Also 
from Figs. 4.13-4.16, we see that a.s any of the parametei*.. t and I correspond- 
ing to the growth of the mosquito population density in« rea.se;,. the iiilei tivi* population 


increases. 
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Figure 4.13: Variation of infective population Yi with time for different growth rate 
coefficients of the mosquito population. 



Figure 4.14: Variation of infective population Yi with time for different growth rate 
coefficients of the mosquito population due to environmental discharges. 
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Figure 4.15; Variation of infective population Yi with time for different rates of the 
cumulative discharges. 



Time in years 


Figure 4.16; Variation of infective population Yi with time for different 1. 
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4.4 Conclusions 

In this chapter, an SIS model for malaria is proposed and analyzed by considering en- 
vironmental and demographic effects. The rate of cumulative environmental discharges 
are assumed to be either a constant or human population dependent. The cases of con- 
stant immigration and logistic growth of human population are considered. The growth 
of the mosquito population is assumed to be logistically varying. Threshold conditions 
for spread of malaria are derived in each case. If the threshold is greater than one, 
the nontrivial equilibrium is always feasible. In the case of constant immigration, the 
nontrivial equilibrium is shown to be locally asymptotically stable when the rate of cu- 
mulative environmental discharges is a constant and the same result is also true when 
the rate of cumulative environmental discharges is human population dependent, but 
under certain conditions. However, in the case of logistic one, the nontrivial equilibrium 
is locally asymptotically stable to small perturbations in both the cases of environmen- 
tal discharges with under certain conditions. By simulation it is also observed that in 
both types of demographics, these equilibria are in fact globally stable under local sta- 
bility conditions for the set of parameters considered. Thus, it is concluded that due 
to household discharges, the mosquito population can grow very large leading to fast 
spread of malaria. Also, if the rate of immigration increases or the growth rate of human 
population increases, the infective population increases and the disease becomes more 
endemic. 



Modelliiig the Spread of Malaxia with Human Reservoir: Environ- mental and Demographic 
^22^ Effects 

'—X2X2Z1 — (^2X2 + ^2X2Ej 

!> = p2X2Yi + X2X2Zi-S^a2+d2 + {l-a')^^N2]Y2, 

iVj = r2N2 ~~ ^2X2 + S2N2E, 

E = QiNi)- 5 oE, 

Ni — Xi -h Yi + Zi, N2 — X2 + Y2, 

:^i(0) > 0, Fi(0) > 0, Zi{ 0 ) > 0, X2(0) > 0, F20 > 0, E{ 0 ) > 0. 

In model (5.1), A is the constant immigration rate of the human population; rfi is the 
natural death rate constant; Pi is the interaction coefficient of susceptible human with the 
infective mosquito population; p\ is the recovery rate coefficient of the human population; 
Oi is the disease related death coefficient; 5 i is the rate coefficient corresponding to 
movenient of human population from the infective class to the reservoir class; r2 = 62-^2 
is the growth rate coefficient of the mosquito population; where 62 and d2 are birth 
and death rates constant corresponding to the mosquito population; K2 is the carrying 
capacity of the mosquito population in the natural environment; 0:2 is the death rate of 
mosquito population due to control measures, (0:2 < J'a); P2 and A2 are the interaction 
coefficients of susceptible mosquitoes with the infective and reservoir classes of the human 
populations respectively; 62 is the growth rate coefficient of the mosquito population due 
to environmental discharges of the cumulative concentration E\ Q{Ni) is the cumulative 
rate of discharges into the environment which depends upon the density of the human 
population and do is its natural depletion rate; 0 <a' < 1 is a constant, which governs 
the logistic birth and the logistic death of the mosquito population. 

We analyze the model in the following two cases, 

(i) the cumulative rate of environmental discharges <5 is a constant, and 

(ii) Q depends upon density of human population and is considered of the form 
Q(iVl) = Qo + INy. 
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5.2.1 Case I : Q = Qa-, a Constant 

Following the methods as discussed before, the behavior of the system of equation (5.1) 
can be studied by the following system, 


—Yi — Zi)Y 2 — {ui + ai + di 5i)Yi, 

Zi = SiYi — d\Zi, 

Ni = A — diNi — aiYi, (5.2) 

= /52(iV2-y2)n + A2(iV2-y2)^i-{a2 + d2 + (l-a')^iV2}y2, 

where 

N 2 = lim sup N 2 = — -{t 2 — 0-2 + 52 E} and E = lim sup E = 

r2 t -^00 ^ Sq 

We see that the region of attraction is 

Tx = |(yx,yx,iVi,y2) : o < yi + ^1 < Wi < ^, o < y2 < Wsj . 

The result of equilibrium analysis of the system (5.2) is stated in the following theorem. 
The proof is obvious. 


THEOREM 5.1 There exist the following two equilibria corresponding to system (5.2) 
namely, (i) Pi(0, 0,^,0) and (ii) P 2 {Yi,Zi,Ni,Y 2 ). The equilibrium P 2 exists if 


{0:2 + ^2 + (1 ~ + Oil + di + 5i) 


= Rs (say) > 1 


(5.3) 


where 


^ _ /3iA(^2 + ^2^)N2 — d^i/i + Oil + di + <?i){q!2 + ^2 + (1 — a')jf^N2} 

{p2 + A2|J-){A(o:i + <^i + di)N2 + di{ui + O'! + di + 5i)} 

_ A-aiYi ^ ^ (^2 + A2|^jiV2yi 

di ’ di {oi2 + d2 + 0- — a,')^N2} + {^2 + ^2 ^)Yi 

Here Rg is the threshold for the system (5.2). 



' 

No^ . j equilibrium. The local stability results of these equilibria are stated in the 
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stability Analysis 


Jjgcuss linear stability of both the equilibria and nonlinear stability of only the 



noU ^lieorem. 




e? 




^00 5.2 The equilibrium Pi is stable if Rs > 1 and unstable if Rs > 1 and the 
■^jn P2, If It exists, is locally asymptotically stable. 


^ rpjie general variational matrix M corresponding to the system ( 5 . 2 ) with respect 
^1,^1’ ^2) is given by 


to 


[Yv‘ 


+ + + -AF 2 P,Y2 pi{Ni-Y,-Z,)\ 


i 




-Oil 

^2{N2 - Y2) 


-di 0 

0 -di 

HN2 - 5^2) 0 


0 

0 

mu 


J 


ffiii 


~[P2Yi + X2Z1 + {0:2 + ^2 -t- (1 — a')j^N2}]. 


At 


the 


Mr 


^4 


ftie 


equilibrium point Pi the variational matrix Mi is given by 

f-{i'i+o:i + di+Si) 0 0 

5 i —di 0 0 

-ai 0 -di 0 

\ I32N2 X2N2 0 -|a 2 + d 2 + (l-a')^^ 2 }y 

^cteristic polynomial corresponding to above matrix is given by 

(di + + aiiP -f a2ip + 03 } = 0 , 




ai 


a 2 


fl 3 


i/i -f Qi -P 2di -1- di + |a2 + d2 + (1 — ) 

(i^i - 1 - Q:i -f di ■+■ di) di + |a 2 -I- d 2 -I- (1 - a')^N2 
L t K2 

+di |a 2 -t- d 2 -+■ (1 - (Y)~N2^ - ^2^2^—, 

(f/i + ai -h di + di)di {a2 + d2 + (1 - aO^iVs} - (di/52 + Aad^Aa^, 
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and 

0-l0>2 ~ 03 — 


1/1 + ai + 2di + di + |o;2 + <^2 + (1 ~ + ai + di + ( 5 i) 


+ ') CK2 + <^2 + (1 — a') r^N2 

K2 


} 


+[di + |oi2 + ^2 + (1 — Ci')^N2^]di |q;2 + ^2 + (1 — 

^ { 


-I 3 iP 2 ^N 2 

di 


+ 0^1 + di + ( 5 i + ] Q^2 + ^^2 + (1 “ 0!) A^2 

j^2 


-\-^i^N 25 iX 

di 


We note that aia2 - 03 > 0 if < 1 and so Pi is stable. For > 1 , obviously Pi is 
unstable and P2 exists. 

At the equilibrium point P2, the variational matrix M2 is given by 

l-{PiY2Pyi^ai+di + 5i) -AF 2 ^i{Ni-Yi- Zi)\ 

M2 = 


V 


— di 0 0 

— CKi 0 — di 0 

HN2-Y2) \2{N2-Y2) 0 mu 


J 


where mu — ~[P2Yi + A22^i + {0:2 + d2 + (1 — o!)^^N2}\- 
The characteristic polynomial corresponding to above matrix is 

(di + + giip"^ + g2'il) + 53} = 0 , 


where 

5 i = 
92 = 

9 z = 


di + ^il^ + z^i + Oil + di + di + ^ 2 X 1 + ^2.^1 + s 0:2 + d2 + (1 — o!) -j^N2 ? , 

j\2 


{■ 




{P1Y2 + i^i + CKi + 2di + di) 


P2Y1 + A2.^i + + d2 + (1 — 

+di{l3iY2 4- + Q!i + di + di) — /5il^(di + Oil) — Pi/d2{N2 ~ Y2){Ni — Yi — Zi), 

+ Q:i + di + di) ^2Yi + Aa-^i + |q: 2 + d2 + (1 — a')-j^N2'^ 


+PiY2{5i + Oil) 


P2Y1 + A2.^i + |o;2 + d2 + (1 — o!)-^N2 


r2 


-l3i{Ni - Yi - Zi)(iV2 - y2)(5iA2 + diP2)- 

It can be checked that 9192-93 > 0 . Hence by the Routh-Hurwitz criteria the equilibrium 
point P2 is locally asymptotically stable. 


Nonlinear Analysis and Simulation: It is noted that system (5.2) is bounded by the 
system ( 5 . 2 ) with Oi = 0 , 02 = 0 , which on using comparison theorems (Lakshmikantham 
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and Leela 1969) implies that the solution of (5.2) is bounded by the solution of (5.2) with 
ai = 0, for i = 1, 2. Using a heuristic method by setting Zi = the system (5.2) with 
ai = 0, aa = 0 can be shown to be globally stable (see Appendix I). Hence we speculate 
that the nontrivial equilibrium point of the system (5.2) may be globally stable if we start 
away from Pi. To further illustrate this, and to see the effects of various parameters on 
the equilibrium level of the infective population, the system (5.2) is integrated using 
the fourth order Runge-Kutta Method by using following set of parameters (Greenhalgh 
1990); 

Pi = 0.00000022 = p 2 , vi = 0.012, ^i = 0.0005, di = 0.0004, 

5i = 0.00002, A = 10, Aa = 0.00000011, aa = 0.045, ra = 1, 

da = 0.04, a' = 0.999, S 2 = 0.0000002, Qo = 20, 5o = 0.001, K 2 = 1000000. 

The equilibrium values of Yi, Zi, Ni and Y 2 corresponding to Pa are as follows, 

= 701.53, Zi= 35.077, iVi = 24123.07, 72 = 1761.68. 

Simulation of (5.2) is performed for different initial positions 1, 2, 3 and 4 as shown 
in Fig. 5.1. In Fig. 5.1, the infective population is plotted against the susceptible 
population, which indicates that it is plausible for the system to be globally stable about 
this equilibrium point Pa provided we start away from other equilibrium point. Further to 
see the effects of environmental discharges and immigration on the equilibrium level of the 
infective population, in Figs. 5.2 and 5.3, the infective population is plotted against time 
for different Qa, Si and A. From these figures, we conclude that the endemic equilibrium 
level of the infective population increases as any of A, iVa and 61 , increases. 
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Susceptible popualtion 


gure 5.1; 


Variation of infective human population with susceptible human population. 



Figure 5.2: Variation of infective human population with time for different rates of 
immigration and different rates of cumulative environmental discharges. 
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Figure 5.3: Variation of infective human population with time for different rate coeffi- 
cients corresponding to movement of human population from infective class to reservoir 
class. 

5.2.2 Case II: Q=.Qq + INi 


In this case it is sufficient to consider the following subsystem of the system (5.1). 


Fi = p,{N,-Y,-Z^)Y2-{u, + a, + d^ + 6,)Yu 
Zi = 5iYi-diZu 

Ni = A-diNi-aiYu (5.4) 

^2 = 02{N2-Y2)Y, + X2{N2-Y2)Zi-i^a2+d2 + {l-a')^^N2^Y2, 

N2 = r 2 N 2 (l - - 0^2N2 + S 2 N 2 E, 

E = Q^ + in^-SqE. 


In this case the region of attraction is 

T( = [{Yi,ZuN,,Y2,N2,E):0<Y, + Z^<N^<y, 


0 < ^2 < iV2 < ^ f r2 - a2 -1-<52^ 
’’2 \ 5o 


, 0<E< 


Qij) 

5o- 


and the model is well posed. The result of equilibrium analysis of the system (5.4) is 
stated in the following theorem. 
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them and form a separate class of reservoir population which helps in spreading malaria 
without getting affected themselves (Nchinda 1998). 

Although, as described earlier, usual mathematical models for the spread of malaria by 
considering different aspects have been proposed and analyzed (Chapter 4), none of the 
models takes into account the role of the reservoir population as well as environmental 
effects. 

In this chapter, therefore, a mathematical model for the spread of malaria with a reser- 
voir class is proposed and analyzed. The effect of household discharges conducive to 
the growth of the mosquito population on the spread of malaria is also considered in 
the model. As in the previous chapters, here also the following two types of population 
demographics are considered, 

(i) a population with constant immigration, and 

(ii) a population with logistic growth. 


5.2 Malaria Model with Constant Immigration 


We consider here an SIS model with immigration, where the human population density 
Ni{t) is divided into three classes namely, a susceptible class Xi{t), an infective class 
Yi{t) and a reservoir class Zi{t). The mosquito population of density N 2 {t) is divided 
into a susceptible class Xiit) and an infective class Y 2 {t). It is assumed that the rate 
of growth of density of mosquito follows a generalized logistic model (Gao and Hethcote 
1992). In view of the above and by considering the criss-cross interaction of the mosquito 
population with the human population, a model is proposed as follows, 

AA = A-d,Xi-PiXiY2 + u^Y,, 

Yi = P1X1Y2 — {di -\- v\ + a\ -\r ^i)Ti, 

Zx = - dxZx, (5-1) 

Nx = A-dxNx-axYx, 
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> 0 . 


s 

— di |q; 2 + da + (1 — (i^i + ai + di + di) ~ ()^2 + •^2^) 

Hence by the Routh-Hurwitz criteria this equilibrium point is stable under the condition 
R; < 1 and if R; > 1 it is unstable and the third equilibrium point exists. 

The variational matrix M3 at equilibrium point .P3 is given by 

( mn P1Y2 Pi{Ni-Yi-A) 0 0 


M 3 = 


( 5 i -di 0 

— CKi 0 — di 

/ 32 iN 2 -Y 2 ) X 2 {N 2 -Y 2 ) 0 

0 0 0 

0 0 I 


0 

0 

77144 

0 

0 


0 0 

0 0 

77145 0^ 

0 -<5o ) 


where 

TTlii = — (ySlia + ^ 1 + CKi + di + di), 77144 = — [;52^1 + + {<^2 + ^2 + (l — a')^JV2|j, 

= {^2Yi + Aa-^i) - (1 - o/)j^Y2 and 77755 = -{ra - 0-2 + da^}. 

From the characteristic polynomial corresponding to matrix M3, we see that one root of 
the above matrix is — di and other roots are given by the following polynomial, 


where 

a 4 = 

03 = 


7/;® + + a2'ip^ + aixjj + oo = 0 , 

Sq r 2 — oi 2 52E + P1Y2 + + 0^1 + 2 di + + P2Y1 H- X2E1 

+ joa -f da + (1 - Y)-^N2^ , 

^0 [^2 ■“ 0^2 + ^2^ + P 1 Y 2 + + oti + 2di + (5i + j32Yi + X 2 Z 1 

+{o'2 + ^2 + (1 - a')~^iV2} + (r2 - 0^2 + S 2 E) P 1 Y 2 + z/1 + ai + 2di 
^2 J •- 

+<^1 + /32Y1 + X2Z1 + |o'2 + ^2 + (1 — 

■i-{^iY2 + -f- + <Ji) 


a2 = 


P2Y1 + X2Z1 + 1^2 +d2 + {1 - + di 

-^W2Yi + X2Z1 + ^2)^1 +PiY2{Si + ai) + ^i/ 32 {Ni -Yi- Zi){N2 - ^2), 

~ 0:2 + S2E) ^fiiY2 + z/1 + ai + 2di + ^1 + P2Y1 + X2Z1+ 


|a2 + da + (1 — Y)-^N2 ^ + di 


K2 


+ ( 7-2 - 0!2 + hE){0iY2 + vi -f ai + di + (5i) 


02Y1 + Aa^i + {aa + da + (1 — a')-^N2} + di 


+W1Y2 + vi + Q;i + di + ( 5 i) 
+ 


02Y1 + X2Z1 + /aa + da + (1 — a!)~N2 

t iva 


hYi + Aa^i + jaa + da + (1 - aO^iVaj + di 


(di + (Jo) 


I di (do 


+ 7-2 - ora + hE) 
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+d\ 5 Q{P\Y 2 + i/i + Qi + + 5 i) 

+/ 3 i^ 2 (< 5 i + CKl) 


< 5 o + r2 - 0:2 + 62E + P2Y1 + X2Z1 + <j 0:2 + £i2 + (1 - o!)^N2 

-^2 


+/5i(iVi — Yi - Zi){N2 - Y 2 ) 1^1 A 2 + ^ 2(^0 + di + r2 - 0:2 + (52.®)| , 

{‘ 


ai = <5 o(’"2 - ct2 + 52E){piY2 + ui + ai + di+ Si) |;02Yi + A 2 Z 1 + -j 02 + ^2 + (1 - a0^'^2| 


+(r2 - a2 + 52^5) (^1^2 + Z.1 + ai + + Si) |^2Yi + A 2 Z 1 + |a2 + ^2 + (1 - } \ di 


r 2 




+(/ 3 i^ + z/i + O'! + di + ( 5 i) 


I^ 2 Yi + >^2^1 + |^QJ2 + ^2 + (1 — o!)-^N 2 


r 2 


diSo 


+ 


^ 2^1 + ^2^1 + |q' 2 + ^2 + (1 — a')^iV2 


<^0(7^2 0^2 + S2E) 


+(il(^o(^2 — Oi2 62 E){PiY 2 + Ui + (Xi + di + 61) 


+A^ 2 (< 5'1 + Oil) 


{ 


^0 + 02^1 + A2^1 + ^2 + ^2 H" (1 0!) ~^^2 (^2 — 0^2 + 

K2 

^2 


+<^o{/^2^1 + ^ 2^1 4 - Q!2 + <^2 H" (1 

-^2 


+/ 3 i(i'^i — Zi){N 2 ^2) [<^iA2(<^o + r2 ~ 0^2 + S2E) 

+P2 {Sodi + {di + ^o)(^2 -0^2 + hE)J^ , 

<^o = ^0(^2 — a2 + S2E){piY2 + i^i + c^i + di + 61) 

P2Y1 + A2Z1 + |a2 + d2 + (1 - aO;^^2| 


/?2n + A2Y1 + \a 2 + d 2 + {l- a')^N 2 

K2 


di 


+^1^2 (<^1 + 0 'i)(Jo (^2 “ a2 4 “ 52E) 

+^l(iVi — Yi — ^i)(iV2 — 3^2) (^2 “ Ck:2 + <^2 -E)^o(<^1 A 2 + J02<^l) 


+^2-^2 


hYi + A2Y1 - (1 - o.')^Y 2 


PiiNi -Yi- Zi)ail. 


Using the Routh-Hurwitz criteria, we see that this equilibrium point P3 is locally asymp- 
totically stable if the following conditions are satisfied, 


04 > 0, 


(14 CL 2 
1 03 


> 0 , 


04 02 Oo 
1 03 

0 O4 02 


> 0 , 


04 02 oo 0 

1 03 oi 0 

0 04 02 Oo 

0 1 03 oi 


> 0 , 


04 02 oq 0 0 

1 03 oi 0 0 

0 O4 O2 Oq 0 

0 1 03 Oi 0 

0 0 04 02 Oo 


> 0 . 


The first two conditions are obviously true. If the next two are satisfied then so is 
the fifth. Hence the equilibrium .P3 is locally asymptotically stable under the conditions 
mentioned in the theorem. 


Nonlinear Analysis and Simulation: We note that the system ( 5 . 4 ) is bounded by its 
corresponding system when ai = 0 in T[. As in the earlier case using heuristic approach 
the system (5.4) with ai = 0 can be shown to be globally stable about the non-trivial 
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equilibrium point provided we do not start from the other equilibria (see Appendix II). 
Hence we speculate that P3 may be globally stable provided that we start away from 
other equilibria. 

To show this, the system (5.4) is integrated using the fourth order Runge-Kutta method 
using previous values of parameters with Qo = Qa ^.nd an additional parameter I = 
0.0005, which satisfy the local stability conditions mentioned in Theorem 5.4. 

The equilibrium values of li, Zi, Ni, Y 2 ,N 2 and E are as follows, 

Yi = 4181.8, Zi = 209.2, Ni = 19771.8, Y 2 = 14051.1, N 2 = 960977.1, E = 29885.6. 

Simulation of the system (5.4) has been performed for different initial positions 1, 2, 3 
and 4 as shown in Fig. 5.4. In this figure, the infective population is plotted against 
the susceptible population and from the solution curves, we conclude that the system is 
globally stable about this equilibrium point P3 under the set of parameters considered. 
Also from Figs. 5. 5-5. 9, where the infective population is plotted with time for different 
^2) <^2) Qo, I, and (5i, we note that the infective population increases as any of these 
parameters increases, which correspond to the growth of the mosquito population and 
environmental discharges. 
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gure 5.4: Variation of infective human population with susceptible human population. 



Figure 5.5: Variation of infective human population with time for different intrinsic 
growth rates of mosquito population. 
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Figure 5.6: Variation of infective human population with time for different growth rate 
coefficients of mosquito population due to environmental discharges. 



Figure 5.7: Variation of infective human population with time for different rates of 
cumulative environmental discharges. 
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Figure 5.9: Variation of infective human population with time for diflFerent rate coeffi- 
cients corresponding to movement of the humam population from the infective class to 
the reservoir class. 
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5.3 Malaria Model with Logistic Population Growth 

Now let us consider an SIS model where the human population growth is logistic so 
that both the birth and death rates are density dependent in such a manner that the 
birth rate decreases and death rate increases as the population size increases towards its 
carrying capacity. Keeping in view the other considerations discussed in previous section, 
a mathematical model is proposed as follows: 

= \hi-a^N^N,-\k^{l-a)^N^Xx-piXiY-i + uiYu 
L Ai J L J 

Ki = /5 iXiF2- L +ai + di + (5i + (l-a)^lyi, 

L ill J 

iVi = ri(l - - aiKi, (5 9) 

^2 = (^>2-a'^iV2)iV2-{d2 + (l-a0^iV2}X2-;d2^2Ki -A 2 X 2 Z 1 -a2>2, 

Y^2 = P 2 X 2 Y 1 -h ^2^22^1 - {^2 + c!2 + (1 — a')-^N2}Y2, 

K 2 

No 

^2 = r2N2{l — ~^) — 012^2 ^2N2E , 

K 2 

E = Q{Ni)-5oE, 

Ni = Xi+Yi + Zi, N2 = X2 + Y2,r2 > a2, 5i < di, 0 < a < 1, 0 < a' < 1, 

Xi (0) = Xio > 0, Yi (0) = 7io > 0, ^i(O) = Zio > 0 , X2(0} = X 20 > 0, ^20 > 0, E{0) = Eq > 0. 

Here bi and di are the natural birth and death rates, Ti = bi — di is the growth rate 
constant; ri and Ki are the logistic growth rate and carrying capacity of the environment 
corresponding to human population; r 2 and K 2 are the logistic growth rate and carrying 
capacity of the mosquito population in the environment. All other parameters are defined 
in the previous section. For 0 < a < l, the birth rate decreases and the death rate 
increases as N, increases to its carrying capacity K^. When a = 1, the model could be 
called simply a logistic birth model as all of the restricted growth is due to a decreasing 
birth rate and the death rate is constant. Similarly, when o = 0, it could be called a 
logistic death model as all of the restricted growth is due to an increasing death rate and 
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the birth rate is constant. 

As in the previous section, here also the following two cases are considered, 

(i) the rate of cumulative environmental discharges Q is a constant, and 

(ii) Q is human population density dependent. 


5.3.1 Case I: Q = a Constant 


In this case since X\ + yi + Zi = N\ and ^2 + 12 = Z^-, the behaviour of the system (5.9) 
can be studied by the following system. 


5^1 = ^i[N,-Y,-Z,]Y2 
Zi 
Ni 
Y2 


+ Q!i + di + (5i + (1 - a) 


niVi 


Ki J 


5iY,-i^dr + {l-a)^N,^Zu 

Ni 

nN,{l - - a,Y^, 

MNi - n)i'i + A2(iV2 - n)Zi - (02 + * + (1 - o') 


Yi, 




(5.10) 


where as in the previous section, 


A'2 

N 2 = lim sup N 2 = — {r2 — 0:2 4- SqE} and E = lim sup E = 

t^OO T2 t-^OQ 


<5o' 


The results of an equilibrium analysis of the system (5.10) are stated in the following 
theorem. 


THEOREM 5.5 There exist the following three equilibria corresponding to system (5.10) 
namely, (i) ^1(0, 0, 0, 0), (ii) ^2(0,0, A'i,0) and (Hi) E 3 {Y(,Zl,N(,Yf), which exists if 


Pi {^2-^ ^2 


61 


di + (1 - a)ri 


N2Kr 


{ 


a2 + d2 + {I — {*^1 + 0^1 + di + + (1 — a)^i} 


= Rs > 1. 


This equilibrium is unique if > 0 for all Ni > 0, and when N(q > y, where Yi is 
given in equation (5.15) of the proof of this theorem. 
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under condition (5.16a), where D* = {di -I- (1 - 

Now plotting (5.15) and (5.11) in the N-Y plane (Fig. 5.10), we get an intersecting point 
(N^,V*) and then using (5.12) and (5.13), and V^* can be determined respectively. 
If the slope ^ of (5.15) is positive, then there exists a unique (N*, Y*) when > £ 
and corresponding to this, a unique positive equilibrium point exists. 


5. 3. 1,1 Stability Analysis 

In the following we present the linear stability results of these equilibria. 


THEOREM 5.6 The equilibrium Ei is unstable, the equilibrium E2 is stable if Rs < 1 


and unstable if Rs > 1 and the equilibrium E3 is locally asymptotically stable Oq > 0 
' 03 Oi 0 I ’ 

> 0, where oq, Ui, 02, and as are given explicitly in 


^3 ^1 

1 a2 


> 0 , 


1 02 Co 

,0 as ai 

the proof of the theorem. 


Proof: The variational matrix M corresponding to the system of equations (5.10) 
^t{Y^,Z,,NuY2) is given by 



/ mu 


-Oi 

\P2{N2-Y2) 


-P1Y2 

~{di + {l-a)j^N,} 
0 

^ 2(^2 — Y2) 


AK2-(l-a)^F, 

0 


Pi{Nr-Y,-Z,)\ 

0 

0 

m.i.i J 


where mu = -{A Fa + + a;i -I- -f dj + (1 - a)^yV, } 

and m44 = — [/JaFi + X2Z1 + 0:2 + ^2 + (1 — of)j^N2]. 

It IS easy to see that the variational matrices corresponding to £'1 (0,0, 0,0) has one 
positive eigenvalue implying instability of this equilibrium. 

The variational matrix corresponding to E, has characteristic equation 


ip^ + + A2‘ip -f ids = 0, 

where 


i^i + ai+di-hdi + (l-a)ri-t-di + (l-a)rl-Ha2-^-d2^-(l-a')— 

/f■; 


No 
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A2 — 


^3 


{vi + ai + + (5i + (1 — o)ri} {di + (1 — a)ri} 

+{di + (1 - 0)7-1} {^2 + ^2 + (1 - aO^iVaj - P1P2K1N2, and 
{ui + Qi + di + (5i + (1 - a)ri} {di + (1 - a)ri} |q :2 + ^2 + (1 - 0 !) 
—P1K151X2N2 ~ Pi^ 2 KiN 2 {di + (1 — a)ri}. 


ll 

Ko 


No 


From above equation it is easy to see that E2 is stable if < 1 and unstable if > 1. 
The variational matrix M3 at equilibrium point £^3(17, N^.Y^) is given by 


/ 


Ms 


m 


11 


m. 


'22 


<^1 

—ay 0 

KI 32 {N 2 -Y^*) X 2 iN 2 -Yi) 


-PiYi l 3 iY^*-{l-a)j^Y,* A(iVr-yi*-^r)\ 


"(1 - 


m 


■33 


0 

0 


m 


'44 


J 


where m*i = - {P1Y2 + Ui + ai + di + 5 i + {1 - - {di + (1 - CL)j^N*Y 


m 


33 


.rL 

Ki 


(2N* — Ki)^ and 77144 — + ^22^1 + |cr2 + (i2 + (1 ~ • 


The characteristic polynomial corresponding to above is given by 


+ Clo — • 0, 


where 


az = -{Ki + ^22 + ^33 + ^44} > 0^ 


a2 




^11 (^22 + ^33 + K4) + ^22 (^33 + ^44) + +^33 ^4 

+ijAn- + “1 {ftn- - (1 - - z;)A (^2 - n-). 


- 771^2 


+ ^22 ^33 77144 


+ ^33 "^44 "^11 


+ 77744 777 i 4 77122] 


+PiY;at{l - a)^Zl - Ks + m-„)5iPiY; 

-a, (ftvy - (1 - [">i + n>«l - - Y; - Zl)(N^ - ij) 

+(my + my)P^(Nl - Y; - Zl)h(N 2 - Y^), 

7*1 

Oo = 77X41 ^22 ”^33 ”^44 ” ^^44 /^1^2**^l(l “ 

+5 i/ 3 iF 2* rn\^ - " Zl)P 2 {N 2 - i?) ^*22 


+ {av; - (1 - a. mi mi + iiAAj m^W - K,* - Zr)(JV2 - 1?) 

- y,- - Zi*)(i - o)^zra:,A2(JV2 - y;). 

ill 
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By the Routh-Hurwitz criteria, conditions for the local stability of the system are the 
following: 


as > 0, 


Os ai 
1 Cl 2 


> 0 , 


as Cl 0 

1 0,2 ao 

0 Os ai 


> 0 , 


Os ai 0 0 

1 a2 ao 0 
0 Os Oi 0 
0 1 02 Oo 


Remark: It is noted that as > 0 and the second inequality is obvious if JV* > ^. So 
if the third inequality is satisfied then the fourth one is obvious provided Oq > 0. Thus, 
under these conditions this equilibrium point is locally asymptotically stable. 


Nonlinear Analysis and Simulation: As before, here also we speculate that the 
equilibrium E 3 is globally stable under local stability conditions. 

To illustrate this, the system (5.10) is integrated by the fourth order Runge-Kutta method 
using the following parameter values in the simulation. 

Pi = 0.00000022 = P 2 , = 0.012, oi = 0.0005, di = 0.0004, 5i = 0.00002, 


a = 0.3, ri = 0.0003, A, = 50000, A 2 = 0.00000011, ^2 = 0.045, a' = 0.999, 

K 2 = 1000000, r2 = 1.0, d2 = 0.04, 02 = 0.045, §2 = 0.0000002, = 20, 5o = 0.001. 

The equilibrium values of Fi, Zi, Ni and F 2 are as follows 

Fi = 6986.664, = 262.425, Ai = 31540.491, F 2 = 17064.251. 

Simulation is performed for different initial positions 1, 2, 3 and 4 shown in Fig. 5.11. In 
this figure, the infected human population is plotted against the susceptible population 
and also against total population N. From the solution curves, we conclude that the 
system is globally stable about this equilibrium point under above conditions provided 
we start away from other equilibria. Also in Figures 5.12 and 5.13, the infected population 
is plotted against time for different ri, Q^, Si and it is noted that with the increase of any 
of these parameters, the infective population increases showing that the disease spreads 
faster due to these effects. 
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Figure 5.12: Variation of infective human population with time for different intrinsic 
growth rates of human population and different rates of cumulative environmental dis- 
charges. 
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Figure 5.13: Variation of infective population with time for different rate coefficients 
corresponding to movement of human population from infective class to reservoir class. 

5.3.2 Case II: Q = Qq-\-INi 


As indicated before, in this case it is sufficient to consider the following equivalent system 
of the system of equations (5.9). 

^'i = /^i(M — Fj — + 0^1 + di -(- -f (1 — K], 

Zi = (5.17) 

Ni = ri (l - iVi - aiFi, 

F 2 = /?2(A2 — F2)yi 4- A2(A’2 — ^)2'i — |q;2 + ^2 + (1 - F 2 , 

1 iv2 J 

N2 = r 2 iV 2 Oi 2 N 2 + 52N2E, 

E = Qo + lNi-SoE. 

The result of an equilibrium analysis is stated in the following theorem. 


THEOREM 5 . 7 There exist the following five equilibria, Tianiely 

(i) E, (0, 0, 0, 0, 0, ^), (H) (0, 0, 0. f {r, - aj + 


(Hi) E, (0, 0, if,, 0, 0, , (iv) E, (0, 0, if,. 0, N^, 2^), 
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Ko { 5i 1 

where N 2 = — < r 2 - 0:2 + -riQo + ^-^ 1 ) 

’'2 t <^o J 


and (v) EsCVi, Zi, Ni,Y 2 i N 2 , E) . E 5 exists if 


a| 

^ R„ A. 

h 1 


1 P2 + ^23^ 
[ di 

+ (1 - a)ri J 

+ <^2 + (1 


{ui + ai + di + 5 i + (1 - a)ri} 


This equilibrium is unique if > 0 and Nio* > where Yi and N*q are given in the 
proof of the theorem. 


Proof: The existence of the first four equilibria is obvious. The existence of E 5 is shown 
by the isocline method. Setting the right hand sides of (5.17) to zero, we have the 
following equations, when Ni ^ 0 , Ni ^ Ki, N 2 ^ 0 , 




(5.18) 


E 


Qq + INi 
^0 


^1 = 




K2r . . ^2. 


, ,f, iV2 = — {r2-Q!2 + ^(Qo + ^A^l)}, 

di + (1 — ojj^Ni . T 2 ^0 



Modelling the Spread of Malaria with Human Reservoir: Environ- mental and Demographic 
148 Effects 


¥> = 


f p \ <^1 1 

[yiiV2 


a2 + d2 + (l-a0;^^2 + ] 

^\P2 + ^ 

, 1 




¥> = 


+ ai -f + (^1 + (1 — a)^Nij Yi 
Pi [iVi - {: 




(5.20) 


(5.21) 


From (5.20) and (5.21), we get 




hPiN^Ni - \ 

1 ^ 0:2 + ^2 + (1 — a0^-^2| 


\vi + Q;i + di + + (1 — 

1 ill J 

PI 

ft] 

fl+ ft 1 

[ di + (1 - a)^N\ J 

1 iV2 + j 

+ cti + di + di + (1 — 



Ko ( 5o 1 

where N 2 = —Ir 2 -a 2 + j-iQo + 


(5.22), 


Clearly in the N-Y plane (5.18) is a parabola with vertex at and passing 

through (0, 0) and (i^, 0). 

From (5.22), the following points are observed. 

(i) For Ni = 0, 


Y 


1 — 


(ui + ai+di + di) 1 

0(2 + ^2 + (1 - a') 1 

(r, - a, + J,| 5 ) 

} 


I 

^Pi 1 

yi 

\ K2 f r Qo\ 1 

1 y2 ~ 0:2 + J “b "b ^1 d" dj + > 


< 0 . 


(ii) For Ni = Ki, 


P1P2 ^2^1 {^2 + (^2 + (1 •“ ^')~^N2}{^i + Oil + dj + - 4 - (1 — a)-^Ki} 

U = £!_>0. 

P2 [A {1 + ^ (1 + {^i + Q:i + di + di + (1 - a)ri}] 

provided 

PiP;*n;ki = 

{0:2 + <^2 + (1 — a )— iV2}{i/i + O'! + di + + (1 — a)ri} 

where ft" = ft -f- A, 5,:^. 

(iii) For Yi = 0, we get the following cubic in Ni 


hiNf + h2Nl + hzNi + ^4 = 0 , 


(5.24) 
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where 


hi 


h 2 


hz 




A {P 2 di + A2(5i)— I + — ^7-2 - 0:2 + P 2{1 - a) 

(< 


^0 7-2 

+ ^ 0:2 + 0^2 + (1 “ o') ^r2 — Q!2 + <^2 
K 2 


n_ 

Ki 


(1 - Q)n 
Ki 


di{l — a')S 2 l 


So 


Qo ^ T*! 


PlW 2 di + A2 ^i)- 


7-2 


7-2 - Q;2 + S 2 


Qo\_ 

SoJ 


(1 ~ + ai + di + 5i)(l — > 


di 


(1 ~ + cti + di + (5i) 


+ (1 — 1*^2 + <^2 + (1 — 0,'){r2 — a2 + 52 ^)| 

-(1 " |o^2 + d2 + (1 - a') (^2 - 0:2 + *^2^^ I (t^i + cti + di + (5i), 

- 0-2 + ^2^) I ( 1^1 +ai + di +(5i). 


-di |q ;2 + ^2 + (1 - o')(r2 


It is noted that under condition (5.23), the coefficient of Nf i.e. hi is positive and the 
constant term /14 is negative in (5.24), thus implying that it has one positive root. If 
the slope is positive for A^i > 0, then under the above condition plotting (5.18) 
and (5.22) in the N-Y plane (Fig. 5.14), we get a unique intersecting point as (yi,iVi) 
provided Y 2 , N 2 -, Zi and E are then easily derived by using (5.19) and (5.20). 


5.3.2. 1 Stability Analysis 


The linear stability results are stated in the following theorem. 


THEOREM 5.8 The equilibria Ei, E 2 andEz are unstable, the equilibrium E 4 , is locally 
asymptotically stable if Rg < if Rs > 1 ^ is unstable and the equilibrium E 5 exists, 
which is locally asymptotically stable provided 


oo > 0 , 


05 03 

1 04 


> 0 , 


as as ai 

1 a4 a2 
0 as as 


> 0 , 


as as ai 0 

1 (X 4 0,2 ^0 

0 ds ^3 

0 1 a4 a 2 


>0, 


as as ai 0 0 

1 a4 a2 ao 0 

0 as as ai 0 

0 1 a4 as ag 

0 0 as as ai 


> 0 , 


where ag, ai, as, as, a4 and as are given explicitly in the proof of the theorem. 
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x(?7i55 - So) -f mss rriu rrisol-So + mn) — rti^o So (rnn 4- 17122) — 5o 77133 77155 

^ y»l ^ 

X (77122 + 77144) — Sonin rrioo (77144 -t- 77155) -I- (^o — '^u ~ ^55) A ^20^1(1 — a)-~Zi 

Ki 

+Si^iY2{—So {nrioz -I- 77144 + 77155) + rrioziniu + 77155) + ntiunfiob] 

-|-Q!i 7?1 i 3{— 5o(77l22 + 77I44 4" 77I55) 4- 77122(77144 4" 77I55) 4“ 77I4477I55} 

-rnu[N 2 - ¥2)^2 (77122WI33 4- 7712277155 4- 7713377155 - <5o(7?i22 4- 77133 4- 77155)) 

-rnu{N2 - y2)A2a;i(l - a)^Zi - rnuiN2 - 5^2)AA2(A - "^33 - "^55), 

Til 

ai = Attiii 77122 (77133 77144 4- 77144 nfl55 4- 77155 77133) 

4-77133 77144 77I55 {(5o(?71il 4- 77122) - TTln 77122} 

4-Q:i(l - a)—Zij 3 iY 2 {7714477155 - Soiniu 4- 77155)} 

Til 

4 -^iA^{A ^44 nio 5 ~~ nriii viob nrisz 4- 77155 rriozSo 4- 77133 A ^44} 

4-o;i 77113 ?7i44 - rn 22 niu nios + A ^55 {17122 + ni^i)} 

+77114(^/2 - >2) A {’^27713377155 - A(’A2’7l33 + mzzrrioo 4- 77I2277I55)} 

-77114 (iV2 - Y 2 )X 20 'lil - a)^Zi{So - 77I55) - AiV277l45771i4toi 

Al 

-<5 iA 2771 i 4(7V2 - ^2){77l33"^55 - <5o(’n33 + 77I55)}, 
ao = SQrn/^4 rrioo (^—mn7ri22nizz + (XiP{Y 2{1 — a)-j^Zi — airn22i7izz ~ SiftiY2Tno^ 
4-771i477l55<5o(-/V2 — Y2) |A’A2’’^33 4" A2Q;i(l — ~ ^lA277l33|' 

4- A-^277l45771i4 77l22/Q;i . 

By the Routh-Hurwitz criteria, if the following conditions are satisfied then this equilib- 
rium point is stable, 


U5 ^ 0, 


0.5 az 
1 04 


> 0 , 


05 03 oi 0 0 

1 O 4 02 Oq 0 

0 05 03 oi 0 

0 1 04 02 Oq 

0 0 05 03 oi 


0.5 az fli 
1 04 02 

0 05 03 


> 0 , 


05 

0^3 

0i 

0 

1 

(X4. 

02 

Oq 

0 

as 

03 

Ol 

0 

1 

04 

O2 


> 0 , 


05 03 Ol 0 0 0 

1 04 02 Oq 0 0 

0 05 03 oi 0 0 

0 1 04 02 Oq 0 

0 0 05 - 03 oi 0 

0 0 1 04 02 Oq 
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Remark: It is noted that > 0 and the second inequality is obvious if Ni > If 
next three inequalities are satisfied then so is the sixth provided ao > 0. So for numerical 
simulation, we choose parameters such that Ni > 

Nonlinear Analysis and Simulation: Here also the global stability of is speculated 
under local stability conditions. To show this, the system (5.17) is integrated by the 
fourth order Runge-Kutta method using the same parameter values as in the previous 
subsection with Qo = Qa and an additional parameter I = 0.00005. The equilibrium 
values of Yi, Yi, Ni, Y 2 , and E are as follows: 

Y = 7014.7, = 263.9, iVi = 31359.4, 72 = 17281.5, iVs = 962135.9, E = 35679.7. 

Simulation is performed for different initial positions 1, 2, 3 and 4 in the interior of 
the region of attraction as shown in Fig. 5.15. In this figure, infected population is 
plotted against susceptible population. From the solution curves, we conclude that the 
system is globally stable about this equilibrium provided that we start away from other 
equilibria. Also from Figs. 5.16-5.20, where the infective population is plotted with time 
for different r^, <52, Qo, I and (5i, we note that the infective population increases as any 
of these parameters increases showing the effect of these parameters on the spread of 
disease. 
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Figure 5.15; Variation of infective human population with susceptible human population. 



Figure 5.16: Variation of infective human population with time for different intrinsic 
growth rates of mosquito population. 
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Figure 5.17: Variation of infective human population with time for different growth rate 
coefficients of mosquito population due to environmental discharges. 



Figure 5.18: Variation of infective human population with time for different rates of 
cumulative environmental discharges. 
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Figure 5.19: Variation of infective population with time for different 1 . 



Figure 5.20: Variation of infective population with time for different rate coefficients 
corresponding to movement of human population from infective class to reservoir class. 






In this chapter an SIS model for malaria with human population as reservoir is pro- 
posed and analyzed by considering environmental and demographic effects. The cases of 
constant rate of cumulative environmental discharges as well as a population dependent 
rate of cumulative environmental discharges are considered. The following two types of 
demographics for human population are considered, (i) constant immigration and (ii) 
logistic population growth. The threshold condition for spread of malaria is derived in 
each case. It is shown in each case that if the threshold is greater than one, the nontrivial 
equilibrium is always feasible and is locally asymptotically stable to small perturbations 
under certain conditions. If a public health measure is such that it keeps threshold less 
than one then disease will die out, i.e. the disease free equilibrium will be stable. By 
simulation it is observed that under local stability conditions the nontrivial equilibria 
seem to be globally stable, provided that we start away from the other equilibria. It 
is shown that due to environmental discharges, the mosquito population can grows very 
large leading to increased spread of malaria. Also if the rate of immigration or the growth 
rate of human population increases, the infective population increases implying further 
increase in the spread of malaria and the disease becomes more endemic. 



Chapter 6 


Modelling the Spread of a Carrier- 
Dependent Infectious Disease in 
Two Neighboring Habitats with 
Migration in Between 


6.1 Introduction 


As pointed out in Chapter 2, various types of carriers such as flies, ticks, mites, snails and 
so on, are responsible for the spread of infectious diseases in human populations who live 
in the region where densities of these carriers increase due to household discharges into the 
environment (Harold 1960, Harry and John 1962, Harry and Kent 1961). Some examples 
of such diseases are dysentery, gastroenteritis, diarrhea, cholera, measles, tuberculosis 
(Cairnoss and Feachem 1983, Taylor and Knowelden 1964). The migration of populations 
from environmentally degraded regions to cleaner regions also plays a very important 

role in the spread of infectious diseases as infected persons act as carriers or reservoirs of 
infection. 

Migration of the population is a common phenomenon in the same region, country or 
even outside the country due to considerations which may be economic, social, religious, 
political, environmental and so on. The migrating population carries with it all of its 
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Region with better sanitary Region with bad sanitary conditions 

conditions, where rich people live where poor people live 






^2 


Figure 6.1: Spread of disease due to immigration. 

traditional values and cultural heritage, including any diseases that are present in the 
population. After movement of the population into the new habitat, the susceptibles 
join the new susceptible class and the infectives join the new infective class and the usual 
spread of infection begins. 

The effect of migration of the population on the spread of infectious diseases through 
demographic changes has been studied in some cases (Bailey 1980). But the effect of 
migration of the population from an environmentally degraded region or habitat to an 
environmentally cleaner region has not been modeled and analyzed so far. 

In this chapter, therefore, we propose and analyze a model for the spread of an infectious 
disease between two socially structured populations (rich and poor). The rich population 
lives in a region with better sanitary conditions whilst the poor population lives in a region 
with bad sanitary conditions which is further affected by various household discharges 
conducive to the growth of a carrier population. The effect of a carrier on the spread of 
the disease in the region is also included in the proposed model as discussed in Chapter 2. 


6.2 Mathematical Model 

We consider that the population density N\ of the rich class is divided into two subclasses 
of susceptibles Xi and infectives Yi, who live in a region with better sanitary conditions 
of the habitat or city. The population of poor people with density iVj having susceptibles 
X 2 and infectives Y 2 live in a degraded region of the habitat with bad sanitary conditions 
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initial conditions for our analysis. 

Yy = ^iYy{Ni-Yy)-{vx + aiYdi)Yx + 9Y2, 

Nx = Ax — dxNx — 0!iyi + 9N2, 

Y 2 = {p2Y2 + X2CrAiN2-Y2)-{u2 + a2 + d2)Y2-9Y2, (6.3) 

N 2 — -^2 ~ d2N2 — 0 ! 2 ^ — 9 N 2 , 

where Cm = 7(5 - so + SxEm} and Em = ^- 

We note that the system (6.3) has only the nontrivial equilibrium, which is stated in the 
following theorem. 


THEOREM 6.1 There exists an unique nontrivial (endemic) equilibrium point, namely, 
Pi{Yx,NuY 2,N2), where 

Y _ Y \j Ex + 4(^2 + 9){dx + ax)Pidx{d 2 + 9 ) 9 y 2 

2 {d 2 + 9){di + ax)l3i 

Bx = [((^2 + 6){^iAx — di(i'i + Q!i + dx)} + 9A2^i — ^x9oi2Y2\, 

- B A xj B^ + 4/02(0^2 + 0^2 + 9)X2C„A2 
2^2 ( 0^2 A d 2 9) 

B = {P 2 A 2 — {6,2 + 9){X2Cm + 1^2 + 0:2 + ^2 + ^) — Q!2'^2C'm}j 

- {62 + 9){Ai - aiyi) -t- 9 {A 2 - Q! 2 ^) j {A 2 - o:2Y2) 

iVi — — — — and N2 = — T-, t:;: — . 

dx{d2 + d) (62 + 9) 


Proof: Setting the right hand side of the equations in (6.3) to zero and with some 
manipulation, we get, 


N = (^2 + 9){Ax - oiFi) + 9{A2 - 02 ^ 2 ) 
d\{d2 + 9) 


(6.4) 


JV,= 


(A 2 — 0^2^) 

((^2 + 0) 


(6.5) 


{d2+9){dx+ax)6xY^+[{d2+9){dx{i'x-\-ax+dx)-6xAx}-eA20x+piea2Y2]Yi -di{d2+e)eY2 = 0 

( 6 . 6 ) 
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and p2{a2+d2+6)Yi + [{d2+9){X2Cm+U2+a2+d2+d)+a2X2Cm-p2A2]Y2-X2CmA2 = 0, 

(6.7) 


which gives only one positive root Y 2 as follows, 


^ + 4^2(oi2 + (^2 + 0)X2CmJ^2 

2^2 (0:2 + ^2 + 0) 

where B = {p 2-'‘^2 ~ (da + 9){X2Cm + 1^2 + Oi 2 + d 2 + 6) — aaAaCm}. 

Using this value of Y 2 in equation (6.6), we get corresponding positive root Yi as follows, 

Y = -Si + yjB^ + AMi0{d2 + 9Y{di + ai)Y2 

2 Pi{d2 + 0){di + ai) 


where B\ = [(da + ~ di[i>i + ai + di)} + 9 A 2 P 1 — ,5i0a2^]. 

Then Ni and N 2 can be found by putting values of and Y 2 in equations (6.4) and (6.5). 
For Ni and N 2 to be positive, we should have conditions as Yi < ^ + (d2%) ) 

and Fa < ^5 which are satisfied in view of (6.6) and (6.7). Thus we have only one 
equilibrium point namely Pi(Yi, Ai, Fa, -A^a). 


Remark: We note from (6.7) that Fa increases as Cm increases. Also from (6.6), we 
note that (5^) is positive only when Fi < which implies that an increase in Fa 

causes an increase in Fi provided that the above condition is satisfied (for aa very very 
small it seems to be satisfied). 


6.3.1 Stability Analysis 

Now we present the stability analysis of this equilibrium. The local stability result of 
this equilibrium is stated in the following theorem. 


THEOREM 6.2 The equUibrium Pi{Yi,Ni,Y 2 ,N 2 ) is locally asymptotically stable. 
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Proof: In order to show the stability of Pi, we require the variational matrix M corre- 
sponding to the system (6.3) given by 

/mil PiYi e 0 \ 

M = ^ ^ 

0 0 77233 I 32 Y 2 + ^2Cm 

\ 0 0 -~0'2 —'(^2 + ^) / 


whciG 777-11 = PiNi — (2^iYi + 2^1 + Odi + rfi) 3.1ld 
^33 “ P 2 N 2 ~ ( 2 /? 2 ^ + ^2Cm + 2^2 + ^2 + + ^). 

Now the variational matrix Mi at the equilibrium point Pi(li, TVi, ¥ 2 ,^ 2 ) is given by, 


+ 0iYi 


Ml = 


V 


e 0 

-O'! —di 0 9 

0 0 -(^2^2 + ^^) y025^2+A2C„ 

0 0 —02 


-(d2 + ^) / 

The characteristic polynomial corresponding to the above matrix is given by 


evo 


BYo 


"0 + {^{Yi -t- -f- di)'0 -h di{^iYi -1- -^) -1- aiPiYi / s + {hY2 + X^C, 


Fi 


Fi 




N. 


+(^2 + 9)ljj -h {P 2 Y 2 + A2Cm'^)(d2 + ^) + CX2ip2Y2 + A2Cto) J = 0. 


Since the coefficients in both the quadratics are positive, the local stability of this equi- 
librium point is guaranteed by using the Routh-Hurwitz criteria. 


Nonlinear Analysis and Simulation: 

We first analyze the model (6.3) for = 02 = 0, i.e. the disease related deaths in both 
the populations are zero. In this case, the model (6.3) with the same initial condition 
reduces to the following. 

Yi = PiYi{Ni - Yi) - {ui + di)Yi + 9 Y 2 , 

Ni = Ai-diNi + 9N2, 

Y2 = (,02^2 + X2Cm){N2 - Y2) - (1/2 + ^2)12 - 0F2, (6.8) 

N 2 = A 2 — (^2^2 — dN2, 

where - sq + siEm} and = 

The equilibrium point corresponding to the above system (6.8) is given by i5'*((Ki*, N^, Y^,^^), 
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where 




^2 j,r* _ + 6N^ 

{d2 + ey d, ’ 


_ {/^2iV2 - (^2Cm + Z/2 + Ct2 + 0)} + - {X2Cm + 1^2 + d2 + 6)Y + AhX2CmNi 

■X 0 — — — — — 

2/32 


^ {Aivr - (t^i + rfi)} + - (i/i + di)}2 + 

1 ^ • 

Here we can check that the equilibrium point E*{Y^ ,Ny) is globally stable by 
using the following Liapunov function. 

i' = \(yi - y;? + + k2^(y2 - n*)" + ksi(N2 - JVJ)^ (6.9) 


where ki, k 2 and kz are chosen as follows, 


ki 


AlPlY{ 

dY^d\ ’ 


k2 > 


eYi* 

A2C^iV|’ 


(6.10) 


kz > max 


kiO ^2^2* I 

IT’ ~(d^T&)hPJ^ y 


( 6 . 11 ) 


We further note that the system (6.3) is bounded by the system (6.8), and hence using 
comparison theorems (Lakshmikantham and Leela, 1969), the solution of (6.3) will be 
bounded by the solution of (6.8). Hence we speculate that the nontrivial equilibrium 
point Pi of the system (6.3) may be globally stable. Keeping this in view, we illustrate 
this result by simulation. The system (6.3) is integrated using the fourth order Runge- 
Kutta method and using the following set of parameters in the simulation (Greenhalgh 


1990, 1992). 

/3i = 0.00000051 = /32, i'l = 0.012, di = 0.0004 = dz, 0 = 0.0001, ai = 0.0005, 
Ai = 10 = >l2, A2 = 0.000000021, ^2 = 0.00052, uz = 0.011, Cm = 100000. 
The equilibrium values of Yi, Ni, Yz, and Nz have been found as, 

Yi = 1816.779, Ni = 26672.209, Yz = 4064.687, Nz = 15772.727. 

Simulation is performed for different initial positions, 

In A, ^10 = 6981, iV,o = 20377, 720 = 7301, iVao = 17150. 

In 5,710 = 100, iVio = 10000, 72o = 400, iVzo = 21000. 
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Figure 6.2: Variation of the total population with the infective population in the richer 
class. 



Figure 6.3. Variation of the total population with the infective population in the poorer 
class. 
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population. 

We may note that the last four equations of the model (6.1) govern the growth of the 
infectious disease in the reservoir poor population (as discussed in Chapter 2) affecting 
the rich population in the habitat. Our main aim here is to study the effect of migration, 
particularly 9 on the spread of the disease in the rich population. 

We see that the region of attraction 


I di[d2 + 9) 


0<Y2<N2< 0 < C < 0<E< 




where Cmax = 7 (s — sq + Si^max) and Emax = — is positively invariant and all 
solutions starting in this region T stay in T. The continuity of the right hand sides of 
(6.1) and their derivatives imply that a unique solution exists (Hale 1969). 


6.3 Case I: Q is a Constant Qa 

In this case, we consider that the rate of cumulative environmental discharges Q is in- 
dependent of the population density N 2 . Since Xi + Yi = Ni and X2 -1- T2 = N 2 , it 
is sufficient to consider the following subsystem of model (6.1) with the same initial 
conditions for its analysis. 

n = PiYi{Ni-Yi)-iiy,+ar+di)Yi + eY2, 

Ni = Ai-diNi-aiYi+eN2, 

Y 2 = {l32Y2 + X2C){N2-Y2)-{iy2 + a2 + d2)Y2-eY2, ( 6 . 2 ) 

N2 = A2 ~ d2'^2 “ <^2^^ ~ 

^ — soC" -f- S\ECi 

E = Qa-SoE. 

As discussed in Chapter 2, we use the asymptotic values of C and E in the rest of the 
equations of the system (6.2) and consider further the following subsystem with the same 
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Figure 6.6: Variation of the infective class in the poorer class with time. 



Figure 6.7: Variation of the total population in the poorer class with time. 
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In C, Yio = 5500, Nio = 29000, Fso = 2000, iVao = 12400. 

In A ^10 = 130, Ario = 29100, Fao = 100, A^ao = 15000. 

In Figs. 6.2 and 6.3, we have plotted the total population against the infective population 
of the respective classes. Prom the solution curves, we conclude that the system appears 
to be globally stable about the endemic equilibrium point (Yi, iVi, Ya, N 2 ) for this set 
of parameters. Also in Figs. 6. 4-6. 7 the effect of migration is shown on populations in 
the two regions. It is shown that as 9 increases, the total population and the infective 
population of rich people increase and the total population and the infective population 
of poor people decrease. 

Remark: Mena-Lorca and Hethcote (1992) determined a threshold cr for a single pop- 

ulation without carrier as disease grows in the population only if a > 1, 

otherwise if c < 1 disease dies out and the disease free equilibrium is globally stable. 

Pi 

We note that for above set of parameters, — < 1. Here it is also observed that 
in absence of 9, i.c. (9 = 0, disease would die out in richer class. But for 9^0, disease 
persists in the richer class and regulates the population density. 

6.4 Case II: Q is a Variable 

We consider the case, when the rate of cumulative environmental discharges is dependent 
on the population density of the environmentally degraded region. Then in this case, it 
is sufficient to study the following subsystem of system (6.1) having the same initial 
conditions, 

Yi = l3iYi{Ni-Yi)-{iyi+ai + di)Yi+9Y2, 

Ni = Ai-diNi-aiYi+9N2, 

Y 2 = {p2Y2 + X2C){N2-Y2)-{l^2 + OC2+d2)Y2-9Y2, 

N 2 = A2-d2N2-OC2Y2-9N2, 

C = sC{1-y)-sqC + siEC, 

Ju 

E = Q„+IN2-SoE. 


( 6 . 12 ) 
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In such a case equation (6.22) reduces to 


(q!2 + d2 + $2) (a - A2-C? ) Yi - [P2A2 - {d2 + 0) {^2 Y (^2 Yd2 + e) 


■(<3^2 + + ^)P GA2 Y2 — X2A2 — F — 0 


We assume that ^2 > >^2^0. For / = 0 this equation should reduce to the previous case. 
Since the constant term in the above quadratic is negative, this implies existence of one 
and only one positive root, say ¥2- 


Under the condition, P2 > ^2^, from (6.15) we find the corresponding value of Yi and 
from (6.13) and (6.14), Ni and N2 are calculated. Thus we get a nontrivial equilibrium 
point E4(fi,iVi, 1^2, iV2, (?,£:) provided /32 > ^G. 

Hence we have four equilibria namely, Ei, E2, E3 and E4, where E2 exists if 
^i{ 9A2 + {d2 + 6)Ai) > di{d2 + 0){ui + 0:1 + di), Ez exists if ^ land 

E4 exists provided /?2 > ^2^. 


6.4.1 Stability Analysis 

We now present the stability analysis of these equilibria in the following. The linear 
stability results of these equilibria are stated in the following theorem. 

THEOREM 6.4 The equilibria Ei, E2and E^ are unstable and the equilibrium E4 is 
stable if 03(0402 — O0O3) — of >0, where Oi, a2and 03 are given explicitly in the proof 
of the theorem. 


Proof: The variational matrix A/ corresponding to the system of equations (6.12) at 
{Y„NuY2,N2,C,E)is^vonhy 

f mu ^lYi 9 0 0 ON 

-oi -di 0 9 0 0 

j[^ — ® ”^33 P2Y2 + X2C ^2(^/2 ~ I2) 0 

0 0 -02 -(^2 + ^) 0 0 ’ 

0 0 0 0 s-so + Si£^-^C SiC 

V 0 0 0 / 0 -60) 
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where mu = {N^ -2Yx)- + di), mga = p 2 {N 2 - 2 Y 2 ) - X 2 C - {v 2 + ^2 + + 0)- 

The variational matrix Mi at equilibrium point 


0 , 


Al (^2 + d)+ 9A2 ^ A2 ^ Qo + l-^g 


di{d2 + ff) 


. 0 , 


6,2 9 


, 0 , 


is given by 


Ml = 


where 
mu = /5i 


( mu 

-tti 
0 
0 
0 

V 0 


0 

-di 

0 

0 

0 

0 


9 

0 

mzz 

-a2 

0 

0 


0 

9 

0 

0 {s — 

1 


0 

0 

\. A, 

^ d2+0 
0 


u 

So + U-(Qo + 

0 


0 \ 
0 
0 
0 
0 

— J 


Ai{d 2 + 0 ) + 9A2 


I - (z/i + Qi + di) and mas = ^2-^ - {y2 + 012 + ^2 + 0). 


^ di (^2 + ^) 

It is easy to see that the matrix Mi has one positive eigenvalue given by 


{s — So + Hence Ei is unstable. 

The variational matrix M 2 at the equilibrium point E 2 is given by 


f-PiY* 


M 2 = 


V 


-q;i 

0 

0 

0 

0 


PiY{ 

-di 

0 

0 

0 

0 


9 0 

0 9 

m^s 0 

— q;2 — (<^2 + ^) 

0 0 

0 I 


0 

0 


\ At 
Mzrr 


0 


s — 


So + §^(Qo + I 

0 



0 \ 

0 

0 

0 ’ 
0 

— (^0 / 


where m^g = - {^^2 +a 2 + d 2 + 9). The above matrix has one positive eigenvalue 

as {s - So + |^(( 3 o + implying instability of the equilibrium E 2 . 

The variational matrix M 3 at equilibrium point Ez is given by 


Mz 


/mn 

-ai 

0 

0 

0 

V 0 


PiYr 

-di 

0 

0 

0 

0 


9 

0 

-hY2*- 

- 0!2 

0 

0 


0 
9 

P 2 Y 2 ** 

■(<^2 + 9) 

0 s 

1 


0 

0 

x2{Nr - Yr) 
0 

So + SiE** 
0 


0 
0 
0 
0 
0 

-SoJ 


where miJ = — {PiY^* + ^^)- 

Here also one eigenvalue, which is s - sq + s^E**, is positive, so the equilibrium Ez is 
unstable. 
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The variational matrix M4 at equilibrium point £'4 is given by 

/mil PiYi e 0 0 ON 

—O'! —di 0 9 0 0 

^ _ 0 0 m 33 ^2^ + -^ 2 ^ A 2 (iV 2 — 0 

0 0 -a 2 -(^2 + ^) 0 0 ’ 

0 0 0 0 -fC SiC 

V 0 0 0 I 0 -60 J 

where m'u = - and m'33 = - (P2Y2 + . 

Two characteristic roots of the matrix M4 are given by 




3iYi -1- 9-^ + di'j tp -h (piYi -f 


di + (XiPiYi = 0 , 


and the other four characteristic roots are given by 


■ 0 ^ + CL^Ip^ -h ( 220 ^ "f <^10 + Go “ 0 , 


where 


<23 = P 2 Y 2 + X 2 + d2 + 9 + —C + So, 

1 2 ’L/ 


(22 — ( ^2^2 + A2 


CN2\ , , . ^ . s /s 

0 N 


(<^2 + ^ + + ( 5 o) + (^2 + 0 ) 


+—CSo -h a2(/?2i2 + ^2^), 


CN2 


«! - ^^25^2 + ^2-^ j {^2 + 9)^C + {d2 + 9)^C5 q + |-C75o (^2^2 + \2^^ 


Yo 


+< 5 o { ^ 2 ^ + A2 


CN 2 


{d 2 -I- 0) + (X 2 {P 2 Y 2 + X 2 C) + 5o] , 


Go - ^y32y2 + ^2^^^ (^2 + 9)jC5o + a2{02Y2 + \2C)^C6o + a2X2{N2 - Y 2 )siCl > 0. 


Here 03 > 0 and 


03 Ol 
1 C2 


> 0, thus using the Routh-Hurwitz criteria we observe 
that the equilibrium point £4 is locally asymptotically stable if the following conditions 
are satisfied, 


03 > 0, 


<23 <2l 
1 02 


> 0 , 


03 Oi 0 

1 02 Oq 

0 03 Oi 


> 0 , 


03 Oi 0 0 

1 02 flo 0 

0 03 oi 0 

0 1 02 Oq 
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Figure 6.8: Variation of the total population with the infective population in the richer 
class. 


We see that the first two inequalities are obvious. Also if the third inequality is satisfied 
so is the fourth one as Oo > 0. Hence the equilibrium E 4 is locally asymptotically stable 
under the conditions stated in the theorem. 

Nonlinear Analysis and Simulation: 

As before we speculate that the system (6.12) is globally stable about the equilibrium E 4 
under local stability conditions in the interior of the region of attraction. To illustrate 
this, the system (6.16) is integrated using the fourth order Runge-Kutta method using 
the following set of parameters in the simulation (Greenhalgh 1990, 1992), which satisfy 
local stability conditions. 

Pi = 0.00000051 = P 2 , = 0.012, di = 0.0004 = da, 0 = 0.0001, ai = 0.0005, 

Ai = 10 = A 2 , h = 0.000000021, tta = 0.00052, 1^2 = 0.011, L = 100000, 
s = 0.9, So = 0.6, Si = 0.000002, Qo = 20, So = 0.001, I = 0.00005. 

Here all parameters are in per day except the carrying capacity L and 1. 

Simulation is performed for different initial positions. 

In 1, Vio = 2580, iVio = 28000, 720 = 4000, ATao = 10600, Co = 8843, = 473. 

In 2, Yio = 520, Nio = 27590, ^20 = 3100, N 20 = 23015, Co = 700, Eo = 3000. 
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Figure 6.9: Variation of the total population with the infective population in the poorer 
class. 



Figure 6.10: Variation of the infective population in the richer class with time. 
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Figure 6.11: Variation of the total population in the richer class with time. 



Time in years ^ 


Figure 6.12: Variation of the infective population in the poorer class with time. 
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Figure 6.13: Variation of the total population in the poorer class with time. 


In 3, Vio = 2520, Mo = 24902, F 20 = 410, Mo = 7155, Co = 1230, M = 500. 

In 4, Yio = 1000, Mo = 24020, F 20 = 500, Mo = 19200, Co = 300, Eq = 3400. 

The equilibrium values of M, M, Y 2 , M, C and E have been found as 
Yi = 1910.05, M = 26937.81,^2 = 2594.69, M = 17301.51, C = 37970.01, .C = 20865.07. 
In Figs. 6.8 and 6.9, we have plotted the total population against the infective population 
of the respective classes. From the solution curves, we conclude that the system appears 
to be globally stable about the endemic equilibrium point (Yi, M, Y 2 , M, C, E) 
for the considered set of parameters, provided that we start away from other equilibria. 
In Figs. 6.10-6.13, the effects of 0 on the infectives and the total population of the 
respective classes are shown. It is noted that as 9 increases both the total population 
and the infectives in the environmentally better habitat increase while these decrease in 
environmentally degraded region, as expected. 

Remark 1: The case when the migrating poor population works as service providers to the 
rich population and does not immigrate. 

Consider the case when the migrating population living in a environmentally degraded 
region works as a service provider to the rich population living in a environmentally 
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cleaner region. Using the same notation as before, the model can be expressed as follows, 

— -^1 ~ — PiXiYi — P12X1Y2 + I'lFi, 

Yi = piXiYi + 1312X1Y2 — {i/i + Q!i + di)Yi, 

X 2 = A 2 — ( 12 X 2 - {P 2 Y 2 + \2C)X2 + P 2 Y 2 , (6.17) 

Y2 = {P2Y2 + X 2 C)X 2 — (1^2 + 0:2 + ^2)1^2, 

^ — SqC + siEC, 

E = Q{N2)-SoE. 

Since there is no immigration, 9 is taken as zero and P 12 is the rate of interaction with 
the infected poor population. In writing down the model it is assumed that the infected 
rich people do not interact with the service providers. 

The analogy of the model (6.17) is similar to that discussed in Chapter 7, and hence 
it is not included in this chapter. The corresponding results can be visualized from the 
discussions and results presented in Chapter 7. 

Remark 2: The same analysis can be done for the model with logistically growing 
human population. For a single population the logistic case has been done in Chapter 2, 
so one can easily extend this to the case of two populations as has been done in the case 
of constant immigration. 


c = 


6.5 Conclusions 

In this chapter, a non-linear mathematical model is proposed to study the effect of immi- 
gration of population, from an environmentally degraded habitat to the environmentally 
better habitat from sanitation point of view, on the spread of an infectious disease. This 
model is analyzed for two cases, (i) when the cumulative rate of environmental discharge 
is constant and (ii) when this rate is variable. In case I, we have shown that the non- 
trivial equilibrium point is locally asymptotically stable. Also by comparing the model 
of case I to the model when the disease related death rates are zero, we conjecture that 
this nontrivial equilibrium point is globally stable showing the endemic nature of the 
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and poor live in nearby habitats such as in big cities of the third world countries, then 
the poor people work as service providers in the houses of rich people but do not settle 
in the habitat of rich people. This population of poor people, as service providers, play 
an important role in the spread of infectious diseases as they carry pathogens in or on 
their bodies and may also transport disease vectors, such as lice. Organisms that survive 
primarily or entirely on a human host and are spread through sexual contact, droplet 
nuclei and close physical contact, can be readily carried to the susceptible population by 
these service providers, for example, AIDS, tuberculosis, measles, pertussis, diphtheria 
and hepatitis B (Morse 1995). 

In this chapter, therefore, the spread of epidemics in a socially structured population 
(rich and poor) living in two nearby habitats, one cleaner than the other, is modeled 
and analyzed. It is assumed that the poor people suffer from an infectious disease (such 
as typhoid) caused by bacteria, the growth of which increases due to poor sanitation 
and various kinds of household discharges such as waste milk products, food, etc. It is 
assumed that the poor people living in the environmentally degraded habitat spread the 
disease by interacting with the rich population by working in their habitat as service 
providers. It may be noted here that the uncontrolled household discharges into the 
environment, not only increase the population of bacteria, but also affect the human 
population in several ways related to their general state of health, physical vitality, nat- 
ural or acquired immunity and thereby making them more susceptible to various kinds 
of infectious diseases. 


7.2 SIS Model with Immigration 


We consider a general SIS model for infectious diseases caused by bacteria in a socially 
structured population (rich and poor) living in two adjoining habitats or neighborhoods. 
The environment where the rich people live is much cleaner, while the environment in 
the neighboring habitat where poor people live is not so clean and is very conducive 
to the growth of bacteria population caused by uncontrolled household discharges. It is 
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assumed further that infected poor people interact with rich people by working as various 
kinds of service providers and infect them. It is assumed further that the infected and 
sick rich people interact only with susceptibles of the rich class but not with people of the 
poor class. As in the previous chapter, here also the total population density of the rich 
class Ni is divided into the susceptible class Xi and the infective class Yy. Similarly the 
total population density of the poor class is divided into the susceptible class X2 and 
the infective class Y2. It is assumed that all susceptibles X2 of the poor population are 
affected by the bacteria whose density B2{t) grows logistically with given intrinsic growth 
rate and carrying capacity. It is considered that the growth of the bacteria population also 
increases due to increase in the cumulative density of environmental discharges, caused 
by sources independent or dependent on the human population in the environment. It is 
assumed that the bacteria population does not survive in the clean environment of rich 
people and only affects the population in the degraded environment of the poor class. In 
view of this, the model can be written as follows: 

= Ai-diX,-^,X{Yi-XiX,Y 2 + v{Yu 
n = PrX{Yi+X^X^Y2-[vx^a^ + di)Yu 
N\ = Ai — diNi — ojiYx? 

A2 = A2-d2X2- ^2X2Y2-X2X2B2+V2Y2, 

Y2 = P2X2Y2 + X2X2B2 — ( 1^2 + 0:2 + d‘^Y 2 i 

N2 = A2 - ^2^2 - 0:2^2, (7-1) 

+ S2Y2 — S20B2 + S3B2E, 

E = Q{N 2 )- 5 oE = Qo + IN 2 -SoE, 

Xi(0)=Xio > 0, yi(0) = yio>0, Ni{0) = Nio > O, A2(0)=X2o > o, 

^2(0) = y20 > 0, N 2 { 0 ) > 0, .62(0) > 0, E{ 0 ) = Eo> 0 . 

Here E{t) is the cumulative density of household discharges conducive to the growth of 
the bacteria population; Ai and A2 are the immigration rate constants of the human 
population into the rich and poor populations respectively; di and d2 are the natural 
death rate constants corresponding to rich and poor populations; and Ai are the dis- 


B2 


SB2 (^1 


L 
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ease transmission coefficients in the rich population due to the infectives of the rich and 
poor populations respectively; P2 ^ind A2 are the disease transmission coefficients in the 
poor population due to the infectives in the poor population and bacteria respectively; 
tti and 0:2 are the disease related death rate constants corresponding to rich and poor 
populations respectively; s is the intrinsic growth rate of the bacteria population; L is 
the carrying capacity of the bacteria population in the natural environment; S20 is the 
death rate of bacteria due to control measures, where s >820', S2 is the rate of release 
of bacteria from the infective population; S3 is the rate of growth of the bacteria popu- 
lation due to the environmental discharges; Q{N2) is the rate of cumulative environment 
discharges conducive to the growth of bacteria into the poor population which depends 
on the density N2 of the poor population and (5o is the depletion rate coefficient of the 
cumulative environmental discharges. This model is analyzed for two types of environ- 
mental conditions: 

(i) the rate of cumulative environmental discharges Q is a constant Qa and 

(ii) the rate of cumulative environmental discharges Q is a function of the poor population 
density N2, assumed in the form Q = Qo + IN2, I being a positive constant. 


7.2.1 Case I: Q is a Constant Qa 

Since Afi -f Ti = Ni and X2 +Y2 = N2, the system (7.1) is equivalent to 


= PliN^-Y,)Y^ + X^{N^-Y,)Y2-{u^ + a^+d^)Y^, 

N, = Ai-diNi-a^Y^, 

Y 2 = ^2{X2 — Y2)Y2 + X2{N2 — Y2)B2 — (1^2 + CX 2 + (^ 2 ) 12 , 
N 2 — -^2 ~~ d2N2 — OC‘^2i 
B2 — SB2 ^ -b S2I2 ~ S20B2 + S3B2E, 

E = Q^- 5oE. 


(7.2) 
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Prom the last equation of the above system of equations, we note that the asymptotic 
value of E is 


lim E(t) 

t^oo ^ ' 


Qa 

5o' 


Therefore, in this case it is sufficient to study the following subsystem of (7.2) for its 
equilibrium analysis. 


Yi = Pi{Ni-Yi)Yi + Xi{Ni — Yi)Y2- {ui + ai + di)Yi, 
iVi = Ai-diiVi-aiyi, 

Y = ^ 2(^2 — Y2)Y2 + X2{N2 — Y2)B2 — {l>2 + OC2 + d2)Y2, 

N 2 = A 2 — d2N2 — 0:2^2, (7.3) 

B 2 = sB2(i-^'^+S2Y2-S20B2 + Sz^B2. 

The region of attraction T of the above system is 


T = {(Ki, Nu V,, iV2. Bj) : 0 < n < JVi < ^, 0 < y, < ^. 0 < B, < 


where B 2 max = 


£ 

2s 


r 

f , , 

f , Qa\^ 

4sS2^2 


|S-S20 + S3^| 

d2L 


is positively invariant and all solutions starting in this region T stay in T. The continuity 
of the right hand sides of (7.3) and their derivatives imply that a unique solution exists 
(Hale 1969). 

The result of equilibrium analysis is stated in the following theorem. 


THEOREM 7.1 There exist the following three equilibria of the model (7.3), 
(i) £^1(0, ^, 0,^,0), 


(ii) £'2(17, , 0, 0), which exists if > {I'l + cxi + di) , 


where Yf 


“ (r/i + Q!i + di) 

A(l + if) 


and Ni — 


Ai 


-cciYf 
di ’ 


(Hi) Ez{Y,NuY2,N2,B2). 
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Figure 7.1: Existence of equilibrium point. 


Proof : The existence of either of the first two equilibria is obvious. We prove the 
existence of the third equilibrium Ez by the isocline method. Setting the right hand side 
of the equations in (7.3) to zero and with some manipulation, one gets: 


Ai — aiYx A2 — OL2Y2 

/Vj _ , JS2 = 


di 


d. 


Yo = 


— [1 
S2 [ l ' 


— ( S — S20 + 53 


Qa 


Bo 


(7.4) 

(7.5) 


+ 1) ("■ + “■ (1 + 1) = 0. 


(7.6) 

“ (l ‘2 + 0:2 + (h) — A 2 (1 + ^) B 2 I ^2 - = 0. (7.7) 

Now we show the existence of Y 2 and B 2 from (7.5) and (7.7) and the corresponding 
values of Yi,Ni and N 2 can be obtained from (7.4) and (7.6). 

Prom equation (7.5): we have, 

(i) ^2 = 0 or P 2 = — — S 20 + ^3^1 = B 2 > 0, when = 0, 

(ii) the slope of the curve (7.5) is given by 


P2 


■•S)« 


f P2A2 

I d2 
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at (0, 0), the slope of (7.5) i.e. is negative and at (BajO), it is positive and 
from (7.6), the slope of (7.5) is zero at B 2 = i(s - S 20 + Sz^) - 
(iii) the slope is increasing as (g|) > 0. 


Prom (7.7), the following points are observed, 
(i) For B 2 = 0, 


Y 2=0 or y, = = 

/l I rvn \ 


^ 2(1 + %) 


'■ ^2 (say). 


which is positive if ^ 2 ^ > ( 1^2 + 02 + ^ 2 ) and negative otherwise, 
(ii) At (0, 0), the slope of (7.7) h 


i.Q crivpn 


dY^ 

dB2 ^2^ — {y2 + 02 + rfr) ’ 

which is positive or negative depending upon Y 2 being negative or positive respectively. 
At (0, ^ 2 ), the slope is given by 


' ^ 1^2 + 02 + (^2 

dY2 __ ^2 

1^2-^ ~ {^2 + Ol2 + d2) 

0.2 

which is positive or negative depending upon F 2 being positive or negative respectively, 
(iii) Slope ^ is decreasing as (^■) < 0. 

Thus, after plotting Y 2 with B 2 corresponding to (7.5) and (7.7) in Fig 7.1, we see that 
there are two intersecting points (0,0) and (B 2 , F 2 )- After finding ^ 2,^2 and Yi, we 
can calculate Ni and N 2 using (7.6) and (7.4). Also, we note from (7.6) and (7.7) that 
Ti < and Y 2 < Thus corresponding to (152,1^2), we have the third 

nontrivial equilibrium, namely E 3 {Yi,NijY 2 ,N 2 ,B 2 )- 


7. 2. 1.1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results of 
these equilibria are stated in the following theorem. 

THEOREM 7.2 The equilibria Ei and E 2 are unstable and the equilibrium E 3 is locally 
asymptotically stable. 
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Proof : The variational matrix M at (Fi, Ni,Y2,N2, B2) is given by 


M 


where 


and 


/mil 

piYi + XiY2 Ai(iVi-Fi) 

-Oi 

-di 

0 

0 

0 

mzz 

0 

0 

-Q!2 

0 

0 

S2 

1-2F1) 

— {X 1 Y 2 + + ai + di)j TTizz 


0 

0 


0 

0 




P2Y2 + X2B2 M{N 2 — Y2) 


-d 2 

0 


0 

^ 7^55 


J 


^2(^2 ^ 21 ^) ~ {^ 2^2 + + <^2) 


17155 = s - S20 + f ^2-, 

The variational matrix Mi at the equilibrium point Ei, is given by 


Ml = 


/mil 

0 


0 

-Oi 

-di 

0 

0 

0 

0 

mzz 

0 

0 

0 

-0L2 

— <^2 

0 

0 

S2 

0 


0 

0 

0 


S — S 20 «+ ^ 3 ^/ 

It is easy to see that at least one eigenvalue of Mi is positive or has positive real part. 
Hence the equilibrium Ei is unstable. 

The variational matrix M2 at the equilibrium point E2, is given by 


M, 


/-^iF; /?iFi* Ai(iV*-Fi*) 0 

-oi -di 0 0 

0 0 1 ^ 2 ^ — {1^2 + 02 + <^2) 0 

0 0 —02 —d2 

V 0 0 S2 0 


0 

0 




s — S20 + 53^ / 


The characteristic polynomial corresponding to the above matrix is given by 


(^2 + + {di + PiYi)'tp + {di + ai)l3iY { } + hiip + / 12 ] = 0 , 

where 

~ ~ (^2 + 0:2 + ^2) + s — S20 + 53^ I 

Using the Routh-Hurwitz criteria, this equilibrium is stable, provided 
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We note that both the above conditions are simultaneously not satisfied, so this equilib- 
rium is unstable. 


The variational matrix 


M3 = 


/mil P 1 Y 1 + X 1 Y 2 Ai(iVi-yi) 0 0 

-di 0 0 0 

0 0 m33 ^ 7 X 2 + X 7 B 2 A2(iV2 — Yt] 

0 0 — Q;2 — ^2 0 

V 0 0 S2 0 mss 




)■ 


^ ^ 52 0 ?rf55 / 

where m'u = - (A^i + mzz = - [^ 2 % + and m^ = - (|;52 + S2^ 

The characteristic polynomial corresponding to the above matrix is given by 

> , AiiVi^\ , . /^ AiiViyAl r ,3 ,2 , , , ] 

' H + ai'ijr + a 2 'tp + 03] = 0 


I V-" + [d, + ftr, + j V + <i, f Afi + 


(7.9) 


where 

Qi 

02 

as 


d2+/32F2 + ^^^ + fB2+52^ > 0, 

12 ^ B 2 

dj + ;^S 2 + S 2 ^ j + (^*^2 + 

+a2(/?2^ + '^2^2) ■“ ^ 2(^2 ^2)52 > 0 ^ 

d2 (P2Y2 + + 52 ^^ + a2{P2Y2 + X2B2) [jB2 + 52 ^ 

— X2{N2 — Y2)d2S2 > 0 . 


^2' 

jB2 + .22^_ 


In equation (7.9), the coefficients of the quadratic are positive which implies that the 
roots cannot have positive real parts. Also it can be checked that 0102 — 03 > 0, hence 
by the Routh-Hurwitz criteria, the roots of the cubic in equation (7.9) have negative real 
parts. Hence the system is locally asymptotically stable about this equilibrium E 3 . 


Nonlinear Analysis and Simulation: Before we proceed for simulation, it can be 
shown heuristically that the system (7.3) is globally stable about its nontrivial equilib- 
rium, when oi = 0 and 02 = 0 (see Appendix III). Since the system (7.3) is bounded 
by its corresponding system with Oi = 0 and 0:2 = 0, using a comparison theorem (Lak- 
shmikantham and Leela 1969), it is concluded that the solution of the system (7.3) is 
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Figure 7.5: Variation of Fi and Y 2 with time for different rates of release of bacteria from 
the infective population. 



Figure 7.6: Variation of Yi and Y 2 with time for different rates of growth of bacteria 
population due to the environmental discharges. 
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Figure 7.7: Variation of Yi with time for different disease transmission coefficients due 
to infectives of poor population. 

these figures, 

(i) if the bacteria population increases due to household discharges or otherwise, the 
spread of the infectious disease in both the rich and poor populations increases, 

(ii) if infectives in poor population increase, the disease spreads faster in the rich popu- 
lation and 

(iii) if the disease transmission coefficient between the infective poor population and the 
rich population increases then so does the number of infectives in the rich population. 
Therefore, it is in the interest of rich people not only to clean their own habitat but also 
the nearby habitat of poor people. 

7.2.2 Case II: Q is a Function of N 2 

As before, in this case also, it is sufficient to consider following subsystem of the original 
system of equation (7.1) with the same initial conditions, 

Yi = 0i{Nx-Yx)Yi + Xx{Ni-Yi)Y-i-{ui+a^+di)Yu 
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ATi 

% 

N2 

B2 

E 


Ai - diNi - aiFi, 

I32{N2 - ^^2)^2 + ^2(^2 - Y2)B2 - {^2 + Q!2 + C?2)^2, 
A2 — d2N2 — CX2Y2, 

SB2 ^2^2 ~ S20B2 4 " S3B2E, 

Qo + IN2 — SqE. 


(7.10) 


The region of attraction T' 


r = \yi,Ni,Y2,N2,B2 ,E):Q<Yi<Ni<'^, Q<Y 2 <^, 0 <B 2 <B 2rmx, 

Qq + 


Ao 


0<E < 


So 


where Samax = 


£ 


S-S20 + Ss-^j + ^ 


s - S 20 + S 3 ' ■ > + 


, 4552^2 

"X” J d 2 L 


is positively invariant and all solutions starting in this region T' stay in T' . The continuity 
of the right hand side of (7.10) and its derivatives imply that a unique solution exists 
(Hale 1969). 

The result of equilibrium analysis is stated in the following theorem. 


THEOREM 7.3 There exist the following three equilibria, namely 

di d2 do 


(it) E, I y,*, ^,Q, — — , which exults > (t/i+fti+rfi), 

y.* = ^ gi , + . °- + <^0 . 

/3i(di + ai) /Si(o:i+di) 


i) EzfYi, Ni,Y 2 j ^ 2] B2, E), which exists if $2 > > 


, ,5 L , Qo S3IA2 . 

where B2 = —{s - S20 + S3— + 

S do doU2 
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Proof. Existence of either of the equilibria Ex or is easy to show. Hence we only 
prove the existence of Ez by the isocline method. Equating the right hand sides of the 
above equations in ( 7 . 10 ) to zero and after some manipulation we get 


Nx = 7V'2 = “^2 ~ O-TX-i ^ Qo + IN 2 


dx 


do 


( 7 . 11 ) 


^ ~ + Ai (1 + = 0 , 


p2 (1 + Yi - - {u2 + 02 + dz) - A 2 (^1 + Saj F 2 - A2^B2 = 0 (7.13) 


( 7 . 12 ) 






„ „ I Qo S3I yl2\ JO 

S — S20 + S3— h B2 

K do do 02 / 


S 3 ^a 2 „ 
52 - -y-j-Bz 
0 ott 2 


( 7 . 14 ) 


We now show the existence of 1^2 and B2 from ( 7 . 13 ) and ( 7 . 14 ). Then Ti, Ni, Nq, and E 
can be obtained from ( 7 . 11 ) and ( 7 . 12 ). The equation ( 7 . 13 ) is same as (7.7). From 
( 7 . 14 ), we have, 

(i) when ^2 = 0, I2 = 0, 

(ii) when y2 = 0, H2 = 0 or B2 = —{s - S20 + 53% + = B2, 

s do do02 

(hi) the slope of this curve is given by 


2s 


dY 2 

dBo 


(iv) at (0, 0), the slope is — 


B2 — S — S20 + 53 -r h 


Qo , 53^2') , S 3 la 2 


Sq ^ 0^2 j <^od; 


-F2 


( S 3 la 2 

V" 5 od 2 \ 

S - S20 + 53 ^ + 


52 


(v) at the point (B2, 0), the slope is given by 


< 0, 

^B< 


1^2 


_ sslajBi 

*2 Sod 2 


=-, which is positive for S2 > 


We make this assumption because for / = 0, we must get the slope at {B2, 0) positive as 
in Case I. 

As before plotting curves ( 7 . 13 ) and ( 7 . 14 ), we get two intersecting points ( 0 , 0 ) and 
(B2, Y2). Also it can be checked that Yx < and Y2 < 5^ and B2 increases with 
an increase in Qq. Thus corresponding to {Bz^Y^, we have a third equilibrium point as 
E 3 (YuNuY!,,N 2 ,B 2 ,E). 
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7.2.2. 1 Stability Analysis 


In the following we discuss the stability of the system (7.10). The local stability results 
of these equilibria are stated in the following theorem. 


THEOREM 7-4 The equilibria Ei and E^ are unstable and the third equilibrium is 
locally asymptotically stable provided, as(aia 2 - aoUs ) - af > 0 , where ai, 02 and 03 are 
given explicitly in the proof of the theorem. 


Proof: The variational matrix M corresponding to system (7.10) is 


/ ^11 ATi + A1I2 Ai(iVi — li) 


M 


-ai 

0 

0 

0 

V 0 


—di 

0 

0 

0 

0 


0 

^33 

- 0!2 

S2 

0 


0 0 

0 0 

^2^2 + A 2 B 2 X2{N2 — Y 2 ] 

-d2 0. 


0 

I 


where 

mn = PiNi - - AiFa - (i/j + oi + dj, 

^33 = P2N2 - 2P2Y2 - X2B2 - (1^2 + 02 + ^2) and 
m55=S- S20 + S3E - ^B2. 

The variational matrix Mi at the equilibrium point Ei is 


/At 


Ml = 


+ Qfi + dj) 0 

~0;i ~di 

0 0 

0 0 

0 0 

0 0 


+ it 


di 

0 

At ~ (^2 + 02 + ^2) 


-02 

52 

0 


W55 

0 


0 

0 

0 

—d2 

0 


0 \ 
0 
0 
0 

S3B2 

“A J 


0 

0 

d2 

0 

rrtss* 

0 


0 \ 
0 
0 
0 
0 

~SoJ 


where = s - S20 + S3 , 

Now it is noted that at le^t one eigenvalue of Mi is positive or has positive real part, 
implying instability of the equilibrium Ey. 
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The variational matrix at the equilibrium point E2 is 

— Oi —di 

0 0 


M2 = 


where E* = 


0 

0 

0 

Qo + ^ 


XiN* 0 

0 0 

^2^ - {1^2 + 02 + d 2 ) 0 

0 -q ;2 -d 2 

0 S 2 0 

0 0 I 


0 

0 

^2^ 


0 

S - S 20 + S 3 E* 
0 


0 \ 
0 
0 
0 
0 

—^oj 


d 2 


So 


The matrix M2 also has at least one positive root or root with positive real part. So the 
equilibrium E2 is unstable. 

The variational matrix M3 at the equilibrium point Ez is 

/mil ^ 1 X 1 + XiY Ai(iVi — Fi) 0 0 


Mz = 


-Oil 

0 

0 

0 

V 0 


-d, 

0 

0 

0 

0 


0 „ 0 „ „ 
nfss P 2 Y 2 + A2.B2 A2(iV2 — Y^ 

—OL2 — (^2 0 

S 2 0 mss 

0 I 0 


0 \ 
0 
0 
0 ^ 
S3B2 

-Sol 


where m'n — — (^iFi + mss = —{P2Y2 + _I_ S 2 ^). 

The characteristic polynomial corresponding to the above matrix is given by 

{"0^ + (/ 3 iFi H — + di)‘>P + di(^iYi -i — ^ ) + o:i(^iYi + AiF)} 

Yi 11 

X {'0^ + 03 ?/^ + a2iP' + aX + ao} = 0, 


where 


as 


02 


P^Y 2 + ^^^^^+d 2 +So + yB 2 + S 2 ^ > 0 , 

F2 ^ B2 

[P 2 Y 2 + (o{2 + 5o + ^B2 + ^2^) + d2 + J^B2 + ^2^ 


<^1 


(_^2 + S2^)<Jo + Q:2(y^2^ + A 2 .S 2 ) — S2A2(iV2 — ^) > 0, 

+a 2 (; 025^2 + A 2 B 2 ){^O + |;^2 + S 2 ^}-S 2 A 2 (iV 2 -y 2 )(rf 2 + ^o) > 0 
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Figure 7.8: Variation of Ni with Yi. 


Oq — ( P 2 Y 2 + 


\ 2 N 2 B 2 
Y 2 . 




+ Oi2(^2^2 + A2B2)5 o [ + S2“^ 

B 2 . 


— A2(A'2 ~ ^){Q^2^53.B2 + S2<^2<^o} > 0- 


By the Routh-Hurwitz criteria, Ez is locally asymptotically stable if following conditions 
are satisfied. 



1 02 


03 Cl 0 


03 > 0, 

V 

0 

1 O 2 Oo 

0 az Oi 

V 

0 


03 Oi 0 0 

1 02 Oo 0 

0 03 Oi 0 

0 1 O 2 Oq 


> 0. 


The first two inequalities are obvious. If the third is satisfied, then so is the fourth one 
as Oq > 0. Hence Ez is locally asymptotically stable under the condition mentioned in 
the theorem. 


Nonlinear Analysis and Simulation: As before it is speculated that system (7.10) 
may be globally stable, provided that we start away from other equilibria (see Appendix 
III). To show this, the system (7.10) is integrated using the fourth order Runge-Kutta 
method and using the same parameter values as in Case I with Qq = Qa and an additional 
parameter value I = 0.0005. The equilibrium values of FI, Ni, Yz, N 2 , B 2 and E 
have been found as 


Yi = 757.59, Ai = 24053.01, Y 2 = 3013.64, N 2 = 21082.28, B 2 = 2255803.38, E = 30541.10. 
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each case is globally stable under its local stability conditions, provided we start away 
from other equilibria. By simulation, it is observed that if Ai, i.e. the interaction between 
rich and poor people increases, the infective population of the rich class increases as 
expected. Also as other parameters such as the growth rate of the bacteria population 
or the rate of release of bacteria from infectives or the rate of growth of the bacteria 
population increases, infective populations of both the classes increase. This suggests 
that the spread of the infectious disease in rich people living in a better environment 
increases due to the interaction with service providers, who are living in relatively poor 
environmental conditions. This further suggests that rich people must involve themselves 
economically and otherwise to improve the habitat of service providers. 
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Figure 7.13; Variation of Yi with time for different disease transmission coefficients due 
to the infectives of the poor population. 

In Figs. 7.9 and 7.10, we have plotted the infective population against the total popu- 
lation of respective classes for different initial positions 1, 2, 3 and 4. From the solution 
curves, we conclude that the system is globally stable about the endemic equilibrium 
point {Yi, Ni, Y 2 , N 2 , B 2 , E). The effects of various parameters on Yl and Y 2 are 
shown in Figs. 7.11-7.14 and similar results have been found as in Case I. 

7.3 Conclusions 

In this chapter, a nonlinear mathematical model is considered for infectious diseases 
caused by bacteria in a socially structured population (rich and poor) living in two nearby 
habitats, rich class living in a cleaner region and the poor class living in a neighboring 
region which is not so clean due to household discharges of poor people causing the growth 
of bacteria. Two cases are considered, (i) the rate of cumulative environmental discharges 
is a constant and (ii) the rate of cumulative environmental discharges is population 
dependent. In each case the existence of equilibria is shown and their local stability 
results are discussed. It is shown by simulation that the endemic equilibrium point in 
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Using Sylvester criteria we get the following inequalities for negative definiteness of V 
4k2 {I'l + + 5i) (d2 + ^ 




<^1 


t2 


I3i\Ni — [1 + — ]Yi \ + k2\{P2 + A 2 — )(iV 2 - U 2 ) 


di 


and 4^iY2kidi > {^iY 2 f. 


From the equations for the equilibrium of the above system, we get 


jftiVi - ft(l + ^)yA |(& + >'1^)(N2 - n)! = {vi + d, + 


Thus the first inequality reduces to 


r 

4(z/i + dll + Si)k2P2^2~^y^l > 

di 


- ^ 1(1 + |-)yi| - k2 Up2 + ^2^){N2 - Y2) 


l2 


— /5i(l + 

Now choosing k 2 = 7 z — the last inequality gives Yi > 0. 

W2 + X2^){N2-Y2) 

Further taking the maximum of the right hand side of the second inequality gives 


max j 

ki > — - — di. 

Thus choosing ki and k 2 in the above manner, we get V negative definite implying global 
stability of the reformulated system in the interior of the region of attraction. Hence the 
global stability of the original system is expected. 

As mentioned earlier the global stability is considered with respect to the interior of the 
region of attraction. 
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Appendix II: Global stability proof of the system (5.4) for cii = 0 


The proof of the global stability of the system (5.4) when ai = 0 in the interior of the 
region of attraction is carried out for a reduced system only. By considering the case 
when Oi\ = 0, the system (5.4) can be written as 

— yi~ Zi)Y 2 — {ui + di + ^i)!! 

= 5iYi-diZi 

Ni = A-d^Ni (1) 

^2 = ^2(iV2-F2)n + A2(iV2-y2)^i-|a2 + d2 + (l-a')^iV2}l2 

N 2 — T 2 N 2 ( 1 ) — Oi2N2 + 52 N 2 E 

E 


\ K2J 
Qo + INi — SqE 


Using the asymptotic values Ni, N 2 and E respectively as follows: 

N 2 = ^ (^2 - 0:2 + S2Ey E = -^ (Qo + y) ) 
and Zi = f^Yi (obtained by putting Zi = 0), we get the following reduced two dimen- 
sional system: 


^1 


/?! I A^i — ^1 + Y^Y2 — {yi + di <5i)yi 
5x' 


(2) 


Y 2 = \p2-^X2J^j{N2-Y2)Y^-d2Y2 

where (J 2 = ^2 + <^2 + (1 - o!)^^N 2 - Now it is easy to observe that the system (2) is 
globally stable by taking following Liapunov function: 

V = \{Y^-Yxf + ^{Y2-Y2)\ 

where . _ , c \ - -i 

_ {a-^ 1 - ft (i + i) u} 

(ft + A2U(^2-^2) ■’ 

and Yi & 1^2 are the equilibrium values corresponding to the system (2), which are the 
same as the nontrivial equilibrium values corresponding to the original system (1). 

[V = - {/3i (1 + I) F 2 + i/i + di + 5i} (5^1 - Zif + [{A^i - a (1 + I)-) ^ 1 } 

+ki (P 2 + A 2 I) {N 2 - Y 2 )] (Fi - Yi){Y2 - Y 2 ) -ki{ [P 2 + X2^) Fi + J 2 } (^2 - y2)^] - 
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Appendix III; Global stability proof of the systems (7 .3) and (7.10) 
for ai = 0, a 2 = 0) 


The system (7.3) is as follows: 


= Pi{Ni-Y,)Y, + X,{Ni-Yi)Y2-{ui+di)Yi 

N^ = A^-d^N, (1) 

Y = l32{N2-Y2)Y2 + X2{N2-Y2)B2-{U2 + d2)Y2 
N 2 = A 2 ^2-^2 

B 2 = 5^2 +32^2-520^2 + 53^52 

Again using asymptotic values of Ni and N 2 bs ^ and ^ in the last system we get 
following reduced system as follows: 


ii = A(^-n)Ki + Ai(^-yi)y2-(i/i + di)Fi 

T2 = ^2(^-T2)r2 + A2(^-y2)52-(l/2 + ^i2)T2 (2) 

B 2 = SB 2 ^1 + 32^2 — S 20 B 2 + 33^52 

Now taking the following Liapunov function corresponding to the system (2) 

V = 2^^^ ~ — B 2 — ) 

we get derivative of F as 


V 


-kr 


m + x, 

A2T2 + A2 


diS 

di Yi 

A 2 B 2 


(n - y,)^ + Ai (^ - y.) (y. - y)(y2 - y) 

(y2 - yf + h\2 (^ - y) (y - y)(B2 - B2) 

^ 2 y ^7 . 1 S2 


-k2r(B2 - B2f - k2^{B2 - B2f + yf (B, - Bj)(y - y) 

^ D2£>2 Bo 


Using Sylvester criteria, V will be negative definite if following inequalities are satisfied: 


(i) 


diYr 2 d2 y. 
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f--\ Al 1 -^2-42 - 62 , Is 

(11) 4fcx--— > 
^ ^2 Y 2 ^ ^ 




(ui) 

4 <*2 y, 2Z, V B2y 

Taking maximum of right hand side of first two inequalities and simplifying all three 
inequalities, we get 


ki > 


AiAiy'id2 


(3) 

(4) 

(5) 

'■2 X2A2B2 S 

We can choose ki satisfying (3). Prom (4) and (5) we get upper and lower limit of so 

2 

if > S 2 , the choice of k 2 is always possible such that (4) and (5) both are satisfied. 

2 Y 2 L 

Thus under above condition the system is globally stable. 


2 dxA 2 .<42.62 

.62 5^2 ki 

2 Y 2 LX 2 A 2 ^ h" 

h 2d2Y2Lsl 
k2^ ' 


Global stability proof of the system (7.10) with aj = 0 and 02 = 0 
In this case also taking asymptotic value of Ni, N 2 before and asymptotic value of 
E as we get system (2) with ^ replaced with Thus similar argument is 

valid for the global stability of the system (7.10) with oi = 0 and 0:2 = 0. 



Appendix I; Global stability proof of the system (5.2) for ai = 0, 
0^2 = 0 


It may be pointed out here that the statement regarding the global stability of the non- 
trivial equilibrium point has been mentioned as a motivation for the numerical simulation 
carried out subsequently. To prove this, we reformulate the system by using the value of 
= ^Yi (obtained by putting Zi = 0) in the rest of the equations of the system. Thus 
we consider the following reduced system of equations for global analysis: 


Ni 

% 


/3i |iVi - (^1 + 72 - [vi + di + 5i)Y, 

A - diNi 

h{N2 - 12)71 + X2{N2 - l2)x^l - 

di 

(ft + Xij] NiYi - ift + Y1Y2 - ftft 


( 1 ) 


where d 2 = d 2 + (1 - a')^N 2 . The equilibrium point (7i,iVi,72) of this system corre- 
sponds to the equilibrium point of the original system. 

To check the stability of (7i, i^i, 7) we use the following Liapunov function: 



Thus V using the above system (1) can be written as follows, 

V — ■|y3i(l -f + cli + ^ (li ■“ li)^ ~ kidi{Ni — Ni) 


+ [|/5iiVi — Pi{l + ^)^i| + ^2 |(/52 + -^2 ^) (.^2 — ^ 2 ) 


( 7 i - 7 )( 7 - 1 ' 2 ) 


-7 1 W 2 + X2^)Y^ + d2| (1^2 - Yf + /?17(7 - 7)(iVi - iVi). 
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We have been guided by the fact that in most tropical third world countries, where pop- 
ulation size changes due to immigration, births and deaths, the population density does 
not remain constant and generally the contact rate is proportional to the size of suscep- 
tibles and infectives (Anderson and May 1979, de Jonge et al. 1995). To make this case 
very specific, one may consider labor colonies in cities like Bombay, Calcutta, Kanpur, 
etc. In these colonies, the area of which may be several square kilometers, working class 
people live with very little municipal facilities such as clean water, bathrooms, toilets, 
etc. Further these people due to lack of civic facilities, throw household wastes in the 
open space on or around the road. A fraction of the garbage in these dumps becomes 
the breeding ground for carriers, vectors or some harmful bacteria. This unhy genic sit- 
uation enhances the growth of carriers, vectors and bacteria, causing the fast spread 
of infectious disease. The population density per square kilometer of such colonies can 
vary from about 50,000 to 1,50,000. In joint families, grown up youngsters get married 
and children continue to live with their parents, thus increasing the population density. 
Further these colonies also work as immigration centers. Laborers from the countryside 
migrate to these colonies and rarely return home. This further increases the population 
density of these colonies. Hence we feel in such situations simple mass action incidence 
may be more appropriate while modelling the spread of infectious diseases, where rate of 
contact may be directly proportional to population density (Mena-Lorca and Hethcote 
1992). 

To be specific, the following problems have been considered selectively to provide sam- 
ples of modelling of environmental and demographic effects on the spread of infectious 
diseases: 

(i) Modelling the Spread of Carrier-Dependent Infectious Diseases with Environmental 
and Demographic Effects, 

(ii) Modelling Bacterial Disease with Environmental and Demographic Effects, 

(iii) Modelling the Spread of Malaria: Environmental and Demographic Effects, 

(iv) Modelling the Spread of Malaria with Human Reservoir; Environmental and Deino- 
graphic Effects, 
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where N = 
and (Hi) P3 


+ (l^ + Q: + d) 

«l+l) 

(y, N, C, B). 


a 

d 


and 


y _ - (^ + 

^(1 + t ) 


> 0 , 


Proof: The existence of equilibrium Pi or P2 is obvious. We prove the existence of P3 
as follows. The equilibrium point P3(y, N, C, E) is given as the solution of the following 
system of equations, 


^Qo + IN 
So 



|s — 5 + Si 


Qo + lN \ 
So j ’ 


(2.7) 


and 



a a 


■ (2.8) 


p Y^-m y+(i/+a+ci)y+A- Is - 5 + ^ ^''^ 1 y-A-{s-6+Si '^° ^ ‘^ }N = 0 . 

5 [ Oq J 5 do 

(2.9) 

Clearly, in the N-Y plane (2.8) is a straight line and (2.9) is a hyperbola, the positive 
branch of which lies in first, second and third quadrants. 

Thus plotting the equations (2.8) and (2.9) in first quadrant (fig. 2.1), we get a branch 
of hyperbola passing through origin and intersecting the line (2.8), which gives Y and N. 
Then C and E can be found by (2.7) and we get the third equilibrium point P3. 
Remark 1: From (2.9), the slope is given by 


(dY\ Y{0Y + XC + ^iN-Y)} 
[dNJ {PY^ + XCN) 


> 0, for F > 0, > 0. 


Also 



^{s-S + Si 
r/ + a + d + ^ {s - 5 + Si 


> 0 . 


It is noted that increases as Qo increases for s — (5 + Si^ > 0. 

Remark 2: When A = 0, i.e. in the absence of carrier population for existence of an 

M 

endemic equilibrium, i.e. for Y to be positive, we must have a threshold as > 1 

(from (2.9)), which is same as mentioned in Gao and Hethcote (1992). 
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Thus the condition for the existence of second equilibrium point (Y,N), is proved. 

Remark 1. Prom (2.13), (div)(oo) ” i^+a+d+ACm ^ (2.15) 

which increases as Cm increases. 

Remark 2: It is also noted from (2.13) that > 0 for iV > y, which implies that 
equilibrium infective density increases as Cm increases. 

Remark 3: It is noted that if ^ = (1 — a)-^, then (2.13) gives a parabola and in this 
case also there exists a unique positive root N in (0, K). 

2. 3. 1.1 Stability Analysis 

Now we present the stability analysis of these equilibria. The local stability results are 
stated in the following theorem. 


THEOREM 2.6 The equilibrium jE'i(0, 0) is unstable and the equilibrium E 2 {Y,N) is 
locally asymptotically stable provided 

(by + ^{2N -K) + a{B -{I -a) ~}Y + AC^c. > 0. 


Proof: The variational matrix Mi at £?i(0,0) corresponding to the system of equations 
(2.11) is given by 

/ — A Cm — {jy -\r oc d) A Cj, 


Ml 


-a 


)■ 


Since one eigenvalue of Mi is positive, Ei is unstable. 

The variational matrix M at E 2 (Y, N) corresponding to system of equations (2.11) is 
given by 


=( 


/jy + AC„-(l-a)r 
2N) 


z. 

K 


-a 


The characteristic polynomial is given by 

t/," + {;8 F + + L(2JV - K)}^ + + ^^) ^(2JV - K) 



Infective population 2^ l3? Infective population 
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gure 2.15: Variation of infective population with time for different carrying capacities 
carrier population. 



Figure 2.16: Variation of infective population with time for different intrinsic growth 
rates of carrier population. 
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Figure 3.11: Variation of infective population with time for different growth rate of 
bacteria population due to infective human population. 



Figure 3.12: Variation of infective population with time for different growth rate of 
bacteria population corresponding to environmental discharges. 
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Susceptible population 


Figure 3.18: Variation of infective population with susceptible population, 
where pi > max | and P2 = ■ 

As the system (3.18) is bounded by its corresponding system with o; = 0, which is globally 
stable, one may speculate that the equilibrium P 4 may be globally stable. The system 
(3.18) is integrated by considering the following set of parameters, which satisfy the local 
stability condition. 

/? = 0.00000031, A = 0.00000000021, z/ = 0.012, a = 0.0005, 

d = 0.0004, a = 0.3, r = 0.0003, K = 40000, s = 1 , L = 5000000, 

So = 0.65, Si = 10, 5 = 0.000002, Qo = 20, So = 0.001. 

The equilibrium values for this set of parameters are obtained as 

F = 3358.354, V = 33270.632, 5 = 2032611.913. 

Simulation is performed for different initial positions 1 , 2 , 3 and 4 as shown in Fig. 3.18. 
From this figure, it is clear that this equilibrium may be globally stable provided that we 
start away from the other equilibria. Also in Figs. 3.19-3.23, the variation of infective 
population is shown for different s, Si, 5, L and A respectively. It is concluded that with 
the increase of these parameters, the infective population increases. 
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Figure 3.26: Variation of infective population with susceptible population. 

same set of parameters as given in Case I, with Qq = Qa, and an additional parameter 
I = 0.000005, which satisfies the local stability conditions. 

The equilibrium values for this set of parameters are determined as: 

Y = 3359.500, N = 33267.756, B = 2034237.432, E = 20166.338. 

Simulation is performed for different initial positions 1, 2, 3 and 4 shown in Fig. 3.26. 
From this figure, it is clear that this equilibrium is globally stable provided that we start 
away from the other equilibria. Also effects of various parameters such as s, Si, d, L, I 
and r, on the infective population are shown in Figs. 3.27-3.32. As before it is concluded 
that due to increase in any of these parameters, the infective population increases as 
expected. 
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Figure 4.6: Variation of infective population Yj with time for different rates of cumulative 
environmental discharges. 



Figure 4.7; Variation of infective population Yi with time for different 1 . 
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Figure 4.11; Existence of equilibrium point. 


where Rqq is the threshold of the system (4.18). Hence (4.20) intersects the Ni axis at 
least once. 

(iii) For Yi = 0, we get the following quadratic 


P1P2 


^2 


- (1 - a)(l - a')-^| ^^ 1 + 0102 — {r 2- 02 + ^2 

Ki ) do L ^2 I Oq J 


— (1 — Q')'^ |o !2 + 0^2 + (1 — Cl'){r 2 — Q ;2 + 52^)| ~ (1 ~ + Q!i + 

— (t'l + CKi + di) ■|q! 2 + d2 + (1 “ <^0(^2 — 0:2 + *^2^“)|' = 0. 


iVi 


(4.22) 


It is noted that under condition (4.21), the coefficient of is positive and the constant 
term in (4.22) is negative which implies it has one positive and one negative root. Also, if 
> 0 is positive, then from Fig. 4.11, it is clear that there exists a unique intersection 
point (Fi, Ni) satisfying both (4.19) and (4.20) which corresponds to the fifth equilibrium 
point Es under condition (4.21) and > 0. is completely determined using (Al) 
and (A2). 
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whorO TTlii (^1^ "I” "i" Q!i cJi 4 ' 5 i), 77144 — — [^ 92 ^ 4 ” 4 " ^>12 4 “ <^2 4 “ (1 — 

rriis = ^ 2^1 4- A 2 Z 1 - (1 - a')j^Y 2 and 97755 = r 2 - 0:2 4- 52 E - ^iVa- 
The variational matrix Mi at Pi is given by 



' 4- Oil 4- di 4- <5i) 

0 

0 

§iA 

dl 

0 

0 ^ 




-di 

0 

0 

0 

0 


Ml = 

-Oil 

0 

0 

0 

-di 

0 

0 

— (oi 2 4 - (^ 2 ) 

0 

0 

0 

0 

) 


0 

0 

0 

0 

r2 — Q!2 + S 2 E* 

0 



^ 0 

0 

1 

0 

0 

-SqJ 


Clearly one root is positive (as r 2 > 0 : 2 ) implying instability of this equilibrium point 

The variational matrix M 2 at equilibrium point P 2 is 




f — (j^i 4- Oil 4 - di 4- <5i) 

0 

0 

iiA 

d\ 

0 

0 \ 




-di 

0 

0 

0 

0 


Ms = 

-Oil 

^2iV| 

0 ■ 
X2N^ 

-dl 

0 

0 

77744 

0 

0 

0 

0 

? 


0 

0 

0 

0 -(r2 

— Oi2 + 52 E*) 

S2N^ 



0 

0 

1 

0 

0 

—do J 



where 77744 = -{ 0:2 4- ^2 4- (1 — 

The characteristic equation corresponding to matrix M 2 is given by 


(di 4--0)(r2 - a2 + 62 E* 4- V')(<5o 4- + 52‘0 4- £^ 3 } = 0, 


where 


g-^ z= l/i + 0^1 + 2di 4 - + 0^2 + (1 ” ^ 

g2 = {I'l + Oil + di + 61) + |o ;2 + cJ 2 + (1 I 

+di |a2 + d2 + (1 - 

53 = di |q'2 4- ^2 4- (1 — ®0^-A^2 I + CKl 4- di + 5l) - PiAN2{P2 4- ^2^) • 

Clearly three roots of the above polynomial are negative, while other roots are given by 
the cubic equation. For no positive root to exist, we must have positive, which gives 
R'g <1. Further, 


9i92 — 9z = 


4- Oil 4" 4" i5i 4“ ■|^Qi2 4” ^2 4* (1 0, )^~‘^2 

X ^di 


4- Oil 4“ di 4- 4" {o!2 4- ^2 4- (1 — 


} 


+ |q; 2 4- ^2 4- (1 — a0^^2 1 4- 0:1 4- di 4- 5i) — Pi^2-^J^^ 
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Clearly in the N-Y plane (5.11) is a parabola with vertex at ^4^) passing 

through (0, 0) and (ilT, 0). 

From (5.15), the following points are noted. 

(1^1 + CKi + + (Ji) /q !2 + ^2 + (1 - cl) ^N 2 \ 

(i)ForiVi = 0, Yl = p-_ < 0. 

{P 2 + -^2^) |/0i(l + '^)-^2 + i-'i + ai 4- di + (5i > 


(ii) For Ni = Ku 


AA*^2/^l-{a2+d2 + (l-, 


{i^x + ax + dx + A + (1 - 

n 

A* 

A h 1 ' 

1 iV2 + {vi + ax + dx + A + (1 CL)r x} 

^ u, 

r\ ^dx + (l-a)rxj 


provided 


A {^2 + ^2 




dx + (1 — cl)ti 


}N2Kr 


|q:2 + a + (1 ~ + CKl + di + A + (1 ~ 0)^1} 

where A* = {A + 


= Rs > 1, 


(5.16a) 


(iii) Also Yi = 0, gives the following quadratic in A'l, 


where 


/i = (1 - a) 


rj_ 

Ki 


/liVi +/2iVl+/3 = 0, 


P 1 P 2 N 2 — ■! a2 + A + (1 “ <20^-^2| (1 ~ 


> 0 . 


(5.166) 


(5.16c) 


The condition (5.16c) is true as we assume that (5.16a) holds for di = 0. 

/2 = A (Adi + A2A)-^2 ~ |q; 2 + d2 + (1 — {c'l + ai + 2di 4- A)(l — 


and 


/s = — |q; 2 + d2 + (1 — + 0:1 + di + A)di < 0. 

Clearly the quadratic (5.16b) has one negative root and one positive root, say N*q. 

(iv) The slope of (5.15) at {Niq, 0) is 


UnJ 


W 2 N 2 - {^2 + d2 + (1 - a')|^-^2}(l - a)^] D*^ + AA2AdiJV2 


n X t ->*2 


'K 2 


Ki 


D- + A2A] 


A (£>• + A)^2 + + ax + di + A + (1 - a)^^NiQ} D* 


> 0 
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By the Routh-Hurwitz criteria the equilibrium E4 is locally asymptotically stable if 02, 
fls positive and also aia2 — 03 > 0, which are satisfied if Ms < 1, otherwise if .Rj > 1 
this is unstable and the fifth equilibrium exists. 

Now the variational matrix M5 corresponding to the fifth equilibrium ^5(^1, Zi, Ni, Y2, N2, E) 
is given by 


57X11 

-P 1 Y 2 

mA3 

m'u 

0 

0 \ 

5i 

m22 

'“(1 ~ 

0 

0 

0 

-O'! 

0 

mzz 

0 

0 

0 

P2{N2 - Y 2 ) 

HN2 - Y 2 ) 

0 

mu 

m45 

0^ 

0 

0 

0 

0 

77 x 55 

S2N2 

0 

0 

1 

0 

0 

— 5o / 


where rriu — —{^{¥2 + + Oi + + (5i + (1 — a)j^Ni), 77x22 = + (1 “ 

rriiz = P1Y2 - (1 - a)j^Yr, mu = A(iVi - - Z,), m'ss = - K,), 

rriii = —{^2^1 + ^2-^1 + a2 + (^2 + (1 — 

^45 = ^2^1 + A2.^i — (1 — a')j^Y2 and 77x55 = 

The characteristic polynomial corresponding to the above matrix is given by 

7/;® + + a4'0^ + + a2'0^ + aixp + Uq = 0, 

where 

<^5 — ^0 ~ (^11 + ^22 + 77x33 + ^44 H" ^55) > 0» 

CX4 = ~5o(^ii + ^22 + ^33 + ^44 + ^^55) + ^11 (^22 + ^33 + ^44 + m55) 

4-m'22 (m33 + niu + m55) + rrizz {mu + m'ss) + ruu rrizz + Q;imi3 
+<^1/^2^ — P 2 mi 4 {N 2 — F2)) 

03 = Soniii {rn 22 + m''33 + m '44 + m' 55 ) - m '22 {rrizs + mu + m^z) 

—77x22 mzz (—<^0 + ^44 + ^55) “■ ^^33 mu (~^o + + m^b) 

—m'44 m 5 s{— 5 o + mu + 77x22) + 77x5550(77x22 + mzs + mu) 

+5om'22 m44 + PiY 20 ii{l — a)—Zi + (5o ~ ^33 “ ’^'^44 ~ mzz)5iPiY2 

El 

+a:i7rii3(5o — m22 — mu — m^z) — m^i 77x3377x55 

—^2mu{N2 - Y 2 ){ 5 o - m22 - m'33 - m'55) - 5iA2771i45o(iV2 - F2)) 

U 2 = — 5o^ll m 22 ( 57 X 33 + 55 X 44 + 57 X 55 ) + TTlii 77X22 57 X 33 ( 77 X 44 + 77 X 55 ) + 77122 77133 77 X 44 
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012 ^ 2 ) (^2 + 0)[P2y2 "I" A 2 CT 2 + ( 2^2 + 0:2 + (^2 + ^)^]- (6.16) 


A. When C=0: 

It is obvious that the following three equilibria exist. 


(i) When C — 0 and 12—0 , we have two equilibrium points namely, 


Aiid2 + 0}-\-OA2 ^ ^2 ^ Qo + l^e\ ^ (,,, ,r* . A2 „ Qo + l^' 






where V* = + (^2 + 0 ){/ 3 iAi - di(i/i + Q!i +cgi)} ^ + z^i + ni + di 

A(^2 + 0)(c^l + ctl) ’ ^ /3i 

E 2 exists if iSi( 0A2 + (d 2 + 0)Ai) > di(d 2 + + ai -h di). 


(ii) When C = 0 , )<2 = V2*, 


we get tlie third equilibrium point as E3 ^Y^** , N** , ¥2** , N2* 
following condition 


Qo + IN2* \ 

5o J 


under the 


P2A2 

d2 ~¥ 0 


> {V2 4- 0:2 + ^2 + 9), 


where Y" = , + + ^1, iv," = 

^2 (^2 + ^ + 0:2) d2-\- 9 

Yj** is the positive root of the quadratic (6.15) and then Nl* and E** are as follows: 


V** = (^2 + ^)(^4i - ttiY;*) + 9{A2 - 0 : 2 ^ 2 **) ^ Qo + lN^\ 

di (<^2 4 " 9) dt) 


B. When C # 0, 

A2 . . 

In this case, we have C = -^(i^ — GY2), which is positive as 12 < ~ for positive N2, 


(using (6.12)), 

, r, Qo S1IA2 , ^ ^ 1^0:2 

where F — s — So 4" ■Sj-jr — h . / , ; and G x f/i ^ o\' 

So So{d2 + 9) 6o{d2 + 9) 



Chapter 7 


Modelling the Spread of Bacterial 
Disease in a Population: Effect of 
Service Providers from an 
Environmentally Degraded Region 


7.1 Introduction 


Infectious diseases are the world’s biggest killer of people and account for more than a 
dozen million deaths per year. The situation has become worse by uncontrolled increase 
in population due to a high birth rate as well as migration/movement of populations over 
the past several decades. Further the poor environmental conditions already existing in 
the densely populated cities of the third world countries have the greatest impact on the 
spread of bacterial disease such as tuberculosis, typhoid and so on. If the environment is 
conducive to the growth of the bacteria population, then this further helps in the spread 
of infectious diseases (Cairnross and Feachem 1983, Harold 1960, Taylor and Knowelden 
1964). 

In Chapter 6, we have modeled and analyzed the effect of immigration of a population 
from an environmentally degraded region on the spread of an infectious disease in a 
population living in a much cleaner environment. It may be noted here that when rich 



7.2 SIS Model with Immigration 
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Figure 7.3: Variation of N 2 with Y 2 . 



Figure 7.4: Variation of Yi and Y 2 with time for different intrinsic growth rates of the 
bacterial population. 




